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Time—  resolved  x—  ray  scattering  has  been  used  to  study  ordering  kinetics  in 
single  crystal  bulk  Cu«Au,  as  well  as  in  sputtered  and  molecular  beam  epitaxy 
grown  films.   After  annealing  at  high  temperatures  the  sample  is  rapidly  quenched 
to  fixed  temperatures  below  the  order—  disorder  transition  temperature.   The 
development  of  order  is  monitored  in  real  time  using  scattering  techniques. 

The  bulk  sample  clearly  showed  three  regimes:  nucleation,  ordering,  and 
coarsening.   The  anisotropic  superlattice  peaks  that  reflect  the  domains  structure 
are  investigated  in  connection  with  the  ordering  kinetics.   The  line  shape  of  the 
scattering  function  exhibits  a  crossover  from  gaussian  to  lorentzian—  squared  as  the 
system  goes  from  the  ordering  regime  to  the  coarsening  regime.   Coarsening  in 
Cu«Au  is  consistent  with  curvature  driven  growth. 

Domain  coarsening  in  stoichiometric  sputtered  films  is  also  consistent  with 
curvature  driven  growth.   However,  coarsening  in  copper  rich  films  proceeds  much 


more  slowly.   The  results  suggest  the  extra  copper  affects  the  ordering  kinetics  in 
the  same  way  diffusive  impurities  would,  resulting  in  a  logarithmic  like  time 
dependence. 

The  M.B.E.  films  show  a  slowing  of  the  growth  at  late  times.   The  450oX 
film  starts  out  with  curvature  driven  growth  but  then  continuously  slows  down  as 
the  domains  grow.   The  710a  film  shows  an  interesting  temperature  dependence  for 
the  growth,  in  such  a  way  that  at  temperatures  close  to  the  transition,  the  domain 
growth  almost  freezes  at  late  times.   The  dominate  factor  is  probably  strain,  all  of 
the  trends  for  slower  growth  are  consistent  with  greater  strain.   The  dimensionality 
in  the  M.B.E.  film  systems  is  considered.   The  scaling  in  the  4500X  and  710A  films 
is  clearly  three  dimensional.   However,  the  dimension  of  the  scaling  in  the  260X  film 
is  unclear. 
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CHAPTER  1 
BACKGROUND 

Introduction 


The  kinetics  of  phase  transitions  in  systems  that  exhibit  ordering  after  a 
quench  from  above  a  critical  ordering  temperature  T    has  received  much  recent 
interest.^       ''   One  focus  of  this  interest  has  been  on  whether  universality  classes  for 
the  kinetics  of  these  systems  exist  as  they  do  for  ordinary  static  critical 
phenomena^   '.   Conservation  laws  are  known  to  greatly  affect  the  growth  kinetics. 
However,  the  effects  of  symmetry,  structure  and  sharpness  of  domain  walls,  order  of 
the  degeneracy  of  the  ground  state,  impurities,  and  stoichiometry  are  not  well 
understood.   Examples  of  kinetic  behavior  can  be  found  in  many  systems  including 
magnetic  materials,  copolymers,  and  alloys. 

This  dissertation  focuses  on  a  careful  study  of  the  ordering  kinetics  in  the 
alloy  Cu„Au.   Cu^Au  is  a  prototypical  system  with  a  first  order  transition,  and  a 
nontrivial  asymmetry.   The  elements  Cu  and  Au  have  considerably  different  atomic 
numbers  (Z),  making  x— ray  measurements  of  superlattice  peaks  easier.   Moreover, 
the  equilibrium  properties  of  Cu^Au  are  well  studied^  ^   Two  important  properties 
of  Cu«  Au  are  its  four— fold  degenerate  ground  state  and  the  existence  of  two  types  of 
domain  walls  that  separate  the  four  ordered  phases.   These  two  differences  make 
Cu^Au  a  good  choice  for  furthering  our  understanding  of  the  non— equilibrium 
kinetics  of  first  order  phase  transitions.   The  non— equilibrium  kinetics  are  less  well 
understood,  particularly  in  the  context  of  the  ideas  of  modern  statistical  mechanics. 


Interesting  effects  in  the  non— equilibrium  kinetics  are  introduced  by  the  existence  of 
different  domain  wall  structures  and  a  corresponding  anisotropy  of  the  correlations. 

Alloys  like  Cu^Au  exhibit  a  change  in  their  atomic  arrangement  as  a 
function  of  temperature.   Such  alloys  have  markedly  different  equilibrium  properties 
as  a  function  of  temperature.   However,  the  present  interest  is  focused  on  the 
approach  to  equilibrium. 

The  crystalline  form  of  Cu„Au  is  face  centered  cubic  (fee),  and  below  the 
critical  temperature  T    (T  ~  390  C)  has  an  equilibrium  state  where  the  Au  atoms 
preferentially  reside  at  cube  corner  sites  and  the  Cu  at  face  center  sites,  (see 
figure— 1).   Above  T    the  atoms  are  randomly  distributed,  but  short  range  order  stiU 
persists.   The  transition  from  the  ordered  state  to  the  disordered  state  is  first  order, 
and  like  so  many  things  in  nature  the  change  from  order  to  disorder  takes  place 
much  faster  than  the  change  back  to  order. 

The  experiments  discussed  here  are  all  concerned  with  the  change  from 
disordered  to  ordered  states.   In  a  typical  experiment  the  sample  is  held  at  a 
temperature  well  above  T   in  order  to  insure  that  the  sample  is  in  a  disordered 
state,  from  which  the  temperature  is  rapidly  lowered  to  a  point  below  T    and  held 
there  for  several  days.  X— ray  scattering  scans  are  collected  as  a  function  of  time. 
The  peak  grows  with  time,  becoming  narrower  and  more  intense.   Moreover,  the 
shape  of  this  x-ray  peak  is  important.   Many  theories  predict  this  shape  and  hence, 
a  careful  investigation  is  important  to  understanding  the  underlying  physics. 
Surprisingly,  there  has  been  little  experimental  work  in  this  direction. 
Unfortunately,  it  is  usually  difficult  to  obtain  the  exact  shape  due  to  effects  from 
instrumental  resolution,  finite  size,  strain,  and  mosaic  spread. 

Experiments  have  been  performed  on  a  bulk  single  crystal,  sputtered  films, 
including  one  Cu  rich  film,  M.B.E.  (molecular  beam  epitaxy)  grown  films  from  260 


Cu^Au 
ordered   cube 


FIG.  1.  Schematic  of  (FCC)  structure  for  ordered  Cu^Au. 

of  a  domain  increases.   Larger  domains  grow  at  the 
expense  of  their  smaller  neighbors. 


to  4500  angstroms  thick.   For  each  sample  the  experiment  was  repeated  for  a  series 
of  temperatures. 

From  the  experiments  on  the  bulk  single  crystal  it  is  evident  that  three 
distinct  mechanisms  are  involved,  nucleation,  ordering,  and  coarsening  as  defined 
below.   Nucleation,  as  the  name  suggests,  refers  to  a  delay  between  the  quench  and 
the  start  of  the  development  of  long  range  order.   Although  in  several  quenches  a 
nucleation  period  is  observed,  the  time  resolution  in  these  experiments  was  such 
that  this  regime  was  not  well  characterized.   Others  using  a  different  experimental 
set  up  and  a  more  powerful  x— ray  source  have  probed  this  aspect  of  the  transition  in 
more  details  ''  Nucleation  refers  to  a  delay  between  the  quench  and  the  start  of  the 
development  of  long  range  order.   Ordering  (or  growth)  is  the  process  by  which  the 
disordered  state  condenses  onto  the  various  nucleated  droplets  of  the  ordered  phase. 
One  of  our  results  is  that  the  ordering  process  can  be  characterized  by  a  change  in 
the  functional  form  of  the  scattering  associated  with  the  structure  factor.   The 
movement  of  domain  boundaries  that  separate  the  different  ground  states  is  known 
as  coarsening.  This  name  refers  to  the  fact  that  as  time  progresses  the  average  size 
of  the  domains  grow  as  the  larger  domains  convert  their  smaller  neighbors  to  the 
same  phase. 

Another  sample  is  an  off-stoichiometric  sputtered  film  of  Cu«Au.   It  is 
interesting  that  an  off— stoichiometric  sample  exhibits  a  different  type  of  growth 
than  its  stoichiometric  counterpart.   Moreover,  the  off— stoichiometric  film  is 
probably  in  a  different  class  altogether.   The  extra  Cu  in  the  sputtered  film  acts  as  a 
diffusive  impurity.   However,  the  stoichiometric  sputtered  films  (all  over  7000 
angstroms)  exhibited  curvature  driven  domain  coarsening,  where  the  size  of  an 
ordered  domain  is  proportional  to  the  square  root  of  time.   The  differences  in  the 
time  dependence  usually  involve  a  lowering  of  the  power  of  the  time  dependence 


from  1/2  (square  root)  to  a  value  less  than  1/2.   However,  the  growth  in  the 
off— stoichiometric  film  is  also  consistent  with  a  logarithmic  time  dependence. 

The  last  samples  are  single  crystal  films  grown  by  molecular  beam  epitaxy 
(M.B.E.),  which  are  reviewed  in  a  contexts  similar  to  the  off-stoichiometric 
sputtered  film.   Both  the  off-stoichiometric  sputtered  film  and  the  M.B.E.  films 
display  a  slower  type  of  time  dependence  for  growth.   However,  the  cause  of  this 
change  from  the  bulk  behavior  at  critical  composition  arises  from  slightly  different 
causes.   Consequently,  similarities  exist,  but  so  do  differences. 

Previous  Work 

The  earliest  studies  of  the  order  disorder  transition  in  Cu^Au  date  at  least  as 
far  back  as  the  early  1900's.''  ''   However,  the  interest  then  was  how  the  level  of 
order  affects  bulk  properties.   The  present  interest  is  in  the  ordering  itself  and  the 
non— equilibrium  properties  associated  with  it.   At  the  time  of  the  earlier  work  the 
intensity  of  available  x— ray  sources  did  not  permit  the  in— situ  type  of  measurement 
made  in  this  work.  The  early  experiments  measured  only  bulk  effects  like  resistance 
as  a  function  of  temperature,  time,  or  order.  The  only  x-ray  work  that  could  be 
done  was  static  measurements.   An  extensive  review  of  order— disorder  transitions  in 
binary  alloys  prior  to  the  late  1930s  can  be  found  in  an  article  by  F.  C.  Nix  and  W. 
Shockely.^     ''   Two  experimental  papers  both  from  1936  are  worth  noting,  one  by 
Sykes  and  Evans^     ''  and  the  second  by  Sykes  and  Jones. ^  *"''   These  two  works 
describe  the  electrical  resistance,  specific  heat,  and  the  crystal  structure  as  related 
to  the  order— disorder  transition.   A  later  paper  by  Jones  and  Sykes^     -'  along  with 
their  earlier  work  are  the  earliest  references  relating  the  width  of  the  superlattice 
peak  to  domain  size. 


Bragg  and  Williams^    ^^^  in  1934  and  35  put  forth  the  first  equilibrium 
theory  from  which  the  qualitative  behavior  of  bulk  properties  such  as  the  specific 
heat,  the  level  of  long  range  order,  and  resistance  could  be  derived.   Formally  this 
theory  is  analogous  to  the  Weiss  theory  for  ferromagnetismA  ■'  This  theory  worked 
only  below  the  transition  temperature  T  .  The  Bragg— Williams  theory  is  a  static 
mean  field  approximation  that  does  not  consider  fluctuations.   The  Bragg— Williams 
theory  considers  the  long  range  interactions  between  atoms  and  therefore,  does  not 

explain  short— range  order  and  only  qualitatively  describes  the  equilibrium 

fl7l 
properties  near  T  .   The  theory  of  Bethe*-     '  considers  nearest  neighbors.   As  in  the 

Bragg-Williams  case  the  energy  just  below  T    is  the  same  as  just  above  T  . 

However,  in  Bethe's  theory  additional  energy  is  required  to  destroy  short  range 

order  above  T  .   This  gives  rise  to  an  anomaly  in  the  specific  heat  in  agreement 

with  experiment.   Due  to  the  complicated  nearest  neighbor  configuration  and 

asymmetry  in  composition  found  in  A„B  systems  like  Cu^Au,  many  of  Bethe's 

assumptions  regarding  the  structure  do  not  apply  to  our  case.   Fortunately 

Peierls'^     '  has  applied  Bethe's  theory  to  this  case.   A  more  formal  approach  to  the 

theory  of  order-disorder  kinetics  appears  in  an  article  from  1956  by  G. 

(19) 
Vineyard'^     ' .   His  theory  starts  with  a  set  of  multi  particle  distribution  functions 

and  equations  of  motion.   The  simplest  approximation  for  these  equations  yield  the 

Bragg-Williams  result.   The  Bethe— Peierls  approximation  can  also  be  derived  from 

this  theory  if  only  nearest  neighbor  interactions  are  considered. 

In  the  70s  interest  in  the  kinetics  of  Cu«Au  was  renewed  by  the  development 

of  theories  that  considered  the  movement  of  domain  walls^  '  -'  as  well  as  the 

availability  of  more  intense  x-ray  sources.   Allen  and  Cahn^     ''  derived  the  time 

dependence  of  the  size  of  an  ordered  domain.   Their  theory  is  often  referred  to  as 

curvature  driven  (or  controlled)  growth  because  the  driving  force  goes  as  the 


curvature  of  the  wall.   The  final  result  is  that  the  size  of  an  ordered  domain  L(t)  is 
proportional  to  t  '     (the  square  root  of  time).   Earlier  theories*-  '  predicted  that  the 
velocity  of  the  wall  would  go  as  the  product  of  the  mean  curvature  and  the  interface 
free  energy.    Therefore,  the  earlier  versions  predicted  critical  slowing  down  of  the 
wall  as  T   is  approached  from  below.   In  the  Cahn— Allen  case  the  diffusional 
dissipation  of  free  energy  of  a  domain  wall  exactly  equals  the  available  driving  force. 
Hence,  the  velocity  of  a  domain  wall  is  not  a  function  of  the  driving  force,''     '  and 
the  power  law  that  describes  the  time  dependence  is  not  a  function  of  temperature. 
The  temperature  dependence  is  only  in  a  proportionality  term  that  precedes  the 
power  law  time  dependence.   This  renewed  interest  is  actually  part  of  a  growing 
interest  in  universality  classes  and  non— equilibrium  phenomena  which  occur  in  such 
transitions. 

The  kinetics  of  transitions  have  been  of  great  interest  in  the  past  decade. 
Theoretical,^ ^'^'^^""^^^  experimental, ('^°~^'^)  and  computer  simulation^'^''^^'^^"''^) 
studies  on  first  order  phase  transitions  abound.   A  variety  of  systems,^     '  binary 
alloys,  binary  fluids,  magnetic  materials, block  copolymers,  glasses,  ferroelectric 
crystals,  and  others  can  be  used  to  study  this  phenomenon.   The  problem  is  often 
discussed  in  terms  of  an  Ising  model, ^     ^  partly  because  the  order  parameter  is  often 

simply  the  magnetization,  an  easily  measured  quantity,  and  partly  because  Ising 

(3) 
like  systems  are  easily  classified.   Ferromagnets  with  spin  exchange  (Kawasaki)^  ' 

dynamics  are  systems  with  a  conserved  order  parameter  (COP),  ferromagnets  with 

spin  flip  (Glauber)^  ^  dynamics  have  a  non-conserved  order  parameter  (NCOP), 

and  antiferromagnets  with  spin  exchange  dynamics^     '  have  a  NCOP  (sub  lattice 

magnetization)  even  though  the  net  magnetization  is  conserved.   The  reason  being 

easily  classified  is  so  important  is  that  the  nature  of  universality  classes^     ^  for 

growth  is  not  fully  understood. 


Binary  alloys  are  also  useful  systems  for  kinetic  studies  for  similar  reasons. 
One  obvious  division  in  classes  is  whether  the  order  parameter  is  conserved  or 
non— conserved.   A  recent  computer  simulation*-     '  on  a  three— dimensional  binary 
alloy  model  (of  the  AB  type)  with  a  NCOP  shows  qualitatively  different  ordering 

kinetics  than  the  equivalent  spinodal  decomposition  found  in  binary  alloys  with  a 

1  /2 
COP.   This  simulation  possesses  a  characteristic  length  that  increases  as  t  ',    and  a 

corresponding  scaling  of  the  structure  factor  of  L  ,  where  d  is  the  dimension. 

Phase  separation  is  an  example  of  a  process  with  a  COP  that  has  been 
extensively  studied  experimentally^     '"'''.   Phase  separation  exhibits  diffusion 
controlled  grouth*-     ''  where  L(t)  is  proportional  to  t  '    as  compared  to  the  NCOP 
process  of  curvature  controlled  grouth^     '     -'  where  L(t)  is  proportional  to  t  '   . 
There  are  also  systems  that  exhibit  interface  controlled  growth^  ''  where  L(t)  is 
proportional  to  time.  Hence,  the  power  of  the  time  dependence  seems  to  be  an 
important  difference  between  classes^     ^   Another  class  may  be  found  in  systems 
with  impurities''  -'  where  L(t)  is  proportional  to  ln(t)"^.   It  has  been  suggested 

this  class  can  be  further  divided  by  the  value  of  the  power  m^     '',  possibly 
indicating  whether  the  impurities  are  quenched*-  ''  or  diffusive*-     '     ^   This 

type  of  growth  might  also  be  found  for  critical  quenches  (close  to  T  )  involving 
spinodal  decomposition*-     '. 

One  purpose  of  this  work  is  to  distinguish  different  universality  classes, 
specifically  the  effect  of  dimensionality,  composition,  the  order  of  the  ground  state 
degeneracy,  and  the  nature  of  the  scaling  form  for  the  structure 
factor*-  '  '    ~    '     /.   There  is  evidence*-     '  that  the  mechanism  of  domain  growth 
will  change  for  thin  films.   Also,  there  was  some  question  as  to  whether  the  order  of 
the  degeneracy  lowers  the  power  of  the  time  dependence*-     ''''•'. 
Another  reason  to  look  at  Cu^Au  is  that  it  has  two  types  of  domain  walls*-     '^^'^^J. 


The  first  of  these  walls  conserves  the  nearest  neighbor  configuration  and  due  to  the 
energies  involved*'  ^  is  easily  formed  compared  to  the  second.   A  comparison  of  the 
growth  of  these  two  very  different  wails  should  provide  further  insight  into  the 
study  of  order  disorder  kineticsA     ■'   Moreover,  less  work  has  been  done  on  systems 

involving  a  NCOP^''''*^'^'^'^'^^  than  COP.   After  a  thermal  quench  both  types  of 

(S  94  95) 
systems  exhibit  three  distinct  regimes^  '     '     '  nucleation,  ordering,  and 

coarsening.   The  functional  form  of  the  structure  factor  can  be  derived  from  most 

applicable  theories,^  .     >     i     i     ;  however,  instrumental  resolution,  finite  size,  and 

strain  often  have  sizable  contributions  to  the  shape  of  the  scattering,  making  it 

difficult  to  discern  the  exact  shape  of  the  structure  factor  experimentally. 

Determining  this  shape  is  another  purpose  of  this  work. 

The  rest  of  this  dissertation  is  arranged  as  follows:  Chapter  2  discusses  the 

theories  of  domain  coarsening  and  nonequilibrium  scaling.   The  first  part  of  chapter 

3  is  on  basic  x— ray  diffraction  and  how  the  level  of  order  is  related  to  the  diffraction 

pattern.   A  description  of  the  different  types  of  domain  walls  in  Cu^Au  can  be  found 

at  the  end  of  chapter  3.  The  experimental  analysis  and  results  are  discussed  in 

chapter  4.   Conclusions  are  presented  in  chapter  5. 


CHAPTER  2 
THE  GROWTH  OF  ORDER 


Cu«Au  is  a  prototypical  system  for  order-disorder  transitions''     K  The 
transition  is  driven  by  minimizing  the  free  energy 

F=U-TS.  (2.1) 

At  low  temperatures  the  free  energy  is  dominated  by  the  internal  energy,  U,  causing 
the  minimum  in  free  energy  to  correspond  to  the  ordered  state.   However,  at  high 
temperature  the  minimum  for  the  free  energy  corresponds  to  a  disordered  state 
because  it  is  dominated  by  the  entropy  term,  (— TS).   The  critical  temperature  is  a 
function  of  the  concentration  of  the  different  atoms.   Figure  2  is  a  generic  phase 
diagram  for  alloy  systems  with  order  disorder  transitions.   Below  the  solid 
coexistence  curve  the  equilibrium  state  corresponds  to  order  and  above  to  disorder. 
The  dashed  curve  in  figure  2  represents  the  classical  spinodal  curve  separating 
regions  of  metastability  and  instability  of  the  disordered  state  after  a  thermal 
quench  to  below  T  .   For  phase  separation  these  two  curves  meet  at  the  critical 
composition.   However,  for  the  process  of  superlattice  formation^  ''  this  is  not  the 
case,  producing  an  extended  nucleation  period  even  for  the  stoichiometric 
composition. 

For  instantaneous  quenches  to  temperatures  just  below  T    the  short— range 
order  fluctuations  associated  with  temperatures  above  T    are  metastable.   This 
metastable  state  then  decays  during  a  finite  incubation  time  towards  a  new 
equilibrium  value*"  ''.   During  the  incubation  time  the  fluctuations  are  such  that  the 
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FIG.  2.  Coexistence  curve  (solid)  and  classical  spinodal  curve  (dashed)  for  a 
first  order  transition  in  a  typical  binary  alloy.  Points  M  and  U  label 
regions  as  metastable  or  unstable. 
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critical  size  for  the  stability  of  an  ordered  droplet  is  not  realized.   The  nucleation 
stage  starts  when  the  fluctuations  have  relaxed  sufficiently  to  allow  a  stable  ordered 
droplet  to  form.  As  quenches  become  deeper  the  fluctuations  associated  with  the 
disordered  state  become  unstable  and  nucleated  ordering  changes  to  continuous  (or 
spinodal)  ordering.   There  is  some  indication  that  this  change  occurs  well  below  the 
classical  spinodal  temperature.^  ^   In  fact  there  is  no  abrupt  change  as  nucleated 
ordering  changes  to  spinodal.   After  nucleation  droplets  of  the  ordered  phase  grow  in 
a  matrix  of  the  disordered  phase  (see  fig  3a).   The  individual  droplets  may  be  in  any 
of  the  p  allowed  ground  states.  Eventually  the  domain  boundaries  meet,  and  the 
system  is  composed  of  domains  in  different  ground  states  separated  by  domain  walls 
(see  fig.  3b, c).   The  coarsening  process  then  begins  as  the  larger  domains  grow  at  the 
expense  of  their  smaller  neighbors. 

In  the  coarsening  regime  the  system  can  be  characterized  by  L(t),  the 
average  size  of  an  ordered  domain  at  time  t.  During  coarsening  it  is  expected  that 
the  system  will  be  self-similar  under  a  rescaling  of  both  space  and  time.   The 
growth  of  the  domains  is  characterized  by  the  time  dependence  of  L(t).   If  all  the 
excess  energy  is  at  the  domain  walls  the  velocity  of  the  wall  is  proportional  to  its 
curvature,  as  discussed  below.  This  leads  to  an  equation  of  motion; 

m'^^'T  (2.2) 

Separating  variables  gives 

r  dr  (X  dt  (2.3) 


Integration  shows 


which  implies 


^r    <xt  (2.4) 


L(t)  «  t^/2  (2.5) 
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FIG.  3  Ordering  and  coarsening  after  a  quench. 

alOrdering  shortly  after  the  on  set  of  critical  nucleation. 

bjEariy  stages  of  domain  coarsening.  c)Late  stages  of  domain  coarsening. 
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This  is  Cahn— Allen  (C— A)  curvature  driven  growth.   For  a  NCOP  system  with  p=2 
and  no  random  impurities  the  C— A  result  has  been  found  to  be  correct.^     '     ' 

Another  growth  mechanism  is  that  of  Lifshitz  and  Slyozov,  where  the  rate  of 
change  in  the  volume  of  a  domain  is  constant.   This  type  of  growth  is  expected  in 
systems  with  a  COP  that  exhibit  phase  separation. 

-TT-  =  constant  (2.6) 

3 
To  convert  2.6  to  a  velocity  is  a  simple  matter.   The  volume  is  proportional  to  r    so, 

dV  «  r^dr.  (2.7) 

Substitution  into  2.6  shows, 

^.V,2.  (2.8) 

The  velocity  of  a  domain  wall  is  proportional  to  the  inverse  of  the  area  of  the  wall. 
Integration  as  done  above  leads  to, 

L(t)<xt^/^  (2.9) 

a  slower  power  law. 

An  even  slower  power  law  of  t  '    has  been  recently  proposed>     ^   When  this 
new  class  was  first  proposed  it  was  explained  by  the  failure  of  the  assumptions  made 

in  the  C— A  theory  and  hence,  rightfully  came  under  attack.   It  is  only  recently  that 

1/4 
the  correct  reasons  for  this  classes  existence  have  been  revealed.   The  t  '     class  is 

associated  with  the  case  of  zero  curvature  in  the  Cahn— Allen  theory.   This  new  class 

exhibits  competition  between  sharp  and  broad  domain  walls.   It  has  also  been  found 

that  this  class  displays  a  cross  over  from  t  '    to  t  '     as  the  temperature  is 

increased. 

A  fourth  type  of  time  dependence  may  be  found  in  systems  with  energy 

barriers  to  the  growth  which  depend  on  the  characteristic  length  L,*^     '' 

E(L)  =  L^/"^/A  (2.10) 

where  A  is  weakly  temperature  dependent.   If  the  time  necessary  to  over  come  such 
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a  barrier  is  of  the  form, 

t  =  r  exp(E(L)/T)  (2.11) 

then  inverting,  and  substituting  for  E(L)  gives, 

L(t)  =  AT  [ln(t/r)]"^  (2.12) 

a  logarithmic  dependence  on  time,  where  m  is  a  constant,  that  depends  on  the  type 
of  barrier,  and  r  is  a  rate  constant. 

Obtaining  the  time  dependence  from  the  velocity  of  a  domain  wall 
(anti— phase  domain  boundary)  as  done  above  is  simple.   It  is  considerably  more 
difficult  to  show  that  the  velocity  either  goes  as  the  curvature,  or  as  the  inverse  of 
the  surface  area.   The  velocity  in  general  is  proportional  to  the  thermodynamic 
driving  force^       ■'  and  associated  constant  of  proportionality  is  called  the  mobility. 
One  way  to  consider  thermodynamic  nonequilibrium  is  through  a  mean  field 
approach.   Mean  field  theories  to  nonequilibrium  dynamics  of  metastable  states  are 
an  alternative  to  cluster  dynamic  theories,  and  are  of  practical  importance  in 
determining  the  thermodynamic  driving  forces  involved  in  domain  coarsening.^       '' 
Attention  is  focused  on  a  small  set  of  semi— macroscopic  variables  that  define  an 
order  parameter  like  the  ones  discussed  at  the  end  of  chapter  three.   These  variables 
are  considered  to  have  slow  dynamic  evolution  as  compared  to  the  microscopic 
degrees  of  freedom,  that  are  averaged,  or  enter  only  as  random  forces.   The  slowness 
is  important  because  the  derivatives  with  respect  to  time  are  frequently 
approximated  as  zero. 

The  typical  attack  is  to  use  a  Ginzburg— Landau  Hamiltonian^   '  where  the 
free  energy  density  has  a  double  well  structure  below  some  critical  temperature. 
The  double  well  corresponds  to  the  two  different  phases  ordered  or  disordered.   In 
the  COP  process  of  phase  separation  the  order  parameter  is  the  local  concentration 
of  one  species  of  atoms,  where  for  the  NCOP  process  of  superlattice  formation  the 
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order  parameter  is  the  occupation  of  a  species  of  atoms  on  a  superlattice.   This 
difference  results  in  very  different  behavior.   For  the  COP  process  the  structure 
factor  or  scattering  not  only  grows  in  amplitude  as  the  width  decreases,  but  the 
position  moves  from  high  to  low  q,  as  phase  separation  increases.   Moreover, 
because  the  order  parameter  is  conserved  the  integrated  intensity  must  be  a 
constant.   For  the  NCOP  case  intensity  at  a  superlattice  peak  is  a  function  of  the 
order  parameter  as  shown  earlier. 

The  Cahn— AUen  result  starts  by  assuming  the  free  energy  per  unit  volume  of 
a  homogeneous  phase  f   is  a  function  of  the  order  parameter  f .  The  important 
quantity  is  Af  ,  the  difference  between  a  state  of  arbitrary  order  and  an  ordered 


state.   The  equation  of  motion  is  then, 


(2.13) 

2 
where  o  is  a  positive  kinetic  coefficient  and  M=2aA;  is  a  diffusion  coefficient  (m  /s). 

Rewriting  2.13  to  correspond  to  a  direction  g  normal  to  a  surface  of  constant  f 

gives, 

di  1  (2.14) 


=  -a  ^-0  -M 


d"^^ 


where  (K, +K2)  is  the  mean  curvature  which  is  equal  to  the  negative  of  the 
divergence  of  a  unit  vector  normal  to  the  surface.  In  the  original  paper  by  Cahn 
and  Allen  g  was  defined  such  that  for  a  spherical  domain  g  was  directed  towards  the 
interior.  Although  many  reverse  this  convention  it  will  be  followed  here.   Now 
consider  a  gently  curved  domain  wall.   A  domain  wall  is  a  surface  of  constant  f . 


The  velocity  of  the  wall  is  just  5g/5t  of  a  constant  f  surface. 


V  = 


^ 


^^ 


d^ 
M 


I 


t 


but  from  2.13  for  this  surface, 


M 


g 


=  -M(K^+K2) 


51' 


t 


(2.15) 


(2.16) 
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therefore,  from  2.14  and  2.15, 

V  =  M(K^+K2).  (2.17) 

As  previously  discussed  this  result  leads  to  the  time  dependence  given  in  equation 
2.5.   Also  this  result  was  derived  in  the  contexts  of  a  second  order  transition  which 
is  being  applied  to  the  late  time  development  of  a  first  order  transition. 

The  degree  of  order  and  microscopic  structure  of  a  system  can  be 
characterized  by  the  order  parameter  correlation  function.   The  fourier  transform  of 
the  correlation  function,  S(Q),  is  called  the  structure  factor.    As  discussed  in 
chapter  3,  S(Q)  is  directly  measured  in  x-ray  experiments.   In  a  nonequilibrium 
situation,  such  as  that  following  a  quench  into  a  metastable  or  unstable  state,  the 
structure  factor  may  be  generalized  and  written  as  S(Q,t),  where  t  is  the  time  after 
a  quench. 

Self  similarity  of  the  ordering  system  under  a  rescaling  of  space  and  time  for 

a  COP  process  is  exhibited  by  the  scaling  form  of  the  structure  factor. 

S(Q,t)  =  [q(t)]^S;[q/q(t)]  (2.18) 

max  max 


For  a  NCOP  system, 

S(Q,t)  =  k(t)]^S;[(q)/.(t)]  (2.19) 

where  q=Q— G,  G  is  the  ordering  wave  vector,  d  is  the  spatial  dimensionality,  and 

e=q/q.   In  the  COP  case  q         is  a  function  of  time  and  in  the  NCOP  case  it  is  a 
^  ^  ^max 

constant.   Although  the  two  types  of  scaling  are  different  they  can  be  made  to  look 
similar.   In  the  COP  case  the  scattering  arises  from  the  domains  forming  an 
imperfect  lattice.   Hence,  the  position  [q        (t)  «  1/L(t)]  is  a  function  of  the 
distance  between  domains,  which  on  average  is  the  size  of  the  domains.   In  contrast 
the  length  scale  in  the  NOCP  case  is  inversely  proportional  to  the  peak  width,  a{t). 
Substitution  of  L(t)  into  2.18  and  2.19  gives  the  identical  result. 

S(Q,t)  =  [L(t)]^S:[qL(t)]  (2.20) 
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Scattering  is  defined  by  the  structure  factor,  therefore,  this  type  of  scaling 
implies  that  the  functional  form  of  the  scattering  remains  constant.   If  scaling  is 
valid,  only  the  amplitude  and  width  change  for  the  NCOP  system.   Equations 
2.18—2.20  also  describes  how  the  dimensionality  of  a  system  can  be  detected.   In 
general  the  scaling  function  S  is  anisotropic,  reflecting  the  structure  of  both  the 
domains  and  their  walls. ^     ■' 

A  recent  proposal''     ''  divides  growth  kinetics  into  four  classes  characterized 
by  different  low  temperature  behaviors.   For  class  1  systems  the  temperature  is  not 
important,  such  as  a  simple  Ising  ferromagnet  with  Glauber  dynamics.   Class  1 
systems  obey  power  law  domafn  growth  at  all  temperatures.   Cu^Au  probably 
belongs  to  class  two.   Class  two  differs  from  class  one  only  in  that  members  have 
local  defects  that  cause  freezing  of  the  domain  growth  at  absolute  zero.   Such  local 
defects  are  expected  to  have  activation  energies  that  are  independent  of  domain  size 
L.   Classes  three  and  four,  on  the  other  hand,  have  defects  with  L(t)  dependent 
activation  energies,  leading  to  logarithmic  growth  which  is  discussed  in  connection 
to  equation  2.12.   Type  Three  and  four  systems  are  differentiated  by  the  power  m  of 
the  logarithmic  growth.   Class  three  has  m=l,  and  is  likely  to  contain  the  random 
field  Ising  model.   Class  four  has  m^l,  and  may  include  dilute  ferromagnets  and  spin 
glasses. 

The  integrated  intensity  is  independent  of  the  existence  of  anti— phase 
domains.^     ''  Experimentally,  however,  if  the  domains  are  very  small,  it  is  possible 
to  lose  some  intensity  in  the  tails  of  the  reflection.   The  width  of  a  superlattice  peak 
is  of  course  greatly  affected  by  the  existence  of  domains.   Broadening  from  the  finite 
size  of  domains,  in  the  Cu^Au  case,  has  a  dependence  on  (h,k,l)''     '  as  will  be  shown 
in  the  next  chapter.   An  effective  size  can  be  defined  but  it  is  purely  formal  since 
there  are  corrections. 
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Lg£f=2T/<r  (2.21) 

where  a  is  the  width  in  units  of  inverse  distance. 


CHAPTER  3 
X-RAY  DIFFRACTION  BY  CugAu 


Basic  X-ray  Scattering 
The  diffraction  of  x— rays  depends  on  the  crystal  structure  and  the 
wavelength.   X-rays  are  scattered  by  the  electric  field  created  by  the  electrons. 
Hence,  the  diffraction  of  x— ray  photons  contains  information  concerning  electron 
density. 

The  position  of  the  diffraction  peak  can  be  easily  derived  by  considering  the 
difference  in  path  length  of  two  beams  reflected  off  two  parallel  crystalline  planes  a 
distance  d  apart. ^     ''  If  the  difference  in  path  length  is  a  multiple  of  the  wavelength 
of  the  radiation  used,  then  constructive  interference  forms  a  peak.   This  logic  leads 
to  Bragg's  Law, 

2dsin(^)=nA  (3.1) 

where  9  is  measured  from  the  plane,  and  n  is  an  integer.   Bragg's  Law  is  valid,  if 
A<2d,  which  is  why  x— rays  are  used  to  probe  crystals.   Bragg's  Law  considers  the 
periodicity  of  the  lattice.   However,  it  does  not  consider  the  structure  or  material 
that  make  up  the  lattice.   Consequently,  Bragg's  Law  does  not  contain  any 
information  about  the  intensity  of  the  diffraction. 

Before  considering  the  electron  density,  which  truly  defines  the  diffraction,  it 
is  convenient  to  define  a  reciprocal  lattice.  Given  three  primitive  lattice  vectors 
(a,  ,a2,a2),  the  primitive  reciprocal  lattice  vectors  are  defined  as; 
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The  primitive  reciprocal  lattice  vectors  are  constructed  so  that  each  is 

orthogonal  to  two  of  the  three  crystal  primitive  lattice  vectors.  Thus, 

b.-a.  =  2t(J..  (3.3) 

I     J  ij 

where  5  ■  ■  is  the  Kronecker  delta  function.  The  factor  of  2t  is  omitted  by 
crystallographers  but  are  used  for  convenience  in  physics.   The  reciprocal  lattice  is 
the  set  of  points  defined  by  the  reciprocal  lattice  vector, 

G  =  v^b^+V2b2+V3b3  (3.4) 

where  v  •  is  an  integer.   All  crystal  structures  have  two  interdependent  lattices,  the 
crystal  lattice  and  the  reciprocal  lattice.   For  example,  the  crystal  lattice  for  CugAu 
is  face  centered  cubic  (fee),  it  is  a  simple  mathematical  exercise  to  show  that  the 
reciprocal  lattice  associated  with  an  fee  crystal  lattice  is  a  body  centered  cubic  (bcc) 
lattice.(«°.^8) 

The  properties  of  a  crystal  are  invariant  under  a  transformation  between 
equivalent  lattice  sites.   The  electron  number  density  n(r)  must,  therefore,  be 
periodic  in  r  with  periods  a.^,X2,a^  in  the  direcrions  of  the  crystal  axes.*-     ''   This  is 
expressed  mathematically  as, 

n(r+T)  =  n(r)  (3.5) 

where  T  is  the  crystal  lattice  vector, 

T  =  n^a^+n2a2+n3a3  (3.6) 

where  n  •  is  an  integer. 

Applying  fourier  analysis  to  the  periodic  function  n(r)  is  equivalent  to  finding 
a  set  of  vectors  J  such  that  n(r)  is  invariant  under  translations  of  the  crystal  lattice 
vector; 

n(r)  =  Sn,.exp(zJ-r)  =  n(r+T).  (3.7) 

J    *' 

By  Bloch's  theorem  G  exactly  satisfies  the  condition  expressed  in  3.7  for  J. 


22 


The  fourier  coefficients  n^  determine  the  intensity  of  the  scattered  radiation, 

and  that  the  Bragg  condition  is  equivalent  to  the  difference  between  the  incident 

and  scattered  wave  vectors  being  G.   As  seen  in  figure  4  the  difference  in  phase 

factor  between  the  incident  (k)  and  scattered  (k')  wave  vectors  is  exp[i(k— k')-r]. 

The  scattering  amplitude  is  an  integral  over  space  of  the  electron  concentration 

times  this  phase  factor; 

F  =  J*dV  n(r)  exp[i(k-k')-r].  (3.8) 

Repeatedly  using  the  invariance  of  n(r)  by  translations  T,  the  scattering  amplitude 

F  can  be  written  as, 

F  =  S  JdV  n^  exp[i(G-Ak).r]  (3.9) 

G 

where  k+Ak  =  k' .   For  (G  =  Ak)  F  is  a  maximum  and  F  =  V  S  n^,   showing  that 

G    '^ 

the  scattering  amplitude  is  a  function  of  the  fourier  coefficients  n^.   The  width  of  a 
peak  for  a  one  dimensional  lattice  goes  as  the  inverse  of  the  number  of  atoms,''     '' 

(r  =  2T/M  (3.10) 

where  M  is  the  number  of  atoms. 

It  is  convenient  to  rewrite  3.9  at  G=Ak  as  the  integral  over  one  cell 
multiplied  by  the  number  of  cells.  The  scattering  amplitude  for  one  cell  is  known  as 
the  structure  factor  F^.^     ^ 

Fq  =  f     dV  n(r)  exp(-iG-r)  (3.11) 

The  electron  concentration  function  can  also  be  written  as  the  superposition  of 
electron  concentrations  n-  Gf  each  atom  of  a  cell.''     '' 

n(r)=S    n.(r-r.)  (3.12) 

j=l    J       ^ 

where  s  is  the  number  of  atoms  in  a  cell.   The  scattering  amplitude  then  becomes, 

Fq  =  S  JdV  n.(r-r.)exp(-zG-r)  (3.13a) 

=  S  exp(-iGT.)  f  dV  n.(p)exp(-iG-p)  (3.13b) 

j  J-'  J 
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^Crystal 


Ikon 


Incident 


Outgoing 


FIG.  4.  Difference  in  path  length  between  two  wave  vectors  scattered  from 
points  0  and  r.   The  total  difference  in  phase  angle  is  (k— k')-r. 
The  corresponding  difference  in  phase  between  wave  vectors  scattered 
froni  dV  at  r  and  0  is  exp[z(k— k')-r]. 
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where  s  =  r— r..   Defining  the  atomic  form  factor  f ., 

f  =  JdV  n.(s)exp(-zG-8)  (3.14) 

rent  that  f  • 
becomes, 


it  is  apparent  that  f.  contains  the  same  information  as  n-  and  the  structure  factor 


FQ  =  Sf.exp(-iGT.).  (3.15) 

The  scattering  amplitude  is  usually  written  in  terms  of  the  integers  v,  ,v„,v„, 
replaced  by  h,k,l  and  r-  in  units  of  the  crystal  lattice  {r-  =  xa, +ya^+za«}, 

Fq  =  S  f.exp[-2Tz(hx.+ky.+lz.)].  (3.16) 

The  intensity  is  proportional  to  F  F,  where  F    denotes  the  complex  conjugate  of  F. 
Therefore,  Fp  does  not  need  to  be  real,  and  hence,  is  not  directly  observable.   In 
addition,  F--,  is  defined  in  terms  of  a  cell  or  basis,  where  there  is  freedom  in  the 
choice  of  a  cell. 

The  scattered  intensity  is  more  naturally  described  in  terms  of  a  correlation 
function.   The  time  independent  pair  correlation  function  is  defined  in  real  space  as, 

G(r)=^/<Kr'-r)p(r')>d^r'  (3.17) 

where  p{t')  is  the  charge  density  at  position  r',  and  <  >  represent  a 
thermodynamic  average.   G(r)  is  an  equilibrium  quantity.   The  fourier  transform  of 
G(r)  is  commonly  called  the  structure  factor,  S(Q). 

S{Q)=;^f^M^Q-^)GW^^  (3.18) 

If  the  system  of  interest  is  in  a  non-equilibrium  state  the  structure  factor  becomes  a 
function  of  both  position  and  time,  S(Q,t). 

Some  confusion  may  arise  between  structure  factor  S(Q,t)  and  the  time 
dependent  equilibrium  structure  factor  normally  written  S(Q,t)  which  refers  to 
fluctuations  about  equilibrium  as  opposed  to  the  approach  to  equilibrium.   Here 
S(Q,t)  is  a  non— equilibrium  S(Q)  at  a  time  t,  the  time  is  added  to  emphasize  the 
fact  that  we  are  measuring  a  non— equilibrium  function.   Another  confusion  in 
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definitions  comes  from  the  scattering  amplitude  F  historically  being  called  a 
structure  factor.   In  older  text  books  on  x— ray  scattering  F  is  commonly  called  a 
structure  factor.   Although  F  F  and  S(Q)  are  interchangeable,  physicists  generally 
use  S.   The  newer  theories  use  S  and  the  older  theories  often  use  F.   Since  1954, 
following  a  paper  published  by  Leon  Van  Hove,  physicists  have  preferred  to  use 
S(Q).  In  Van  Hove's  paper  equation  3.18  is  derived  in  the  Born  approximation.^     ' 
In  the  above  discussion  and  throughout  this  work  S(Q,t)  is  a  non-equilibrium 
quantity,  and  F  will  be  call  the  scattering  amplitude. 


Diffraction  with  Cu»Au  and  its  Relation  to  the  Order  Parameter 


For  an  fee  lattice  the  usual  choice  of  a  basis  is  one  corner  atom  at  [0,0,0],  and 

three  face  centers  at  [O.^io])  [^''^'d'  ^^'^  [2'2''^1'   ^°^  *^^  choice, 

F(h,k,l)  =  f^^^^{l+exp[-Tz(k+l)]+exp[-Tz(h+l)]+exp[-Ti(h+k)]}.   (3.19) 
If  the  indices  (h,k,l)  are  unmixed  (all  even  or  all  odd)  ^n^^^^ir.m-  ^^  *^^  indices 
are  mixed  (2  even  1  odd,  or  1  even  2  odd)  two  of  the  exponents  will  be  odd  multiples 
of  — TZ  and  one  even,  so  Fp  will  vanish. 

The  above  assumes  that  only  one  type  of  atom  makes  up  the  crystal.  This  is 
not  the  case  for  Cu«Au.  The  values  of  f  a     and  fp    can  be  found  in  the 
international  tables  for  x-ray  crystallography  Vol.  3.''       ''  Since  a  real  sample  of 
Cu«Au  is  never  fully  ordered,  a  partially  ordered  sample  should  be  considered.   The 
following  discussion  is  a  long— range  order  theory  only,   and  is  equivalent  to  the 
Bragg— Williams  theory>  ''   Starting  with  the  fraction  of  Cu  and  Au  atoms  located 
at  corner  and  face  center  sites  a  definition  of  long  range  order  [L]  can  be 
introduced.^     ''  This  long  range  order  variable  L  is  not  the  same  as  a  long  range 
order  parameter  used  in  statistical  mechanics  even  though  it  was  originally  named 
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as  such.   Here  long  range  order  L  is  defined  in  terms  of  the  occupancy  of  lattice  sites. 

r    =  fraction  of  Au  atoms  on  a  cube  corner 
w    =  fraction  of  Cu  atoms  on  a  cube  corner 
r,  =  fraction  of  Cu  atoms  on  a  face  center 
Wr  =  fraction  of  Au  atoms  on  a  face  center 

where  r  stands  for  right  type  of  atom,  w  for  wrong  type,  c  for  corner,  and  f  for  face 
center.   Of  these  four  variables  three  are  dependent,  hence,  a  few  results  are 
immediately  obvious, 

r(,+w^  =  1  ;       r^+w^  =  1  (3.20) 

1       3 

'i^c'^t^{~  ^An    the  fraction    of  (3.21) 

gold  atoms 

3      1 

jr^+^w^  =  ip      the  fraction  of  (3.22) 

copper  atoms 

1         3 
where  the  fractions  j  and  j  refer  to  the  concentration  of  corner  and  face  center  sites 

respectively.   Long  range  order  can  now  be  defined, 

L  =  r^+rpl  =  r^-Wf  =  r^^w^  (3.23) 

The  order  variable  defined  this  way  is  one,  if  the  sample  is  fully  ordered,  and  zero,  if 

the  sample  is  fully  disordered.  Using  3.19  or  3.20  the  last  dependent  variable  can  be 

eliminated, 

^  =  (V  ^Cu)/  I  =  (^r  ^Au)/  \  ■  (3-24) 

The  scattering  amplitude  for  crystalline  Cu^Au  is  a  function  of  the  occupation  of 

sites  and  the  atomic  form  factors  for  Cu  and  Au  as  defined  in  equation  3.19. 

Therefore,  the  relation  between  the  scattering  amplitude  and  the  order  variable  is 

easily  derived.''     -'  Dividing  the  sum  in  3.18  up  into  sums  over  cube  corners  and 

face  centers  gives, 

corners  ^  ^^      c  ^..u  c        c      c 


+  S  (  rffj.^+w^f^^)exp[2T2(hXf+kyf+lZf)].    (3.25) 

I  aC6S 
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The  only  corner  is  at  [0,0,0].  Substituting  the  corner  position  for  [h,k,l]  as  well  as 
the  face  center  positions  as  done  to  get  3.16  it  is  found  that  for  unmixed  (h,k,l), 

F  =  (^c^Au+^c^Cu)  +  3(rffc^+Wff^^)  (3.26a) 

=  (^c+3Wf)fAu+  (^c+3rf)fcu  (3.26b) 

=  4(i^c+^f)fAu+  '(^Ai^^Cn  (3-26C) 

=  4(^AufAu+^Cufcu)  (3.26d) 

and  for  mixed  (h,k,l), 

F  =  (^cUu+Vcu)-(¥cu+Vau)  (3.27a) 

=  ( V^f)fAu-  (^r-c)fcu  (3-27b) 

=  Z(fAu-fcu)  (3-27C) 

where  fxr  is  the  atomic  form  factor  for  N,  and  x^  is  the  atomic  fraction  of  atom  N. 
For  stoichiometric  Cu^Au,  x^^  =  1/4,  and  x^^  =  3/4.   The  unmixed  indices  or 
fundamental  peaks  are  just  those  for  an  fee  crystal  with  an  amplitude  corresponding 
to  scattering  from  an  average  atom  [  i^^Q^+^p^^)/"^  ^^^  stoichiometric  Cu^Au].   The 
scattered  intensity  is  the  structure  factor  multiplied  by  its  complex  conjugate. 
Therefore,  the  intensity  of  a  mixed  indices  or  superlattice  reflection  is  proportional 
to  the  square  of  the  difference  between  the  atomic  form  factors  for  Cu  and  Au, 
multiplied  by  the  square  of  the  order  variable.   Hence,  from  a  comparison  of  the 
intensity  of  a  superlattice  peak  relative  to  a  fundamental,  the  long  range  order 
variable  (Z)  can  be  theoretically  determined. 

In  real  applications  there  are  several  effects  which  can  produce  sizable 
errors.'  ^   Extinction  effects  from  different  grains  can  reduce  the  intensities  of 
fundamentals,  hence,  high-order  peaks  should  be  used  in  the  calculation. 
Nevertheless,  the  intensity  of  a  superlattice  reflection  is  proportional  to  the  square 
of  the  long  range  order  variable. 
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In  the  modern  theories  a  proper  order  parameter  must  be  used.   Often  the 
systems  investigated  are  equivalent  to  Ising  models,  and  therefore,  have  a  one 
component  order  parameter.   For  the  COP  process  of  phase  separation  the  order 
parameter  is  the  local  concentration  of  one  of  the  two  types  of  atoms.   This  is 
equivalent  to  an  Ising  model  with  Kawasaki  (spin  exchange)  dynamics  with  the 
order  parameter  defined  as  the  total  magnetization.   For  order-disorder  transitions 
the  analogous  Ising  model  also  uses  Kawasaki  dynamics  only  now  the  order 
parameter  is  the  magnetization  on  a  sub-lattice,  which  is  not  conserved.   Cu,Au 

can  be  described  by  a  three  component  order  parameter.^^^'^^^)  In  a  detailed 

(92) 
theory  of  Lai^     ^  the  three  component  order  parameter  is  defined  in  terms  of  the 

occupation  of  sites  on  one  of  four  interpenetrating  simple  cubic  sublattices  described 
by  Ising  like  occupation  numbers  s  •, 

g  _ +1  if  occupied  by  Cu 
i     —1  if  occupied  by  Au  (3.28) 

where  i  =  1,  2,  3,  and  4  corresponding  to  the  four  equivalent  ways  of  defining  a 
cube  corner  in  the  FCC  lattice  or  equivalently  the  four  cubic  sublattices.   On  the 
basis  of  symmetry  Lai  constructs  the  order  parameter  according  to  the  Landau 
Lifshitz  rules, 

fl(R)=^[Sj(R)-^2(I^H3(R)+S4(R)]  (3.29a) 

f2(R)=2o[ii(R)-^2(^)+S3(RH4(R)]  (3.29b) 

f3(R)=^[Sj(R)+i2(RH3(RH4(R)].  (3.29c) 

Using  conservation  of  atoms  a  conserved  concentration  f  can  be  written, 

f  =  2o[^i(R)+S2(R)+S3(R)+S4(R)-2].  (3.30) 

From  here  an  equation  of  motion  is  found  using  a  Ginzburg-Landau  potential.   The 

square  modulus  of  the  order  parameter  used  in  modern  theories  is  proportional  to 

the  intensity  of  a  superlattice  peak. 
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When,  as  is  the  case  for  binary  alloys,  the  order  parameter  is  proportional  to 
the  charge  density  (atomic  concentration),  the  charge  density  correlation  function  is 
directly  proportional  to  the  order  parameter  correlation  function.   Therefore,  the 
structure  factor  measured  in  the  scattering  experiment  is  linearly  proportional  to 
the  order  parameter  correlation  function  which  is  normally  modeled  theoretically. 

Anisotropy  in   Cu^Au 

The  solid  Cu„Au  solution  forms  a  fee  crystal  and  below  T    z  390°C  orders  in 
the  LI2  structure  (Au  on  the  corners,  Cu  on  the  face  centers  as  in  figure  1)  with  a 
ground— state  degeneracy  of  p=4.   The  four  differe'nt  phases  correspond  to  redefining 
the  cube  corners  as  one  of  its  three  nearest  face  centers,  and  of  course  the  original 
phase  one  started  with.   Mathematically  one  shifts  the  position  of  the  ordered  state 
by  0,  (al+a2)/2,  (aH-a3)/2  or  (a2+a3)/2.   The  presence  of  domain  walls  leads  to 

anisotropic  line  shapes  in  superlattice  peaks  as  discussed  below. 

(21) 
Two  types  of  domain  walls  exist  in  Cu„Au.   Type— 1  walls^     ''  are  easily 

formed  due  to  their  low  energy.   These  walls  can  be  formed  by  displacing  the  atoms 

on  one  side  of  a  plane  perpendicular  to  the  aS  axis  by  (aH-a2)/2  as  in  figure  5. 

Type— 1  walls  can  also  be  formed  by  any  cyclic  permutation  of  (al,a2,a3)  in  the 

above  method.   A  Type-2  wall  can  be  formed  by  displacing  the  atoms  on  one  side  of 

a  plane  perpendicular  to  a3  by  (al+a3)/2,  as  in  figure  6.   The  important  difference 

between  type— 1  and  type— 2  walls  is  that  type— 1  walls  conserve  the  nearest  neighbor 

configuration  where  type— 2  walls  do  not.   The  type— 2  wall  in  figure  3  has  gold 

atoms  as  nearest  neighbors,  one  could  form  a  slightly  different  type— 2  wall  say 

type— 2b  by  redefining  all  the  Au  atoms  to  be  Cu  atoms  and  all  the  in  plane  Cu 

atoms  to  be  Au.   Type— 2b  walls  have  Cu  atoms  as  nearest  neighbors  in  the  wall 
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Type-1   Donaln   Wall 
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FIG.  5.  Domains  of  different  phase  separated  by  a  type-  1  domain  wall  as 
described  in  the  text.   The  line  marks  the  position  of  the  wall. 
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Type-2   Domain   Wall 
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FIG.  6.  Domains  of  different  phase  separated  by  a  type—  2  domain  wall  as 
described  in  the  text. 
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instead  of  Au  atoms  as  in  type— 2a  walls.   Moreover,  if  one  looks  in  a  direction 
perpendicular  to  a  type-2  wall  one  finds  planes  of  all  Cu  alternating  with  planes  of 
half  Au  and  half  Cu,  the  difference  between  type— 2a  and  type— 2b  walls  is  just  which 
of  the  two  types  of  planes  is  repeated  to  create  the  wall.   In  order  to  move  a  type-2 
wall  there  must  be  some  diffusion  of  atoms  due  to  the  off  stoichiometry  of  the  wall. 
On  the  other  hand  a  type— 1  wall  can  be  moved  by  shifting  the  wall  up  or  down  by 
(aH-a2)/2,  hence,  it  is  often  called  a  half-diagonal  glide  domain  wall.  This 
difference  posed  the  question  of  whether  the'two  types  of  walls  grow  at  the  same 
rate  or  not.  If  the  two  walls  do  not  grow  at  the  same  rate  then  scaling  would  not  be 
valid. 

One  way  to  measure  the  anisotropic  profile  is  to  look  in  both  the  radial  and 
transverse  directions  during  coarsening.   This  corresponds  to  observing  the  motion 
of  both  types  of  walls.   At  the  [1,0,0],  the  radial  direction  probes  type  1  walls  and 
the  transverse  type  2.   At  the  [1,-1,0],  the  situation  is  almost  reversed,  the  radial 
direction  probes  type  2  walls  and  a  specific  transverse  direction  probes  type— 1.   So 
another  way  to  consider  both  types  of  walls  is  to  observe  the  [1,0,0]  and  the  [1,-1,0] 
in  the  radial  (or  transverse)  direction.   In  this  study  the  validity  of  scaling  during 
coarsening  in  an  anisotropic  system  is  established.   The  fact  that  both  walls  behave 
in  the  same  way  suggest  it  is  the  same  mechanism  driving  the  growth  in  all 
directions.  In  figures  5  and  6  straight  domain  walls  have  been  shown,  however,  the 
growth  is  curvature  driven  suggesting  that  a  curved  wall  should  be  considered.   A 
curved  wall  is  shown  in  figure  7,  in  it  one  can  find  elements  for  all  types  of  walls. 
Although  it  is  possible  to  envision  this  wall  growing  with  different  velocities  in 
different  directions,  that  has  been  found  not  to  be  the  case. 
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Curved  Domain   Walls 
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FIG.  7.  Domains  of  three  different  phases  separated  by  curved  domain  walls 
constructed  from  a  combination  of  both  type-  1  and  type—  2  walls. 
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The  Effect  of  Type   1  Walls  on  Scattering. 

The  width  of  a  superlattice  peak  is  of  course  greatly  affected  by  the  existence 
of  domain  walls.   In  the  following  consideration  of  the  scattering  a  number  of 
assumptions  are  made:  1.  type  2  walls  are  not  present,  2.  the  half  diagonal  glide 
domain  boundaries  perpendicular  to  each  of  the  axes  are  completely  independent 
from  each  other,  and  3.  the  intensity  drops  off  sharply  away  from  the  peak  (a  small 
q  approximation).   For  convenience  the  convention  that  all  superlattice  peaks  be 
written  such  that  h  and  1  have  the  same  parity  has  been  adopted.   The  shape  of  the 

peak,  as  shown  in  appendix  A,  can  then  be  written  as  (in  reduced  units), 

2 
cy     N,a         sin  (rN^h,,)        N„fl 

I(h  h  h  )  =  F^^-^ 2 -1 2-^ 2    (3-31) 

where  h,    h^  and  h^  represent  the  departure  of  the  diffraction  vector  from  the 
[h,k,l]  peak  position  [ie.,q  ={2T(l^)'{)x^,]x^,h^),  Q  =G+Q,  G  =(2T/a)-(h,k,l)],  and  a 
is  the  probability  of  crossing  an  anti— phase  domain  boundary  in  a  distance  a.   The 
intensity  for  a  superlattice  peak  as  given  above  is  roughly  a  thin  disk  lying  in  the  h, 
hg  plane  with  a  thickness  determined  by  the  crystal  dimension  N„a^.   type  2  walls 
would,  of  course  thicken  the  disk.   Figure  8  shows  the  orientations  of  theses  disks  for 
one  plane  of  reciprocal  space.   Also  this  equation  shows  that  the  radius  of  the  disk 
increases  with  the  probability  a  of  anti— phase  boundaries. 

Normally,  the  condition  I(h,h2h2)  equal  to  (l/2)I(h^h^h„)  would  be  used  to 
define  the  width  of  the  peak  in  terms  of  L.   However,  the  approximation  of  h,h^h„ 
being  small  lead  to  problems.   Therefore,  one  starts  with  a  more  general  equation 
for  the  width.(^°) 

CD 

/   I    dq 

FWHM  «  ~^ (3.32) 
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FIG.  8.   Shapes  of  the  reflections  in  reciprocal  space  for  type  1  walls 

normal  to  the  three  axes.   The  fundamental  reflections  are  spheres, 
and  the  superlattice  reflections  are  disk-shaped. 
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Also  in  appendix  A,  the  relation  between  domain  size  and  superlattice  peak  width  is 

shown  from  3.32  to  be,^^°^ 

/  I  I.  I  I  I  1.  I  \ 

(3.33) 


FWHM  ot  (2r/L)     (  I  ^  I  +  I  ^  I  ) 
I      *) 5 T 


This  result  is  for  type  1  walls,  however,  type  2  walls  also  produce  broadening  with  a 
2t/L  dependence.  The  correction  of  2.21  for  the  [0,1,0]  type  peak  is  one,  and  for  the 
[1,0,T]  is  \fl.  For  a  given  functional  form  equation  3.32  can  be  used  to  relate  the 


corresponding  width  to  L.  For  a  gaussian  with  width  c,  defined  as, 

F(x)=A.exp[[(x-x^)/2^]2^ 
equation  3.33  gives, 


<T  =  (V^/L). 


(ih_i  +  m) 
/h'^k'^r^ 


(3.34) 


(3.35) 


where  as  for  a  lorentzian  squared  with  width  T,  defined  as, 

-2 


F(x)= 
equation  3.33  gives, 


^-^ 


(x-xj/r 


.2 


r  =  (v^/L)   (MMKm2 


(3.35) 


(3.36) 


The  reader  should  be  aware  of  the  above  corrections  for  type  1  walls.   However, 
throughout  this  work  it  is  assumed  that  2.21  is  correct  for  both  types  of  walls. 


CHAPTER  4 
EXPERIMENTS 


The  Bulk  Single  Crystal 
Experimental  Details 

The  crystalline  sample  of  Cu«Au  used  in  these  experiments  was  cut  at  Oak 
Ridge  National  Laboratories  to  expose  the  [h,0,0]  face.   In  order  to  remove  any 
oxidation  the  surface  was  etched  in  nitric  acid  and  then  annealed  for  12  hours  at 
700°C  under  an  Ar  atmosphere. 

The  furnace  used  in  these  experiments  is  composed  of  an  evacuable  chamber 
and  a  heating  element.   Heating  was  performed  by  running  AC  current  through  a 
0.003mm  diameter  tungsten  coil  wrapped  around  a  9mm  diameter  copper  rod  that 
had  been  covered  with  a  1.5mm  thick  layer  of  electrically  insulating  boron  nitride. 
The  tungsten  coil  consisted  of  20  to  25  loops  covering  a  length  of  35mm.  The  sample 
was  mounted  to  the  rod  as  shown  in  figure  9.   This  heating  rod  is  placed  in  the 
chamber  which  is  equipped  with  a  360°  Beryllium  window.   Beryllium  is  a  light 
metal  that  is  almost  transparent  to  8  KeV  x-rays.   The  tantalum  shown  in  figure  9 
was  used  to  keep  Cu  from  diffusing  from  the  mount  into  the  sample.  Tantalum  foil 
was  also  placed  above  the  copper  rod  for  the  same  reason. 

A  proportional  relay  (Douglas  Randall  model  RD04A)  was  used  to  control 
the  current  to  the  heating  coil,  and  hence,  the  temperature.  The  proportional  relay 
was  connected  to  a  model  828D  Micristar  temperature  controller  which  used  a 
K-thermocouple  placed  at  the  top  of  the  copper  rod  to  get  a  reference  temperature. 
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FIG.  9.  Heating  element  and  sample  mount  used  in  connection  with  the  bulk 
sample,  the  sputtered  films,  and  the  4500  A  (M.B.E.)  film. 
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In  the  initial  experiments  quenches  were  done  by  putting  the  Micristax  on  manual 
which  fixed  the  current,  then  the  Micristar's  temperature  set  point  was  changed  to  2 
to  5  C  above  the  final  post  temperature  desired.   Next  the  current  was  turned  off  at 
the  source  allowing  the  sample  to  cool.   As  the  post  temperature  reached  the 
Micristar's  new  set  point  the  current  source  was  turned  back  on  and  the  Micristar 
put  on  automatic.   The  system  would  soon  reach  an  almost  constant  temperature  a 
few  degrees  below  the  set  point.   If  left  alone  the  system  temperature  would  slowly 
be  raised  till  it  reached  the  set  point.   Therefore,  the  Micristar  was  temporarily  put 
on  manual  and  the  set  point  lowered  to  the  post  temperature.  The  Quenches  took 
place  in  30  to  90  seconds  depending  on  depth.   Quench  times  were  measured  from 
when  the  temperature  crossed  T    to  when  the  temperature  was  within  one  degree 
Celsius  of  the  final  temperature.   The  quenches  did  not  over  shoot  but  smoothly 
went  to  the  final  temperature.   The  quenches  were  significantly  improved  in  later 
experiments  as  will  be  disused  in  connection  with  the  M.B.E.  films. 

In  addition  to  the  thermocouple  used  for  temperature  control  there  also  was 
a  K— thermocouple  held  directly  onto  the  face  of  the  sample.   The  temperature 
difference  between  the  sample  and  the  top  of  the  heating  rod  was  about  9  C.  The 
vast  majority  of  this  gradient  existed  in  the  tantalum  plate  protecting  the  sample 
from  the  Cu  plate.   The  gradient  across  the  sample  is  believed  to  be  small  due  to  the 
geometry  of  the  heater.   Moreover,  an  estimate  of  the  gradient  suggests  that  the 
temperature  gradient  across  the  sample  was  O.OlC/mm.   This  estimate  was  based  on 
the  thermal  conductivity  of  Cu  and  Au  and  the  Stefan— Boltzmann  law,  weighted  by 
the  emissivity  of  Cu  and  Au.   The  x— ray  spot  size  on  the  sample  was  roughly 
lx6mm  ,  consequently  the  temperature  difference  within  the  area  of  interest  was 
about  0.06°C.   The  temperature  control  was  good  to  ±0.3  C  over  several  hours  and 
±0.5°C  over  several  days.   Therefore, the  gradient  was  negligible  relative  to  the 
temperature  stability. 
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In  order  to  determine  the  transition  temperature  the  intensity  of  the  [0,1,0] 
peak  was  recorded  as  a  function  of  time  for  several  temperatures  after  raising  the 
temperature  from  below  T    to  near  T  .   The  intensity  at  each  temperature  was  then 
plotted  as  a  function  of  inverse  time  and  the  curve  was  extrapolated  to  1/time  =  0. 
At  385.2  C  the  peak  was  observed  for  3.5  hours  and  the  extrapolation  suggested 
that  the  sample  would  almost  disorder  but  not  completely.   On  the  other  hand  at 
385.7  C  the  peak  was  nearly  gone  after  1.5  hours  and  the  extrapolation  clearly 
showed  that  the  sample  was  above  T  .   The  transition  temperature  was  then  taken 
to  be  the  value  of  a  weighted  average  based  on  the  extrapolations,  and  found  to  be 
385.3  C.  This  extrapolation  is  not  very  important  because  the  temperatures  used 
are  only  half  a  degree  apart,  but  it  was  clear  that  the  transition  was  closer  to 
385.2±0.5°C  than  385.7±0.5°C. 

The  wavelength  of  the  x-rays  used  is  named  from  the  core  electron  levels  of 
the.  metal  target  used  in  the  Rigaku  RU300.   The  doublet  line  of  Cu  is  referred  to  as 
CuK*^  (;  =  1.54178X)  or  just  K^.   The  two  components  of  this  doublet  line  are  K? 
(A  =  1.5405lX)  and  K2  (^=1.54433^). 

The  X— rays  were  produced  by  a  rotating  anode.   The  major  components  of  a 
rotating  anode  are,  the  anode  and  cathode  or  filament.  To  produce  x-rays,  electrons 
are  accelerated  by  a  potential  difference  from  the  filament  to  the  anode  or  target. 
In  so  doing  x-rays  are  produced  in  several  ways.   The  first  is  from  bremsstrahlung, 
the  sudden  deceleration  of  the  electrons.   The  second  is  through  the  process  of 
exciting  core  energy  levels  of  the  metal  target,  which  then  decay  producing  x— rays 
of  specific  wavelengths.   Almost  all  of  the  work  done  here  used  a  Cu  target.   The 
most  intense  core  energy  levels  are  labeled  k"  and  K^.   This  method  of  x-ray 
production  is  very  inefficient,  only  about  one  percent  of  the  energy  used  goes  into 
the  production  of  CuK    radiation.   The  majority  of  the  energy  goes  into  heating  the 
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target.   Therefore,  the  target  must  be  cooled.   Tube  sources  have  a  fixed  target  that 
is  typically  cooled  by  running  water.  Tube  sources  at  best  are  run  at  around  3kW. 
In  order  to  increase  the  intensity  the  anode  (or  target)  is  rotated  increasing  the  area 
that  is  heated.   In  such  cases  the  anode  is  in  the  shape  of  a  cylindrical  bowl  that  is 
rotated  about  the  axis  of  the  cylinder  with  the  electrons  bombarding  the  side. 

For  geometrical  reasons  the  resolution  or  spot  size  is  best  at  smaller  angles. 
An  angle  of  six  degrees  as  measured  from  the  surface  of  the  target  was  used  here. 
Six  degrees  gives  a  1  to  10  reduction  of  the  spot  size.   Other  commonly  used  angles 
are  3,  9,  and  12  degrees. 

The  RU300  can  be  operated  in  two  different  geometries  (point  or  line  focus). 
If  the  filament  is  mounted  perpendicular  to  the  scattering  plane  associated  with  the 
Huber  diffractometer,  then  the  anode  must  spin  in  the  scattering  plane.   The 
electrons  must  then  hit  the  anode  in  a  line  perpendicular  to  the  scattering  plane. 
This  produces  a  bar  or  line  of  x— rays  perpendicular  to  the  scattering  plane.   If  the 
situation  is  reversed  then  the  bar  lies  in  the  plane  at  an  angle  of  6  degrees  to  the 
direction  of  the  incident  beam.   The  projection  of  this  bar  onto  a  cross  section  of 
beam  perpendicular  to  its  motion  is  a  point.   For  a  PSD  a  line  focus  has  better 
instrumental  resolution  than  a  point  focus,  but  the  point  focus  generally  gives  higher 
counting  rates,  because  more  of  the  x— rays  are  focused  onto  the  sample. 

The  optics  following  the  x— ray  source  differ  according  to  whether  the 
direction  scanned  in  reciprocal  space  was  parallel  to  the  ordering  vector  Q  (radial 
scan)  or  perpendicular  to  Q  (transverse  scan).   Radial  scans  probe  reciprocal  space 
along  a  line  that  passes  through  the  origin  and  the  peak  in  question.   Hence,  in  a 
radial  scan  the  magnitude  of  Q  is  directly  varied.    Mathematically  any  point  in  a 
radial  scan  is  definable  as  (c-h,c-k,c-l)  where  c  is  a  real  number.   For  the  [0,1,0]  a 
radial  scan  is  along  the  [0,k,0]  direction.   Transverse  scans  in  contrast  probe  a 
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direction  perpendicular  to  the  radial  direction.   For  example,  at  the  [0,1,0]  a  scan 
that  varied  h,  1,  or,  a  linear  combination  of  the  two  would  be  a  transverse  scan.   In 
this  study  the  magnitude  of  Q  is  almost  constant  in  a  transverse  scan.   Hence,  a 
transverse  scan  is  almost  identical  to  a  mosaic  scan  in  which  the  sample  is  rocked 
while  Q  is  held  fixed. 

Recalling  the  disk  shaped  scattering  profiles  for  a  superlattice  peak 
previously  shown  in  figure  8  a  transverse  scan  of  the  [0,1,0]  is  along  the  diameter  of 
the  disk  and  a  radial  scan  is  through  the  disk.   Equation  3.31  implies  that  type  1 
walls  do  not  affect  the  width  in  the  radial  direction  at  the  [0,1,0].   Therefore,  the 
radial  direction  probes  type  2  walls  only.   The  transverse  direction  probes  type  1 
and  2  walls.   However,  the  broadening  in  this  direction  is  dominated  by  type  1 
waUs. 

The  sample  was  placed  at  the  center  of  rotation  of  a  Huber  Eulerian  cradle 
by  adjusting  an  x-y  translation  at  the  base  of  the  furnace.   Figure  10a  and  10b 
shows  schematically  the  apparatus  in  the  two  geometries  used;  10a  was  used  to 
collect  the  radial  data  and  10b  the  transverse  data.   Set  up  10a  uses  a  line  focus  and 
a  Braun  model  OEO— 50M  linear  platinum  wire  based  position  sensitive  detector 
(PSD).  The  Braun  detector's  signal  is  binned  into  1024  channels  by  a  multi  channel 
analyzer.  For  the  sample  to  detector  distance  used  (23.5cm)  each  channel 
corresponds  to  0.0125  degrees  in  29.   The  window  on  the  PSD  was  5cm  long  which 
corresponds  to  a  solid  angle  of  10.3  degrees.  The  transverse  set  up  (figure  10b) 
employed  a  point  focus  and  a  Bicron  model  1XMP040B  scintillation  detector.   The 
Bicron's  signal  was  fed  into  an  E.G.&G.  Ortec  model  590A  single  channel  analyzer 
which  is  connected  to  an  E.G.&G.  model  974  quad  counter/timer.   The 
counter/timer  is  in  turn  connected  to  an  IBM  XT  PC  through  an  I.E.E.E.— 488 
interface  bus.   The  PC  also  directs  the  Huber  through  a  Klinger  Scientific  MC4 
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FIG.  10.  Lay  out  of  equipment:  a)  Radial  setup  b)  Transverse  setup. 
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stepping  motor  controller.   Both  set  ups  use  a  second  Bicron  detector,  interfaced 
with  the  PC  as  an  incident  intensity  monitor,  so  that  fluctuations  in  intensity  can 
be  corrected  for. 

For  the  radial  set-up  over  90%  of  the  K?  (A=1.54433A)  intensity  was  slit 
out  leaving  mostly  K?  (A  =  l. 54051a)  where  the  subscripts  define  the  two 
components  of  the  doublet  line.   The  slitting  out  of  K^  was  done  on  the  [0,2,0]  peak 
with  the  pre— sample  slits  shown  in  figure  10a.   At  the  [0,2,0]  K?  and  K^  where  well 
separated  as  opposed  to  at  the  [0,1,0]  where  the  two  components  were  not 
completely  resolved.   It  was  important  to  slit  out  K^  due  to  the  complications  it 
would  have  added  to  the  resolution  correction.   In  the  transverse  set— up  it  was  not 
necessary  to  do  so  because  of  the  analyzer.  The  resulting  instrumental  resolution  for 
both  set-ups,  in  the  corresponding  directions  of  the  scans  in  reciprocal  space  was 
0.003X~\hWHM). 

In  the  above  discussion,  PSD  scans  have  been  called  radial  scans.   A  radial 
scan  in  real  space  involves  varying  the  detector  angle  29  while  half  angling  the 
sample  (^=fl)-   The  major  difference  in  a  PSD  scan  relative  to  a  radial  scan  is  the 
lack  of  half  angling  of  the  sample.   The  two  are  approximately  the  same  for  small 
distances  about  a  peak.   For  the  sample  to  detector  distance  used  the  1024  PSD 
channels  covered  ±6.4  degrees  in  20.  The  maximum  error  in  the  sample  angle  is 
then  3.2  degrees.   Typically  a  range  less  than  400  channels  (or  2.5  degrees)  was  used 
in  the  analysis.   At  the  [0,1,0]  the  radial  direction  probes  the  thin  part  of  S(Q). 
Therefore  it  is  reasonable  to  assume  that  S(Q)  is  constant  (or  flat)  for  small 
deviations  from  the  radial  direction.   Under  this  assumption,  a  PSD  scan  differs 
from  a  radial  scan  in  reciprocal  space  by  a  displacement  of  90  degrees  from  the 
radial  direction  of  (Q— q)-tan(^— fl).   For  the  range  used  in  the  analysis,  the 
difference  in  the  distances  from  Q  between  the  two  types  of  scans  is  less  than  7.7 
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percent  at  the  end  points  and  3.4  percent  at  the  mid— points.   The  mosaic  spread  of 
the  sample  (0.25  degrees)  would  further  lessen  the  difference  between  the  these  two 
types  of  scans.   This  small  difference  between  a  PSD  scan  and  a  radial  scan  only 
affects  the  results  in  that  the  width  of  a  PSD  scan  is  slightly  dependent  on  type  1 
walls  where  a  true  radial  scan  probes  only  type  2  walls. 

Analysis  of  the  Bulk  Data 

One  of  the  purposes  of  this  work  was  to  consider  the  growth  kinetics  in  an 
anisotropic  system.   Noda,  Nishimura,  and  Takeuchi*-     ''  measured  the  [1,1,0] 
superlattice  peak  in  Cu^Au  by  time  resolved  x-ray  scattering  techniques.   They 
found  evidence  for  scaling  of  S(Q,t),  a  Lorentzian-squared  line  shape,  and 
curvature— driven  growth.   However,  they  incorrectly  assumed  isotropic  scattering 
profiles,  casting  some  doubt  on  the  validity  of  their  conclusions.   It  was  postulated 
that  type  2  walls,  which  were  not  considered  in  this  earlier  work,  would  have  a 
slower  growth  due  to  the  need  to  diffuse  one  type  of  atom  relative  to  the  other  in 
moving  a  type  2  wall.   The  off  stoichiometric  nature  of  type  2  walls  that  might  have 
caused  the  growth  to  become  diffusion  limited  is  discussed  in  chapter  3  section  C. 

In  these  experiments  the  intrinsic  x-ray  peak  width  and  shape  is  of  interest 
for  several  reasons  as  discussed  earlier.   Measuring  the  intrinsic  line  shape  of  S(Q,t) 
is  complicated  by  several  extraneous  contributions  to  the  line  width,  arising  from 
the  effects  of,  instrumental  resolution,  finite  grain  size,  strain,  and  the  mosaic 
spread.   The  correction  for  instrumental  resolution  is  the  most  important.   The 
accuracy  of  all  results  depends  heavily  on  a  careful  handling  of  the  resolution.   The 
radial  data  was  corrected  for  instrumental  resolution  as  well  as  finite  size  and  strain 
as  discussed  below.   The  correction  for  the  transverse  direction  was  further 
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complicated  by  the  mosaic  of  the  sample.  In  both  set  ups  the  instrumental 

resolution  was  measured  by  placing  a  perfect  Ge  single  crystal  at  the  sample 

position  and  observing  the  [1,1,1]  reflection,  which  is  at  an  angle  near  that  of  the 

[0,1,0]  of  our  sample.   The  intrinsic  width  of  a  Ge  single  crystal  peak  is  negligible 

compared  to  the  instrumental  resolution.   Estimates  of  broadening  from  sample 

characteristics  other  than  the  ordered  domain  size  are  made  from  the  broadening  of 

fundamental  Cu^Au  bragg  peaks.   Graphs  of  the  measured  instrumental  resolution 

functions  are  shown  in  figures  11  and  12  for  the  radial  and  transverse  set  ups 

respectively.   A  detailed  study,  involving  both  theoretical  and  experimental 

determinations  of  the  instrumental  resolution  function  show  the  complex  nature  of 

this  function(^^^'''^^^^  In  this  work  the  instrumental  resolution  was  carefully 

measured.  The  validity  of  the  procedure  discussed  here  is  supported  by  the  detailed 

study.  The  contribution  to  the  width  of  the  [0,1,0]  arising  from  finite  size  and  strain 

was  estimated  from  radial  scans  of  the  [0,2,0]  and  [0,4,0].   These  scans  were  fit  to 

gaussian  profiles,  and  the  resulting  widths  in  degrees  were  corrected  for  instrumental 

I — 2 2 

resolution  assuming  the  widths  add  in  quadrature  (""cor""^  Veasured~'^res  ^' 

Ultimately  the  effects  of  finite  grain  size  and  strain  can  be  included  in  an 
effective  resolution  function.  Finite  grain  size  leads  to  an  identical  contribution  to 
the  reciprocal  space  line  width  of  each  bragg  peak  in  the  sample.   Strain  can  be 
viewed  as  providing  a  distribution  of  lattice  constants,  leading  to  an  additional 
width  at  each  peak  that  is  linearly  proportional  to  Q.  To  obtain  the  contribution  to 
the  width  at  the  [0,1,0]  arising  from  finite  size  and  strain  the  [0,2,0]  and  [0,4,0] 
corrected  widths  were  converted  into  units  of  inverse  angstroms  and  fit  to  [a  +  bQ], 
where  a  and  b  are  fit  parameters  and  Q  is  the  magnitude  of  the  corresponding 
scattering  vector.   The  widths  due  to  finite  size  and  strain  correspond  to  a  and  bq 
respectively.  The  result  of  this  fit  suggested  that  the  correction  at  the  [0,1,0]  was 
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FIG.  11.  Scan  of  Ge  [1,1,1]  used  to  measure  instrumental  resolution  in  the 

radial  direction.  The  graph  shows  only  the  central  portion  of  the  scan, 
many  more  points  were  taken  in  the  tails. 
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Transverse  Resolution 
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FIG.  12.  Scan  of  Ge  [1,1,1]  used  to  measure  instrumental  resolution  in  the 
transverse  direction. 
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0.002^      (FWHM),  one  third  that  of  the  instrumental  resolution.   One  could  argue 

2 
that  the  widths  should  be  added  as  a+bQ  ,  not  linearly  as  was  done.  If  the  fits  are 

done  this  way  the  correction  for  finite  size  and  strain  becomes  ten  percent  larger  and 

the  overall  correction  1.4  percent  larger.  The  simplest  correction  is  what  was  done, 

this  method  is  equivalent  to  assuming  the  widths  add  like  lorentzian  widths.  The 

case  for  using  a-l-bQ    assumes  the  widths  add  like  gaussian  widths  and  that  the 

width  from  strain  is  small  compared  to  that  for  finite  size,  both  of  which  are  poor 

assumptions.   In  any  case  the  final  result  is  not  sensitive  to  the  choice.   A  high 

resolution  set— up  similar  to  10b  but  with  a  germanium  monochromator  and 

analyzer  was  used  to  measure  the  mosaic  of  the  Cu«Au  sample  at  the  [0,2,0]. 

Mosaic  spread  is  a  measure  of  the  distribution  of  orientations  of  the  grains  in  the 

sample.   The  mosaic  spread  contributes  to  the  each  peak  a  constant  angular  width 

in  real  space,  and  hence,  in  reciprocal  space  is  proportional  to  Q.   The  transverse 

scans  were  taken  as  a  function  of  Q  so  the  positions  in  the  high  resolution  scan  were 

converted  to  Q  at  the  [0,1,0]  after  a  small  correction  for  finite  size  and  strain.   The 

mosaic  at  the  [0,1,0],  shown  in  figure  13,  has  a  FWHM  of  0.0044  radians 

(O.OOtX"  ),  a  value  that  is  a  little  larger  than  instrumental  resolution. 

Initially  the  data  from  the  bulk  sample  was  fit  to  a  lorentzian  squared  form 

and  then  the  width  were  corrected  by  subtracting  the  effective  resolution^      '.  That 

method,  although  useful  for  preliminary  analysis,  introduces  a  systematic  error. 

Unlike  a  lorentzian  and  a  gaussian  the  convolution  of  two  lorentzian  squared 

functions  is  not  a  similarly  shaped  function,  hence,  there  is  not  a  simple  analytical 

way  to  remove  the  resolution  from  a  lorentzian  squared  fit.  In  order  to  correctly 

correctly  account  for  the  resolution  the  data  was  fit  to  the  convolution  of  the 

resolution  function  and  a  model  functional  form.   Although  this  procedure  was  a 

better  method,  attempts  to  assign  a  functional  form  to  the  resolution  failed  to  model 
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FIG.  13.  Scan  of  the  [0,2,0]  of  Cu^Au  with  the  position  linearly 

corrected  to  correspond  to  the  mosaic  spread  at  the  [0,1,0]. 
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the  data  well.   The  reason  for  this  is  simple:  The  various  contributions  to  the 
resolution  were  in  general  equally  poorly  fit  by  a  gaussian  and/or  a  lorentzian 
squared.  Other  forms  fit  even  worse.  The  instrumental  resolution  had  a  slight 
asymmetry  from  slitting,  and  the  mosaic  had  a  little  structure  in  it  (see  figures 
11—12).   Therefore,  fitting  was  done  by  convolving  the  instrumental  resolution  data 
itself  point  by  point  with  a  narrow  gaussian  representing  finite  size  and  stain  and 
the  model  line  shape.   This  method  worked  very  well  for  the  radial  data.   The  radial 
data  was  collected  with  a  PSD  as  described  earlier  which  took  many  points  in  in  the 
radial  direction  minimizing  any  error  from  the  uncertainty  in  the  shape  of  the  peak. 
Figure  14  is  the  result  of  convolving  the  instrumental  resolution  (figure  11)  with  a 
gaussian  representing  finite  size  and  strain.   Moreover,  the  correction  in  the  radial 
direction  did  not  involve  a  mosaic  spread.   In  the  transverse  direction  this  method 
also  worked  well.   However,  the  uncertainty  in  the  line  width,  which  is  dominated 
by  errors  in  the  resolution  correction,  is  larger  in  the  transverse  data  than  the  radial 
data.   The  convolution  of  the  transverse  instrumental  resolution  (figure  12)  with  a 
gaussian  representing  finite  size  and  strain  as  well  as  the  mosaic  (figure  13)  is  shown 
in  figure  15.   Despite  the  larger  absolute  resolution  correction,  the  raw  transverse 
data  shows  changes  in  the  shape  of  the  peak,  and  hence,  S(Q,t)  more  clearly  than 
the  radial  data.  In  fact  the  resolution  is  not  very  important  to  the  change  in  S(Q,t) 
because  the  change  occurs  early  when  the  resolution  correction  is  small  compared  to 
the  peak  width. 

Nucleation  is  associated  with  very  early  times  and  the  equilibration  of  short 
range  order.   The  short  range  order  peak  intensity  is  very  low,  consequently  only  its 
existence  could  be  seen  by  the  lack  of  a  long  range  order  peak.   However,  changes 
during  the  growth  stage  that  are  revealed  through  changes  in  the  shape  of  the 
scattered  intensity  are  seen.  This  change  is  analyzed  by  fitting  the  data  to  both 
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FIG.  14.  Convolution  of  the  instrumental  resolution  (fig.  11.)  with  a  gaussian 

of  width  0.001  K~^  (.002X~-^FWHM)  (representing  finite  size  and  strain. 
The  convolution  represents  the  correction  in  the  radial  data. 
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Transverse  Resolution  Correction 
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FIG.  15.  Convolution  of  the  instrumental  resolution  (Cg.  12.),  mosaic  spread 

(fig.  13.),  and  a  gaussian  of  width  0.001  A~^  (.002^"^  FWHM) 
representing  finite  size  and  strain.   The  convolution  represents  the 
correction  in  the  transverse  direction. 
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theoretical  and  fundamental  functions.   Another  method  for  seeing  changes  in  line 
shape  is  to  scale  peaks  from  different  times  by  their  amplitude  and  width.  In  such 
scaling  plots,  differences  in  shape  are  apparent.   However,  they  do  not  define  or 
systematize  the  change  as  do  the  fits. 

The  fits  to  the  data  are  used  to  extract  a  reciprocal  space  line  width  which  is 
inversely  proportional  to  an  ordered  domain  size  L(t).   L(t)  is  compared  to  various 
functional  forms,  specifically  a  power  law  (equation  2.5)  and/or  a  logarithmic  law 
(equation-2.12).   The  fits  also  provide  an  amplitude  and  integrated  intensity  (width 
multiplied  by  the  amplitude)  as  functions  of  time  and  temperature.   Scaling  the 
intensity  and  time  gives  information  relating  to  changes,  or  lack  thereof,  in  the 
mechanisms  for  growth. 

Results  from  the  Bulk  Data: 

The  shape  of  our  peak  at  early  times  (  first  half  hour  )  is  gaussian  (GS). 
During  this  time  the  superlattice  peak  arises  from  isolated  ordered  clusters 
embedded  in  a  disordered  matrix.  A  GS  shape  at  early  times  is  not  surprising,  after 
all  it  is  associated  with  finite  size  broadening.   At  later  times  the  line  shape  crosses 
over  to  a  Lorentzian-squared  (LSQ)  form.   The  LSQ  shape  is  characteristic  of  the 
coarsening  process,  and  results  from  a  random  distribution  of  ordered  domains 
separated  by  sharp  domain  walls. *>     ^  Figure  16  compares  the  shape  of  the  peak  at 
early  and  late  times. 

The  simplest  model  for  superlattice  peak  line  shapes  in  Cu«Au  is  a 
Hendrichs-Teller^      '  (HT)  approach,  which  assumes  independent  domain  walls  as 
discussed  in  chapter  3  section  D.  The  predicted  scattering  function^^  in  this 
approximation  is  an  anisotropic  lorentzian  in  the  small  q  limit.   The  function 
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Fig.  16.  Change  in  the  line  shape  after  a  quench.   Transverse  scans  of  the 

[0,1,0]  with  a*=1.666R~^   Solid  (dotted)  line  is  a  GS  (LSQ) 
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2 
Inset:  x    (as  later  defined)  of  the  fit  vs.  time  for  the  GS(diamonds)  and 

LSQ(circles). 
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definition  is  presented  in  appendix  B.   This  HT  approach,  was  used  by  Ludwig  et 
al^  ''  to  parameterize  early— time  kinetic  experiments  on  Cu„Au.   The  small  q  limit 
is  a  modified  Lorentzian  where  the  line  width  depends  on  the  direction  in  k-space. 
Although  this  approach  effectively  parameterized  their  data  it  would  not  do  so  here. 
However,  this  is  not  a  contradiction  because  all  of  their  scans  were  taken  in  the  time 
it  took  to  take  our  first  scan,  the  resolution  they  used  was  very  broad  compared  to 
the  present  work,  and  their  sample  was  a  poly— crystalline  wire.  The  modified 
Lorentzian  produced  poorer  fits  to  our  data  than  a  straight  Lorentzian  which  did 
not  fit  well  at  any  time.   However,  if  the  shape  was  convolved  with  the  resolution 
function  the  fits  were  poor  but  acceptable.   It  is  worth  noting  that  the  background 
for  these  fits  was  constrained  to  be  positive,  if  negative  backgrounds  were  permitted 
considerably  better  fits  were  found  in  the  transverse  scans  which  contained 
considerably  few  points  in  the  tails. 

Several  theoretical  predictions  on  coarsening  suggest*^     '     ''  that  the  tails  of 
S(q,t)  should  decay  as  q     ,  consistent  with  a  LSQ  peak.   The  theory  of  Ohta, 
Jasnow,  and  Kawasaki^     ''  (OJK)  derives  an  equation  for  the  line  shape  that 
involves  an  integration  for  each  q.   Using  numerical  integration  non— convolution  fits 
to  this  line  shape  have  been  done  on  selected  peaks  from  quenches  to  various 
temperatures  and  times  such  that  the  resolution  is  small  compared  to  the  overall 
width.  However,  the  resulting  fits  during  coarsening  were  significantly  worse  than 
the  LSQ  fits  for  the  transverse  scans  but  only  slightly  worse  for  the  radial  scans. 
Definitions  for  this  theory  and  the  one  that  follow  are  summarized  in  appendix  B. 

Recently  Z.  W.  Lai  has  developed  a  coarse— grained  model  of  the  phase 
transition  kinetics  for  a  Cu^Au  type  system^     ' .  The  model  is  based  on  a 
Ginzburg— Landau  hamiltonian  with  the  symmetry  of  Cu„Au.   This  model  exhibits 
distinct  regimes  from  early  to  late  time.   Moreover,  the  line  shape  of  the  structure 
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■factor  starts  off  gaussian  and  then  crosses  over  to  approximately  lorentzian— squared 
in  agreement  with  our  observations  especially  in  the  transverse  direction.   At  late 
times  this  model  produces  a  modified  loren^tian  squared  of  the  form, 

^1^=; 272 2-  (4.1) 

^  (l+bjr)^+2b2k^ 

Fits  of  this  form  to  our  data  do  not  favor  this  form  over  that  of  a  pure 

lorenztian-squared  (the  b=0  case).   However,  Z.  Lai  points  out  that  this  variance 

could  be  due  to  the  size  and  time  used  in  the  computations  for  this  model. 

Moreover,  due  to  the  need  for  a  larger  system,  and  hence,  more  computer  time  a 

smaller  anisotropy  in  the  shape  of  the  domains  was  used  in  the  numerical 

calculations  than  is  found  in  Cu«Au.   Z.  Lai  also  suggested  that  this  model  might 

not  be  "completely  settled  into  the  longest  time  asymptotic  regime".   Recalling  that 

during  the  crossover  there  is  a  period  where  the  peak  is  neither  gaussian  or 

lorenztian— squared,  our  fits  were  checked  to  see  if  b2  could  be  related  to  the 

crossover.   However,  there  was  no  evidence  of  b^  contributing  in  any  systematic 

way.  In  fact  even  at  the  crossover  this  modification  has  very  little  effect  on  the 

shape  of  the  fits.   Nevertheless  this  model  exhibits  the  qualitative  behavior  that  we 

have  observed. 

2 
In  table  1  are  listed  the  reduced  chi-squared  (x  )  values  for  non— convolution 

and  some  resolution  convolution  fits  to  gaussian,  lorentzian,  lorentzian— squared, 

lorentzian  to  a  variable  power,  and  the  two  shapes  based  on  the  discussed  theories 

(see  appendix  B  for  functional  definitions).   In  the  fits,  up  to  a  second  order 

2 
polynomial  background  [a(x— x  )  +b(x— x  )+c]  was  added  to  the  model  function. 

Chi— squared  is  the  goodness  of  fit  parameter  used  in  the  fitting  routine 

2  2 

[x  =E(If..-I         )  /(degrees  of  freedom)]   All  but  the  latest  times  given  in  table  1 

are  such  that  the  resolution  width  is  small  compared  to  the  over  all  width. 

Therefore,  the  resolution  is  not  important  for  the  first  four  times  listed  in  table  1. 
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Table  1. 

Chi-  squared  values  for  selected  scans  as  a 

function  of  time  after  a  quench  on  the  bulk  Cu«Au  sample.  Tables  are 

arranged  by  decreasing  quench  temperature  and  direction  in  reciprocal  space. 
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Table  1.  (cont.) 
Values  of  chi-squared 
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95.2 

79.2 

7.6 

7.7 

4.8 

5.6 

5.7 

4.5 

133.2 

121.1 

95.5 

39.1 

40.8 

28.1 


(convolution  fits) 


Gaussian 

Lorent- 
zian 


LSQ 

Lor.to(m) 
(m)  chi 

O.J.K. 
Theory 


0 
1 
2 
0 
1 
2 
0 
1 
2 
0 
1 
2 
0 
1 
2 


20.3 

23.4 

(139.6) 


4.1 

3.5 

2.2 

103.5 

102.1 

89.1 

20.5 

20.8 

20.7 

3.7 

3.3 

2.2 

111.1 

127.0 

90.8 


(2.6) 

2.7 
2.7 


32. 4[ 

85.51 

139. 9j 

31.5 

84.9 

139.2 

24.5 

64.3 

118.0 

186.9 

260.4 

273.4 

185.9 

259.4 

272.9 

161.1 

210.9 

187.9 

6.9 

4.0 

14.3 

6.9 

4.0 

14.4 

5.9 

4.0 

12.4 

2.8 

'2.0] 

3.9 

'1.61 

6.7 

2.5 

;2.oj 

3.9 

1.6] 

6.7 

2.6 

'2.0 

3.9 

'1.6 

6.8 

207.9 

314.0 

265.9 

212.4 

317.9 

271.7 

159.3 

190.4 

182.1 

1.3] 

'1.3 
'1.3 


58.6' 

58.5 

42.4 

24.7 

24.7 

16.7 

15.1 

15.4 

9.4 

1.5 

1.6 

1.6 

23.1 

23.4 

20.1 
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Table  1.  (cont.) 
Values  of  chi-squared 


\or( 
line         \ 
shape       \ 

Traj 
[non 

ler    j 

of 

bg 

0  1 

1 

2 

0 

1 

2 

0 

1 

2 

0 

1 

2 

0 

1 

2 

0 

1 

2  ! 

[isverse                     T=370.8 
convolution  fits) 

1 

1            8          1           23          1 
minutes    [       minutes    j 

48          1 
minutes    j 

123        1 
minutes    j 

494 
minutes 

Gaussian 

(19.3) 

15.1 

(59.7) 

3.5' 

2.8 

2.1 

127.6 

128.4 

114.4 

25.8 

26.1 

17.9 

6.5 

6.2 

5.6 

145.7 

132.9 

126.6 

6.9 

6.7 

6.41 

;2.9) 
3.0 

;3.i) 

29.8] 

29.4 

22.2 

224.7 

227.2 

204.0 

10.9 

11.0 

8.5 

3.8 

3.4 

3.3 

279.4 

257.7 

213.2 

19.9 

20.2 

17.l! 

;2.2) 

2.2 

;2.5) 

47.9| 

45.1 

31.5 

199.5 

190.1 

159.4 

5.3 

4.7 

4.8 

4.8 

3.9 

3.0 

208.9 

203.8 

169.4 

27.0 

24.9 

17.8! 

;2.o] 

2.0 
;2.0j 

24. 9j 

24.9 

22.0 

85.6 

85.0 

76.8 

3.3 

3.1 

3.0 

3.2 

3.0 

2.9 

93.7 

92.0 

81.6 

17.3 

17.3 

14.6| 

;i.8j 

1.8 
;2.0) 

55.0| 

56.2 

41.7 

Lorent- 
zian 

122.7 
125.1 
107.2 

LSQ 

7.7 
7.8 
5.8 

Lor.to(m) 
(m)  chi 

6.9 
7.1 
5.6 

O.J.K. 
Theory 

161.1 
158.6 
119.2 

HT 
Theory 

38.4 
40.2 
26.8| 

(convolution  fits) 


0  1 

1 

3.7| 

31.2| 

50.0] 

26.0] 

57.9] 

Gaussian 

1 

2.9 

30.0 

46.1 

25.5 

57.9 

2 

0 

2.1 
130.1 

22.2 

219.5 

31.6 
184.7 

22.1 
73.9 

42.2 

Lorent- 

37.1 

zian 

1 

128.1 

217.1 

170.5 

71.4 

37.0 

2 

0 

111.1 
24.9 

189.2 
9.1 

136.0 
4.2 

48.4 
3.2 

25.9 

13.2 

LSQ 

1 

25.3 

9.2 

3.5 

3.0 

13.4 

2 

25.8 

8.3 

3.4 

2.7 

8.3 

Lor.to(m) 

0 

(23.7) 

3.3 

'2.81 

2.5 

'2.01 

4.1 

;i-7i 

2.3 

;i.4i 

2.5 

(m)  chi 

1 

28.4) 

2.6 

2.8) 

2.1 

;2.ij 

3.3 

1.7 

2.1  1 

1.3 

2.5 

2 

(209.4) 

2.1 

3.0) 

2.0 

'2.3l 

2.6 

'1.7 

2.1  ( 

1.4 

2.5 

HT 

0 

143.9 

216.3 

209.3 

58.9 

34.8 

Theory 

1 

129.5 

218.1 

163.6 

57.5 

35.6 

2  1 

1 

125.5! 

196.2! 

127.5| 

51.8! 

29.9! 
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Table  1.  (cent.) 
Values  of  chi-squared 


Transverse 
(non  convolution  fits) 


T=365.8°C 


line 
shape 


\order 
\of 
\bg 


1  ° 

Gaussian       1 
2 

Lorentz-        0 

zian            1 

2 

0 

LSQ            1 

2 

Lor.to(m)       0 

(m)  chi           1 

2 

HT               0 

Theory          1 

2 

O.J.K.           0 

Theory          1 

!    2 

8        ! 

23         1 

48          ! 

123        1 

514 

minutes    [ 

minutes    j 

minutes    ] 

minutes    ] 

minutes 

2.8[ 

24.6J 

35.7J 

25.6{ 

60.5| 

2.4 

24.2 

36.2 

25.4 

61.2 

1.8 

17.2 

26.8 

23.0 

40.2 

129.3 

217.1 

187.9 

97.3 

142.0 

131.0 

219.8 

191.6| 

96.9' 

143.2 

112.6 

197.7 

68.2| 

88.6J 

118.3 

24.6 

13.3 

5.31 

3.0J 

8.5 

25.1 

13.5 

5.4 

2.9 

8.5 

16.4 

10.7 

5.2 

2.9 

6.8 

11(20.0) 

5.5 

'3.21 

4.8 

'2.4] 

3.1 

'2.01 

3.0 

;i.9i 

8.3 

(14.7 

5.5 

'3.2) 

4.6 

2.5 

3.1 

'2.0) 

2.8 

1.9 

8.3 

(66.8) 

5.0 

'3.5) 

4.2 

'2.7) 

2.7 

'2.0) 

2.9 

'2.2) 

6.2 

172.6 

235.6 

207.9 

116.3 

158.9 

161.6 

236.8^ 

!12.< 

'1 
1 

103.5 

168.6 

123.0 

193.8 

163.7| 

94.1 

120.6 

5.9 

19.0 

19.3 

18.1 

42.6 

5.9 

19.6 

19.8 

17.7 

43.1 

II 

5.6| 

16.3J 

15.l| 

15.2' 

27.9J 

(convolution  fits) 


0  1 

1 

2.9| 

25.7' 

37.3' 

26.8! 

62.9] 

Gaussian 

1 

2.5 

24.7 

37.0 

26.0 

62.8 

2 

0 

1.8 
132.1 

17.2 
213.2 

26.8 
174.8 

23.1 
84.6 

41.1 

Lorent- 

45.5 

zian 

1 

130.9 

211.2 

174.3 

81.9 

45.0 

2 

0 

109.5 
24.1 

184.4 
11.5 

147.3 
3.5 

.     71.1 
2.9 

27.7 

13.0 

LSQ 

1 

24.5 

11.6 

3.6 

2.7 

13.1 

2 

0 

1(24.1) 

25.7 
2.6 

'3.01 

10.7 
2.8 

'2.3) 

3.6 

2.3 

'1.8] 

2.7 

2.3 

(1.4] 

7.0 

Lor.to(m) 
(m)chi  [ 

3.0 

1 

(27.4) 

2.3 

3.1] 

2.4 

'2.3) 

2.3 

.^•« 

2.2 

(1.4 
(1.4) 

2.9 

2 
0 

(89.0) 

1.8 
148.1 

:3.4) 

2.1 
222.7 

:2.5) 

2.0 
193.9 

1.8 

2.2 
110.7 

2.8 

HT 

40.2 

Theory 

1 

142.5 

221.6 

198.3 

81.4 

40.6 

2  1 

1 

109. 7| 

178. 0| 

148. 8| 

61. 3| 

29.2] 
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Table  1.  (cont.) 
Values  of  chi-squared 


Radial 
(non  convolution  fits) 


T=382.0''C 


\ order 
line         \  of 
shape       \  bg 


1  ° 

Gaussian       1 
2 

Lorent-          0 

zian             1 

2 

0 
LSQ            1 

2 

Lor.to(m)       0 

(m)  chi          1 

2 

HT               0 

Theory          1 

2 

O.J.K 


8 
minutes 


22 

minutes 


47 
minutes 


123 
minutes 


'1.1) 

1.4 

'l.O 


1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
8.7 
8.6 
2.5 
1.6 


(2.4) 
2.4 
2.9 


2.5 
2.5 
1.3 
2.4 
2.4 
1.3 
1.0 
1.0 
1.0 
0.9 
0.9 
0.9 
3.4 
3.4 
3.2 
3.1 


13.31 


(2.4) 
2.4 
3.0 


13.4' 

6.5 

23.0 

23.1 

11.7 

2.3 

2.3 

2.3 

1.7 

1.7 

1.3 

43.3 

45.2 

34.1 

7.6 


38.6 

38.7 

23.7 

103.1 

103.7 

63.5 

6.3 

6.3 

6.4 

3.9 

3.9 

2.2 

178.5 

168.3 

149.2 

19.2 


485 
minutes 


(2.5) 
2.5 
2.9 


23.0] 

23.1 

16.6 

87.9 

88.4 

62.4 

6.4 

6.5 

6.5 

4.1 

4.2 

2.7 

138.8 

139.8 

122.8 

11.4 


(convolution  fits) 


Gaussian 

Lorent- 
zian 


LSQ 

Lor.to(m) 
(m)  chi 

^T 
Theory 

O.J.K. 


0 
1 
2 
0 
1 
2 
0 
1 
2 
0 
1 
2 
0 
1 
2 
0 


'1.0) 

'l.3 

'l.O 


1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
6.6 
5.3 
1.6 


(2.2) 
2.2 
2.6 


2.5 
2.5 
1.3 
2.4 
2.4 
2.3 
0.9 
0.9 
0.9 
0.9 
0.9 
0.9 
3.1 
3.1 
2.6 


I 


(2.0) 
2.0 

2.4 


11.9 

11.9 

5.6 

20.0 

20.0 

16.7 

1.4 

1.4 

1.3 

1.4 

1.4 

1.2 

31.8 

31.6 

19.8 


'1.7) 

^1.7 

'1.9 


28.9 

28.7 

17.1 

38.4 

38.2 

36.5 

2.8 

2.7 

1.4 

1.9 

1.8 

1.3 

133.6 

132.7 

79.0 


I 


13.31 
13.3 

9.4 

8.8 

8.8 

8.8 

4.5 

4.6 

3.1 

2.4 

2.4 

2.2 

208.9 

208.9 

66.3 

9.3 
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Table  1.  (cont.) 
Values  of  chi-squared 


Radial 
(non  convolution  fits) 


T=375.6°C 


\order 
line         \  of 
shape       \  bg 


1  ° 

Gaussian       1 
2 

Lorent-          0 

zian             1 

2 

0 

LSQ            1 

2 

Lor.to(m)       0 

(m)  chi           1 

2 

HT               0 
Theory          1 

2 

O.J.K. 


8 
minutes 


23 

minutes 


49 
minutes 


126 
minutes 


491 
j      minutes 


1.9 
1.8 
1.2 
2.9 
2.9 
1.4 
1.4 
1.4 
1.2 
1.1 
1.1 
1.1 
9.0 
8.7 
4.4 
2.4 


(3.0) 
3.0 
4.0 


8.1' 

8.1 

3.5 

19.5 

19.6 

8.4 

3.3 

3.3 

2.8 

1.7 

1.7 

1.4 

48.7 

47.3 

31.2 

4.0 


(2.7) 
2.7) 
3.2) 


20.3 

20.3 

10.5 

51.3 

51.6 

25.3 

4.5 

4.5 

4.2 

2.2 

2.1 

1.5 

112.8 

111.9 

79.2 

10.2 


(2.6) 
2.6 
3.2 


42.8 

43.0 

24.0 

150.7 

151.7 

95.8 

11.4 

11.5 

11.5 

5.5 

5.5 

2.2 

244.2 

245.7 

209.3 

19.2 


50.1' 

50.4 

36.0 

215.9 

217.4 

153.3 

15.1 

15.2 

15.3 

7.9 

7.9 

4.7 

347.2 

349.2 

279.8 

26.0 


(convolution  fits) 


0  1 

1 

1.9| 

7.5] 

17.4' 

31. 1[ 

27.5] 

Gaussian 

1 

1.8 

7.5 

17.3 

31.1 

27.4 

2 
0 

1.2 
7.2 

4.7 
31.2 

16.4 
51.8 

63.4 
129.6 

149.8 

Lorent- 

21.6 

zian 

1 

7.0 

31.0 

51.6 

129.4 

21.6 

2 

0 

3.8 
1.3 

21.1 
2.3 

38.9 
1.8 

58.5 
3.0 

33.5 

7.4 

LSQ 

1 

1.3 

2.3 

1.8 

3.0 

7.4 

2 

1.2 

2.1 

1.8 

1.6 

4.4 

Lor.to(m) 

0 

'(3.41 

1.1 

'2.71 

1.7 

'2.21 

1.7 

1.91 

2.8 

'1.51 

3.4 

(m)  chi 

1 

(3.4) 

1.1 

'2.71 

1.6 

'2.2] 

1.6 

1.9) 

2.8 

1-5 

3.4 

2 
0 

(4.7) 

1.1 
7.2 

'3.5 

1.4 
30.7 

:2.5) 

1.3 
61.0 

[2.2) 

1.3 
149.9 

'1.6 

3.0 

HT 

463.8 

Theory 

1 

6.9 

30.4 

59.3 

150.7 

466.8 

2 

0  1 

1 

4.0 

1 

22.2 

1 

47.2 

1 

80.3 

1 

270.6 

O.J.K. 

20.91 
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Table  1.  (cont.) 
Values  of  chi-squared 


Radial 
(non  convolution  fits) 


T=370.8''C 


\ order 
line        \  of 
shape       \  bg 


Gaussian 

Lorent- 
zian 


LSQ 

Lor.to(m) 
(m)  chi 

Theory 
O.J.K. 


0 
1 
2 
0 
1 
2 
0 
1 
2 
0 
1 
2 
0 
1 
2 
0 


8         1 

23 

1 
1 

49         1 

126        1 

487 

minutes    ] 

minutes    j 

minutes    [ 

minutes    ] 

minutes 

1.9| 

8.8j 

19.4] 

49.1| 

51. 7| 

1.9 

8.7 

19.5 

49.4 

52.0 

1.0 

3.1 

9.1 

28.7 

37.2 

3.6 

23.6 

65.4 

170.0 

237.7 

3.7 

23.7 

65.8 

171.1 

239.3 

1.5 

8.9 

30.5 

100.3 

164.9 

1.5 

3.8 

7.5 

12.6 

17.3 

1.5 

3.8 

7.6 

12.6 

17.4 

1.1 

2.8 

6.3 

12.4 

17.4 

II  (3.91 

1.0 

'3.21 

1.5 

'3.0) 

2.5 

;2.7] 

5.4 

'2.7) 

7.9 

(4.0) 

1.0 

'3.2) 

1.4 

'3.0) 

2.5 

'2.7 

5.5 

2.7 

7.9 

(8.9) 

0.9 
17.0 
16.6 

6.3 

:4.5) 

1.1 
69.7 
68.4 
40.2 

;3.8) 

1.7 
144.6 
152.1 
101.6 

[3.2] 

2.8 
317.3 
320.2 
249.2 

3.1 

4.9 
362.8 
366.5 
304.8 

II 

2.5| 

4.0l 

8.9| 

21.7| 

22.7| 

(convolution  fits) 


Gaussian 


Lorent- 
zian 


LSQ 

Lor.to(m) 
(m)  chi 

TT 
Theory 

O.J.K. 


0 
1 
2 
0 
1 
2 
0 
1 
2 
0 
1 
2 
0 
1 
2 
0 


(3.7) 
3.7 
8.4 


1.9 

1.9 

1.0 

14.0 

14.0 

5.5 

1.4 

1.4 

1.4 

1.0 

1.0 

0.9 

13.8 

13.6 

5.7 


(2.91 
2.9 
4.0 


8.3 

8.1 

4.2 

47.3 

47.0 

28.9 

2.6 

2.6 

2.4 

1.4 

1.3 

1.1 

47.0 

46.8 

29.5 


(2.5) 
2.5 
3.2 


16.9 

16.8 

14.7 

81 

81 

58 

3 

3 

3 

2 

2 

1.5 
88.7 
86.3 
66.0 


(2.0) 
2.0 
2.3 


37.0' 

29.3| 

37.0 

29.3 

20.6 

20.4 

101.5 

68.7 

101.4 

68.7 

87.2 

41.4 

2.8 

5.6 

2.8 

5.6 

2.1 

3.4 

2.8 

;i.6) 

3.1 

2.8 

1.6) 

3.1 

1.6 

'1.7 

2.7 

165.3 

442.7 

165.5 

449.2 

109.5 

1 

64.3 
17.1 
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Table  1.  (cont.) 
Values  of  chi-squared 


Radial 
(non  convolution  fits) 


T=365.8"C 


\  order  jj 
line  \  of  " 
shape       \  bg 


1  ° 

Gaussian       1 
2 

Lorent-          0 

zian             1 

2 

0 

LSQ            1 

2 

Lor.to(m)       0 

(m)chi          1 

1    2 

HT               0 

Theory          1 

!    2 

O.J.K 


16         I  27         I 

minutes    |       minutes    ] 


52  I  128        I         496 

minutes    j       minutes    |       minutes 


(3.5) 
3.5 
4.6 


5.0' 

4.9 

1.8 

12.2 

12.2 

3.6 

2.6 

2.5 

1.6 

1.1 

1.0 

0.9 

45.4 

43.3 

21.6 

3.7 


(3.21 
3.2 
4.2 


10.5' 

10.5 

3.8 

26.4 

26.6 

9.6 

4.1 

4.1 

3.0 

1.6 

1.5 

1.2 

80.0 

79.9 

45.1 

5.1 


(2.9) 
2.9 
3.6 


24.1! 

24.3 

10.9 

67.4 

67.9 

31.1 

7.5 

7.6 

6.5 

3.0 

3.0 

2.1 

167.0 

168.1 

110.9 

11.7 


(2.7) 
2.7 
3.2 


51.6' 

51.9 

29.7 

170.0 

171.2 

96.1 

12.6 

12.7 

12.2 

5.3 

5.3 

2.9 

322.0 

323.1 

252.5 

23.2 


55.0' 

55.4 

37.9 

237.7 

239.3 

161.9 

16.3 

16.5 

16.5 

7.7 

7.8 

4.3 

364.0 

366.5 

302.7 

23.9 


(convolution  fits) 


Gaussian 

Lorent- 
zian 


LSQ 

Lor.to(m) 
(m)  chi 

~HT 
Theory 

O.J.K. 


0 
1 
2 
0 
1 
2 
0 
1 
2 
0 
1 
2 
0 
1 
2 
0 


(3.2) 
3.2 
4.2 


4.8| 

9.9] 

21.4] 

40.2] 

32.7] 

4.7 

9.8 

21.4 
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One  would  expect  for  late  times  the  convolution  of  the  resolution  function  with  a 
LSQ  would  fit  better  than  a  straight  LSQ,  as  is  the  case  in  the  radial  data. 
However,  this  is  so  for  the  transverse  data,  implying  that  the  uncertainty  in  the 
transverse  resolution  function  is  too  great  to  adequately  define  its  shape,  which  is 
the  case  due  to  the  uncertainty  in  the  shape  of  the  mosaic. 

The  difference  in  the  chi  squared  values  between  the  LSQ  and  lorentzian  to 
the  m  power  during  coarsening  is  never  so  great  as  to  make  one  believe  that  the 
peak  is  not  a  LSQ  shape.  In  fact,  the  dependence  of  chi  squared  on  m  is  relatively 
flat  between  about  1.3  to  over  3  but  rises  quickly  near  1.0. 

Fits  to  the  convolution  of  the  resolution  function  and  the  first  theoretically 
based  shape  discussed  have  not  been  done  in  the  transverse  direction.   It  is  unlikely 
that  such  a  fit  would  be  an  improvement  over  the  non— convolution  fit,  because  this 
shape  is  similar  to  a  LSQ  in  the  tails  and  the  convolution  fitting  with  a  LSQ 
function  did  not  improve  the  fits  in  the  transverse  direction.   However,  despite  the 
computer  time  needed  for  such  a  convolution  fit  a  few  have  been  done  for  the  radial 
data.  Improvement  in  the  chi  squared  values  in  the  convolution  fits  relative  to  the 
non— convolution  fits  for  this  first  theoretical  shape  were  similar  to  what  was  found 
in  the  LSQ  shape.  The  chi  squared  values  for  the  LSQ  function  at  late  times  are 
approximately  half  those  of  the  OJK  theory.  Keeping  in  mind  that  the  justification 
for  assuming  a  LSQ  peak  has  always  been  that  the  tails  decay  as  q      as  OJK  theory 
predicts,  it  is  some  what  surprising  that  a  LSQ  shape  actually  fits  better  than  the 
full  theory.   However,  several  approximations  were  made  in  the  theory  that  may 
explain  this.   The  theory  assumes  an  isotropic  system,  and  a  gaussian  distribution  of 

fluctuations  of  the  domain  wall  positions.   The  authors  point  out  the  dependence  of 

f35) 
their  results  on  the  choice  of  this  distribution>     ' 
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Another  way  to  see  the  change  in  the  line  shape,  and  hence,  S(q,t)  is  to  scale 
the  peaks.  If  the  peaks  were  the  same  shape,  plots  of  counts  divided  by  the  fitted 
amplitude  vs.  position  divided  by  the  width,  would  lie  on  the  same  curve.  Figure  17 
shows  some  of  the  raw  data  for  the  transverse  365.8  C  quench  with  LSQ  fits  to  all 
but  the  earliest  one  which  is  fit  to  a  GS.  Figure  18  shows  the  data  after  being  scaled, 
the  two  lines  are  a  scaled  LSQ  fit  and  a  scaled  GS  fit.  From  the  scaled  data  it  is 
apparent  that  the  earliest  peak  has  a  different  shape  than  all  the  rest.   This 
demonstrates  that  S(q,t)  can  not  be  scaled  from  the  coarsening  regime  to  the  growth 
regime,  but  does  scale  during  coarsening.   Moreover,  S(q,t)  should  not  scale  in  this 
way  during  growth  when  the  ratio  of  the  volume  of  the  sample  in  one  of  the  ordered 
phases  to  the  volume  in  a  disordered  phase  is  constantly  changing.   In  order  for 
scaling  to  be  valid  the  physical  picture  after  rescaling  must  be  identical  to  an  earlier 
time.   This  is  not  the  case  during  the  growth  stage,  if  the  length  scale  is  reduced 
during  growth  the  percentage  of  the  sample  that  is  in  a  disordered  state  wiU  be 
greater  than  it  was  at  all  earlier  times. 

The  crossover  in  the  line  shape  has  been  observed  in  over  a  dozen  quenches. 
Most  of  these  quenches  were  to  375°C  (T^-10°C).  In  one  quench  to  378C  (T^-7°C) 
the  crossover  is  present,  however,  the  greatest  difference  in  their  respective 
chi— squared  values  while  the  peak  favored  a  GS  line  shape  was  1.37  to  2.21  where 
for  deeper  quenches  ratios  of  1  to  9  where  found.   Moreover,  in  the  382.0  C  quench 
the  first  peak  mildly  favored  a  GS  line  shape  with  chi-squared  values  of  1.9  and  2.2 
for  GS  and  LSQ  fits  respectively.   In  the  deeper  quenches  to  below  376  C  where  the 
crossover  is  pronounced  the  peak  does  not  fit  well  to  either  a  GS  or  LSQ  shape  at 
the  point  the  chi-squared  values  for  the  GS  and  LSQ  fits  cross.  In  fact,  the  peak 
starts  off  GS  and  then  continuously  changes  until  it  is  LSQ,  but  during  the  change 
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FIG.  17.  Raw  data  on  the  [0,1,0]  bulk  superlattice  diffraction  peak.   Solid 
lines  are  fits  to  the  data.   For  the  7.60  min.  scan  a  gaussian  fit  is 
shown,  the  rest  are  lorentzian—  squared  fits. 
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FIG.  18.  Data  shown  in  Cg.  17  scaled  by  amplitude  and  width.   The  solid 

Lines  describe  a  gaussian  and  a  lorentzian—  squared  shape.  The  7.60  min. 
peak  is  gaussian,  the  rest  are  lorentzian—  squared. 
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it  is  neither.   The  rate  of  change  seems  to  be  slower  with  lower  temperature,  in  fact 
for  quenches  above  376°C  it  seems  as  if  the  change  to  LSQ  was  almost  complete  by 
the  end  of  the  first  scan.   In  a  typical  quench  300  second  scans  were  taken  at  the 
start  of  the  experiment,  attempts  with  shorter  scanning  times  only  made  the  effect 
less  clear,  due  to  the  shape  of  the  peak  being  statistically  less  well  defined. 
Moreover,  for  quenches  close  to  T^,  the  intensity  of  the  peak  is  weaker  than  it  is  for 
deeper  quenches,  making  it  more  difficult  to  see  this  change  at  higher  temperatures 
for  the  same  reason  as  above. 

The  change  in  line  shape  also  occurred  earlier  with  each  quench.   This 
suggests  that  the  coarsening  stage  is  reached  sooner  if  the  sample  has  been  ordered 
for  long  times  and  disordered  at  temperatures  within  (T  +40°C)  prior  to  the 
quench.   If  the  sample  is  heated  to  700°C  for  10  or  more  hours  the  change  in  line 
shape  is  again  found  to  take  place  over  the  original  longer  time  period.   Coarsening 
being  reached  sooner  suggest  that  nucleation  and  or  growth  is  faster  in  samples  that 
have  not  been  annealed  at  high  temperatures  than  in  ones  that  have.   This  effect 
may  involve  the  migration  of  vacancies  to  domain  wall  and  grain  boundaries  during 
coarsening. 

The  change  in  line  shape  although  present  in  the  radial  data  is  not  as  clear  as 
in  the  transverse  data.   There  are  two  reasons  for  this.   The  first  is  that  the  radial 
data  was  taken  after  the  transverse  data  and  the  history  dependence  discussed  above 
would  lessen  the  effect.   The  second  is  that  the  width  in  the  radial  direction  is 
narrower  being  due  to  only  type  2  walls.   Therefore,  in  the  radial  direction 
resolution  is  more  important.   Type  2  walls  being  farther  apart  may  also  have  an 
effect  on  the  peak  shape. 

Figure  19  displays  the  integrated  intensities  (fitted  amplitude  multiplied  by 
fitted  width)  versus  time  after  a  thermal  quench  for  the  transverse  directions  at 
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FIG.  19.  Intensity  (width  of  convolution  fits  multiplied  by  the  amplitude  of 

the  fit)  vs.  time.   The  380. 0°C  data  shows  a  delay  of  approximately  1900 
seconds  that  corresponds  to  a  period  of  nucleation.   The  data  shown  was 
taken  in  the  transverse  direction  at  the  [0,1,0]. 
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several  temperatures.   The  380.0°C  transverse  data  shows  a  delay  immediately  after 
the  quench,  similar  delays  have  also  been  observed  in  other  quenches  to 
temperatures  closer  to  T^.   This  delay  is  a  result  of  an  incubation  time  for 
nucleation.   Such  delays  have  been  clearly  observed  in  single  grains  of  Mg  In  by 
Konishi  and  Noda^^^^).   The  delay  time  depends  on  the  sample  history  which  for 
the  bulk  sample  is  shown  in  figure  20.  In  Ught  of  the  history  dependence  discussed 
above  it  is  not  too  surprising  that  the  deepest  quench  in  which  a  delay  is  seen  is  the 
first  one  after  etching  and  annealing  at  700°C.   After  the  delay  a  period  of  rapid 
growth  in  intensity  is  observed  corresponding  to  the  Unear  rise  with  time.   At  late 
times  the  intensity  starts  to  saturate,  asymptotically  approaching  some  maximum 
(I^(T))  that  increases  with  decreasing  temperature.   The  intensity  after  the  delay 
fits  well  to 

I(t)  =  I^[l-exp[-[(t-t^)/r]]].  (4.2) 

Figure  21  shows  the  data  from  19  scaled  by  plotting  I/I^  vs.  t/6  where  b  is  the  time 
for  the  intensity  to  reach  1^/2.   As  expected,^^^)  the  data  faUs  on  a  universal  curve. 
In  theory,  one  could  use  r  from  the  fits  to  equation  4.2,  however,  that  proved  to  be 
less  accurate.  The  reason  for  this  is  simply  that  the  r  in  the  fits  corresponds  to  a 
crossover  time,  which  is  early  and  our  time  resolution,  being  at  best  1/2  the  time  to 
take  a  spectrum,  is  a  sizable  percentage  of  r. 

The  delay  seen  in  the  integrated  intensity  is  also  in  the  amplitude.   A  plot  of 
amplitude  vs.  time  is  shown  in  figure  22.   The  amplitude  used  in  the  plot  is  the 
convolution  fitted  amplitude  and  as  such  is  corrected  for  resolution  effects  in  the 
direction  of  the  scan.   The  amplitude  is  proportional  to  the  square  of  the  order 
parameter  and  should  scale  as  L(t)^..  Therefore,  the  amplitude  should  have  a  time 
dependence  t^/^-   However,  the  data  has  only  been  corrected  for  resolution  effects  in 
one  direction  and  the  measurement  approximately  integrates  over  the  other  two 
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FIG.  20.  Bulk  Cu^Au  sample  history  starting  from  etching  and  annealing  in  an 

Ar  atmosphere  at  700°C  to  end  of  radial  quenches.    Xn  etching  and 
annealing  followed  by  quenches  similar  to  above  was  done  prior  to 
what  is  shown. 
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FIG.  21.  Scaled  intensity  vs.  scaled  time.  The  nucleation  period  for  the  380. OC 
data  was  subtracted  from  the  time  before  scaling.  The  unsealed  data  is 
shown  in  fig.  19. 
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FIG.  22.  Convolution  fitted  amplitude  vs.  time  for  the  temperatures 
shown.  The  delay  seen  in  the  plot  of  the  integrated  intensity 
is  also  seen  here. 


76 


J n 

dimensions,  therefore,  the  measured  amplitude  should  roughly  scale  as  L(t)        and 

1/2 
have  a  time  dependence  of  t  '   .   A  plot  of  Log,Q(amplitude)  vs.  Log,^(time)  is 

shown  in  figure  23.   The  late  time  slope  for  the  deeper  quenches  is  0.46,  a  value  close 

to  the  value  found  for  the  time  dependence  of  L  as  will  shortly  be  discussed.   Trying 

to  extract  the  time  dependence  of  L  from  the  amplitude  would  not  be  a  good 

method  due  to  the  uncertainty  of  the  resolution  effect.   However,  it  is  reassuring 

that  a  value  close  to  the  expected  value  of  1/2  is  found,  indicating  that  the  scaling 

involved  is  as  expected  from  theory.   The  curvature  in  the  Log— Log  plot  is  at  least 

in  part  due  to  the  resolution  function. 

From  a  simplistic  Ising  model  calculation  it  might  be  argued  that  the 

integrated  intensity  should  be  constant  and  that  in  the  previous  graphs  the  rise  in 

integrated  intensity  is  due  to  the  resolution  function  integrating  over  more  of  the 

peak  as  the  peak  becomes  narrower.   This  is  not  the  case.   First  of  all  the  integrated 

intensity  is  calculated  from  the  product  of  the  fitted  width  and  the  fitted  amplitude, 

not  from  the  sum  of  all  the  counts.   Hence,  the  intensity  in  the  tails  is  accounted 

for(at  least  in  the  direction  scanned).   Figure  24  shows  the  integrated  intensity, 

amplitude  and  constant  background  as  a  function  of  time.   In  this  multiple  plot  the 

background  has  been  multiplied  by  1000  to  plot  it  on  the  same  scale  as  the 

amplitude.   The  fact  the  fitted  background  is  1000  times  or  more  smaller  than  the 

fitted  amplitude  implies  that  the  missing  intensity  would  have  to  be  spread  out  in  a 

volume  10    times  larger  than  what  was  effectively  covered.   It  is  obvious  that  some 

of  the  intensity  is  being  falsely  counted  as  the  background  due  to  the  rise  in  the 

background  with  time.   However,  this  change  in  fitted  background  is  in  the  wrong 

direction  because  the  argument  says  intensity  is  missing  at  early  times  not  late. 

This  is  a  different  problem  that  is  discussed  in  connection  with  power  law  growth. 
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FIG.  23.  Log,Q(amplitude)  vs.  Log,Q(time)  for  the 
same  data  at  in  figure  22. 
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FIG.  24.  Plots  of  integrated  intensity,  amplitude  and  fitted  constant 

3 
background.   The  background  has  been  multiplied  by  10    in  order 
to  plot  it  on  the  same  scale  as  the  amplitude.  The  integrated 
intensity  has  also  be  scaled  to  fit  on  the  graph. 
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The  background  should  probably  be  quadratic  due  to  a  diffuse  peak  arising  from 
thermal  fluctuations.  Moreover,  the  fitted  background  for  the  second  and  third 
points  is  slightly  negative,  once  again  indicating  that  if  any  thing  the  change  in 
intensity  is  greater  than  suggested  by  the  plots. 

The  above  argument  shows  that  any  missing  intensity  can  not  be  coming 
from  the  direction  that  was  scanned.  The  resolution  function  in  the  vertical 
direction  is  large  which  strongly  implies  that  any  missing  intensity  is  lost  in  the 
radial  direction.   As  it  happens  several  quenches  with  radial  scans  using  the  graphite 
monochromator  and  analyzer  (what  has  been  called  the  transverse  set  up)  were 
tried.   The  resolution  in  that  direction  turned  out  to  be  too  large;  within  an  hour 
after  a  quench  the  instrument  was  resolution  limited.   The  width  after  six  hours 
suggests  that  the  resolution  was  close  to  .02^"  (FWHM).   Remember  that  at  the 
[0,1,0]  the  radial  direction  probes  the  thin  part  of  the  scattering  disk.   Based  on  this 
value  for  the  radial  resolution  the  intrinsic  peak  width  is  equal  to  the  instrumental 
resolution  width  16  minutes  after  the  quench.  At  this  point  the  peak  should  be 
almost  completely  integrated  over  (for  the  deeper  quenches).   Some  of  the  intensity 
could  still  be  lost  in  the  tails  due  to  differences  in  line  shapes.  The  intrinsic  line 
shape  at  this  time  is  approaching  that  of  a  LSQ  and  the  instrumental  resolution 
(figure  11)  is  some  thing  between  a  gaussian  and  a  LSQ.   A  time  of  16  minutes 
corresponds  to  the  fourth  point  in  figure  24,  which  has  an  intensity  of  35  percent  of 
the  final  intensity.   Hence,  about  65  percent  of  the  change  in  intensity  is  real.   An 
even  stronger  argument  can  be  made  for  the  radial  data.   A  change  in  integrated 
intensity  has  also  been  seen  in  other  systems  by  N.  Wakabayashi''      -' 

Figure  25  shows  a  plot  of  log,Q(L)  vs.  logjQ(t)  whose  slope  is  "a"  in; 

L(t)  «  t^  (4.3) 

In  table  2  are  listed  the  exponents  "a"  for  curvature  driven  growth,  as  a  function  of 
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FIG.  25.  Log,  ^(Length)  vs.  Log,p,(time  after  the  quench).  L  <x  2T/r  ,  where  F  is 

the  LSQ  fitted  width.  Solid  lines  are  linear  fits  whose  slope  corresponds 
to  the  exponent  "a"  for  power  law  growth. 
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temperature  and  type  of  background.   For  the  radial  data  the  order  of  the 
polynomial  background  was  not  important  because  only  the  constant  term  gave  any 
sizable  contribution.  The  exponent  for  the  radial  data  was  found  to  be  .50±.03  and 
temperature  independent.   Here,  the  error  in  "a"  was  estimated  from  the  effect  of 
changing  the  resolution  correction,  as  well  as  the  standard  deviation  as  calculated 
from  values  at  different  temperatures,  which  was  0.014.   An  error  in  the  resolution 
correction  would  contribute  a  systematic  error,  however,  this  correction  as  discussed 
earlier  is  dominated  by  the  instrumental  resolution,  which  is  well  known.   This  is 
not  the  case  in  the  transverse  direction,  where  the  resolution  correction  is  less  well 
defined  and  the  order  of  the  polynomial  background  seems  to  affect  the  result. 

Table  2.  Exponents  "a"  for  power  law  growth  (equation  4.3). 


direction     ] 
background 


temperature 


382. 0°C 
380. 0°C 
375. 6°C 
370.8°C 
365. 8°C 


average 
std.  dev. 


constant 


transverse 
linear 


"a" 


0.45 
0.54 
0.40 
0.42 
0.40 


0.44 
0.06 


"a" 


0.47 
0.56 
0.40 
0.43 
0.43 


0.46 
0.06 


][!  radial 


quadratic     {constant 


"a" 


I  * 


0.49 
0.58 
0.45 
0.46 
0.46 


0.49 

0.51 
0.51 
0.49 


0.49 
0.05 


0.50 
0.01 


It  should  be  noted  that  the  380°C  transverse  data  is  inconsistent  with  all  the 
other  data,  yielding  a  significantly  larger  exponent  than  the  others.   The  reason  for 
this  is  unknown.   However,  this  set  of  data  was  from  the  first  quench  after  etching 
and  annealing  of  the  sample,  and  hence,  can  be  separated  temporally  from  the  rest 
of  the  data.   During  the  experiment  the  filament  in  the  x— ray  source  burnt  out  and 
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was  replaced,  however,  this  happened  late  in  the  series  of  quenches  rather  than  after 
the  first  one.   However,  changes  in  the  filament  before  failing  may  be  important. 
Evidence  for  the  change  in  "a"  being  due  to  the  filament  can  be  found  in  the  fact 
the  exponents  "a"  can  be  correlated  to  the  order  that  the  data  was  collected.   The 
value  obtained  for  "a"  in  the  transverse  direction  decreases  with  time  up  till  the 
filament  was  replaced  then  goes  up.  This  could  be  a  result  of  the  filament  sagging 
with  time,  and  hence,  changing  the  effective  resolution.   Movement  of  the  x— ray 
beam  would  have  at  least  two  effects,  one  it  would  very  slightly  change  the  area  of 
the  sample  considered,  perhaps  having  an  effect  through  the  mosaic,  two  it  would 
misalign  the  sample.   However,  the  most  likely  way  that  the  resolution  would  be 
changed  is  simply  by  increasing  the  spot  size  on  the  anode  that  the  x— rays  are 
emitted  from.   A  lower  value  for  "a"  would  be  consistent  with  an  under  correction 
for  resolution.   Furthermore,  this  change  in  "a"  carmot  be  correlated  with 
temperature.   The  filament  was  replaced  before  the  382°C  quench  and  after  all  the 
others  listed  in  table  1.   Consequently,  the  instrumental  resolution,  which  was 
measured  after  the  series  is  better  known  for  the  382°C  data.   As  can  be  seen  in 
table  one  the  exponents  "a"  in  the  transverse  direction  found  from  assuming  a 
constant  or  linear  background  are  smaller  than  those  from  a  quadratic  background. 
All  values  of  "a"  in  the  transverse  direction  are  less  than  0.5  except  for  those 
connected  with  the  380  C  quench,  which  are  all  above  0.5.   In  any  case,  the  average 
values  for  each  type  of  background  are  within  their  respective  standard  deviations  to 
each  other  and  0.5  in  spite  of  the  systematic  errors  in  determining  "a"  from  the 
transverse  scans.   More  importantly  it  was  type  2  domain  walls  that  were  suspected 
of  not  following  curvature  driven  growth.   Type  2  walls  are  probed  by  radial  scans, 
and  hence,  "a"  is  known  to  a  high  degree  of  accuracy  for  the  domain  walls  in 
question.   This  implies  that  curvature— driven  growth  applies  to  different  types  of 
walls  and  confirms  scaling  in  an  anisotropic  system.^     ^ 
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The  possibility  of  a  lower  exponent  for  the  transverse  data  might  be 

explained  as  being  due  to  impurities  or  incorrect  stoichiometry.   Random 

f 81  8"^  84  1  C\R\ 
imperfections  lead  to  slower  growth,  possible  a  log(t)  growth  lawA     '     '     '      -* 

However,  a  systematic  error  in  the  resolution  correction  is  the  more  likely 

explanation,  especially  in  light  of  the  fact  type  2  walls  are  not  affected. 

The  Experiments  on  Sputtered  Films  and  the  Effect  of  Stoichiometry 

Experimental  Details 

A  single  target  of  Cu^Au  was  used  to  produce  the  fine  grained  sputtered 
films  on  substrates  of  0.13  mm  thick  HN  kapton  or  0.13  mm  thick  ALO^.   A  glass 
slide  was  placed  in  the  chamber  with  the  other  substrates  so  that  film  thickness 
could  be  determined  with  a  stylus  profilometer.   Electron— microprobe  techniques 
verified  the  stoichiometry  and  Rutherford  back  scattering  showed  the  composition 
to  be  uniform  as  a  function  of  depth.   The  film  characteristics  are  independent  of 
which  substrates  was  used.   The  stoichiometric  film  (film  A)  was  deposited  on 
kapton  to  a  thickness  of  7500  ±300  angstroms,  with  a  composition 
Cuq  -r^  ni^^  2'5±  01  ^  determined  by  electron  microprobe  techniques. 

The  nonstoichiometric  film  (film  B)  was  deposited  to  an  initial  thickness  of 
10,000+\^00  angstroms  on  an  ALO„  substrate.   Subsequently  the  film  was  placed 
in  contact  with  a  clean  Cu  block  and  annealed  under  vacuum  for  2.5  hours  at  500  C. 
The  clean  Cu  block  was  the  heating  rod  previously  discussed  in  connection  with  the 
furnace.  The  tantalum  foil  used  in  later  experiments  to  shield  the  sample  was  added 
due  to  experience  with  this  film.   After  the  data  was  collected  and  the  transition 
temperature  was  found  to  be  299  C  the  film  was  recharacterized  by 
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electron-microprobe  techniques,  the  results  showed  the  sample  to  have  a 

homogeneous  concentration  profile  of  Cu^  ^„        An  Ti,io  cf^i^u- 

y      ""  "^^       0.79±.01      0.21±.0r      ^^  stoichiometry  is 

consistent  with  the  transition  temperature. 

Preliminary  x-ray  characterization  of  the  films  found  the  [1,1,1]  direction  to 
be  normal  to  the  substrate  .0  within  a  mosaic  spread  of  approximately  0.2  X"^  or 
3°  half  width  at  half  maximum  (HWHM).   The  in  plane  orientations  were  random. 
The  lattice  constant  at  room  temperature  was  3.74  A.   Prior  to  the  first  quenches 
both  films  were  annealed  at  350°C  for  over  12  hours  in  a  vacuum.   EarHer 
measurements  on  other  films  showed  changes  in  the  relative  intensities  of  the  [1,1,1] 
and  in  plane  fundamentals  as  a  result  of  the  first  annealing,  implying  an 
improvement  in  the  orientation  of  the  grains  along  the  [1,1,1]  direction. 

Inirially  kapton  substrates  were  used.  Kapton  is  transparent  to  x-rays,  easily 
cut  to  fit  any  sample  mount,  and  inexpensive.   However,  it  bows  when  heated  to 
high  temperatures.   Heating  the  kapton  prior  to  film  deposition  may  have  helped, 
but  it  must  be  heated  in  a  vacuum.   Even  in  a  vacuum  of  a  few  millitorr  the  kapton 
would  deteriorate  over  several  weeks  if  held  above  350C,  as  was  done  during  the 
experiments.   Another  problem  with  kapton  substrates  is  its  low  thermal 
conductivity,  most  of  the  thermal  conduction  was  in  the  film  itself.  Substrates  of 
AI2O3  have  better  thermal  properties  than  kapton  for  these  experiments. 

Kapton  was  also  used  for  the  x-ray  windows  on  the  first  version  of  the 
furnace  which  was  used  in  connection  with  all  the  sputtered  films.   Although  the 
vacuum  obtained  with  kapton  windows  is  a  little  poorer  than  with  beryllium 
windows,  they  are  considerably  less  costly  and  worked  very  well  in  the  testing  of  the 
design  of  the  furnace,  prior  to  manufacturing  the  present  version.   Moreover,  kapton 
unlike  beryllium  is  reasonably  transparent  to  visible  light  allowing  visual 
observations  of  the  interior  of  the  furnace  which  helped  pinpoint  weakness  in  the 
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design.   The  only  other  relevant  difference  between  the  furnace  described  for  the 
bulk  sample  and  this  one  is  that  the  sample  mount  was  a  frame  not  a  plate,  allowing 
the  diffraction  peaks  to  be  observed  in  a  transmission  geometry.   Improvements  in 
the  heating  element  are  discussed  in  connection  with  the  experiments  on  the  M.B.E. 
films  where  the  changes  where  implemented. 

The  data  was  collected  in  the  radial  direction  with  a  PSD  as  described  earlier 
and  shown  schematically  in  figure  10a.   The  instrumental  resolution  was  .06° 
(HWHM  in  29)  or  0.0042A~-^  at  the  [1,0,0],  and  0.004lA~-^  at  the  [1,-1,0]. 
Quenches  on  both  the  [1,0,0]  and  [1,-1,0]  were  done,  in  order  to  observe  both  types 
of  walls.   Radial  scans  widths  of  the  [1,0,0]  are  broadened  by  type  2  walls.   The 
[1,0,-1]  peak  width  is  affected  by  both  types  of  walls  but  is  dominated  by  type  1 
walls.  The  notation  [1,-1,0]  instead  of  [1,1,0]  is  used  to  indicate  that  it  is 
perpendicular  to  the  [1,1,1],  or  in  the  plane  of  the  film.   The  transition  temperature 
T   for  each  sample  was  determined  from  the  disappearance  of  superlattice  peaks 
upon  heating.   The  samples  were  heated  well  above  T    to  ensure  complete  disorder, 
then  cooled  to  an  initial  temperature  T-  >  T  .   After  establishing  equilibrium  the 
films  were  quenched  rapidly  to  T,  <  T      X-ray  scattering  scans  were  collected  and 
stored  sequentially  at  fixed  Tr. 

Analysis  for  Data  on  the  Sputtered  Films 

In  analyzing  the  data  a  gaussian  shape  was  assumed  for  all  the  Otted  peaks. 
Superlattice  peaks  from  a  bulk  sample  are  lorentzian  squared  like,  not  gaussian. 
The  shape  of  the  peak  due  to  the  existence  of  anti— phase  domains  is  difficult  to  see 
experimentally  in  films  due  to  the  contribution  from  finite  size.   A  gaussian  shape, 
however,  can  be  supported.   Each  grain  in  the  film  that  contributes  to  the 
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scattering,  will  in  general  be  out  of  phase  with  the  rest,  producing  the  superposition 
of  (possibly  LSQ  shaped)  peaks  with  a  random  distribution  in  phase.   The  resulting 
peak  should  be  well  approximated  by  a  gaussian.   Furthermore,  the  peaks  as 
measured  where  gaussian.  Therefore,  the  data  was  fit  to  gaussians,  and  then 
instrumental  resolution,  finite  size  and  strain  were  subtracted  in  quadrature. 


=  v^m-^  (4-4) 


'fit 

where  (t   is  the  correction  from  instrumental  resolution,  finite  size  and  strain.   The 

values  of  a   are  listed  with  the  results.   Instrumental  resolution  was  measured  by 
c  ■' 

the  width  of  a  Si  [1,1,1]  as  was  done  with  Ge  for  the  bulk  data.  The  contribution 
from  finite  size  and  strain  was  estimated  from  measurements  of  fundamental  Bragg 
peaks  in  the  same  directions  in  reciprocal  space  as  the  superlattice  peak.   The 
measurements  from  the  fundamental  peaks  were  not  sufficient  to  consistently 
determine  the  corrections  for  these  samples.   The  actual  correction  used  was  a 
probable  upper  bound,  which  as  discussed  below,  may  cause  a  cause  a  systematic 
error  in  analyzing  the  data. 

The  time  dependence  is  analyzed  by  fitting  to  functional  forms  as  discussed 
in  chapter  four,  specifically  power  law  (equations  2.5  and  4.3)  and  logarithmic 
(equation  2.12)  time  dependence.   The  comparisons  made  between  samples  are  such 
that  the  quench  depths  are  the  same.  The  quench  depth  being  the  same  means  that 
the  absolute  temperatures  are  187  C  apart.  In  order  to  rule  out  the  possibility  of 
any  effect  being  due  to  the  difference  in  absolute  temperature  one  quench  to  278  C 
on  the  bulk  sample  was  done.   The  results  indicated  curvature  driven  growth  as  is 
found  at  higher  temperatures. 
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Results  from  the  Sputtered  Films 

Some  typical  raw  data  from  a  time  resolved  run  on  the  Cu  rich  film  are 
shown  in  figure  26.   The  data  on  the  stoichiometric  film  is  similar  but  with  better 
signal  to  background.   Figure  27  shows  log— log  plots  of  L  vs  t  for  the  [1,0,0]  and 
[1,-1,0]  superlattice  peaks  in  both  film  A  and  B.   The  quench  depth  (T  -T,)  is 
roughly  the  same  for  each  run  shown.   In  both  films  the  [1,0,0]  peaks  are  much 
narrower  than  the  [1,-1,0]  peaks,  due  to  the  anisotropy  of  the  domains.   These 
results  support,  although  not  as  strongly,  the  isotropic  growth  found  in  the  bulk 
sample. 

Comparing  the  results  from  the  different  samples,  it  is  clear  that  domain 
growth  proceeds  much  faster  in  the  stoichiometric  film.   Fitted  growth  exponents,  a 
in  equation  4.3,  are  listed  in  Table  3.   The  quoted  errors  are  dominated  by  estimates 
of  the  systematic  error  arising  from  uncertainty  in  the  resolution  and  finite  size 
corrections.   The  stoichiometric  film  has  a  a  0.4,  roughly  consistent  (within 
experimental  uncertainties)  with  curvature  driven  growth  (a  =  1/2).   On  the  other 
hand,  values  of  a  «  0.2  are  obtained  for  the  film  with  extra  Cu.   This  cannot  be 
reconciled  with  curvature  driven  growth. 


Table  3.  Quenches  and  fitted  exponents. 


film     peak    Tf(°C)   Tj  (°C)    T^  (°C)  ^rjchannels*)      a  [Eq.  4.3]       m[Eq.  2.12] 


A 

(1,-1,0) 

365 

393 

386 

17.1 

0.42±0.08 

4.9 

A 

1,0,0) 

364 

396 

386 

14.8 

0.41±0.11 

4.2 

B 

1,-1,0) 

278 

320 

299 

13.5 

0.23±0.04 

2.2 

B 

(1,0,0) 

283 

320 

299 

19.0 

G.18±0.06 

1.9 

*  — 4   Q-l 

At  the  (1,-1,0)  peak  position,  1  channel  corresponds  to  dQ=8.51>«10     A     . 
At  the  (1,0,0)  peak  position,  1  channel  corresponds  to  AQ=8. 71x10     A     . 
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FIG.  26.  Typical  [1,-1,0]  profiles  from  film  B  quenched  to  278°C.  For  clarity 
of  presentation,  the  data  points  are  averaged  in  groups  of  ten,  and 
spectra  are  offset  in  increments  of  500  counts.   Counting  times  are 
1000  sec.  for  each  1024  channel  scan.   Only  the  central  500  channels 


are  shown.   Each  channel  corresponds  to  AQ=8. 51x10 


The  solid 


lines  are  gaussian  fits.   Times  after  the  quench:  (a)  1.50x10    sec, 

(b)  3.48x10^  sec,  (c)  6.48x10"^  sec,  (d)  S.eSxlo"^  sec,  (e)  1.70x10^  sec 


89 


2.8 


o<C 


2.4 


CD 
O 


2.0 


1.6 


loq 


FIG.  27.  Log,p,L  (2T/(r  in  A)  vs.  Log(t)  (time  after  quench  in  seconds)  for 

film  A  (upper  curves)  and  film  B  (lower  curves).   The  solid  lines  are 
power  law  fits  to  the  late—  time  data.  The  discontinuity  in  the  film  B 
[1,0,0]  data  is  a  known  experimental  artifact. 
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Figure  28  is  a  plot  of  L  vs  log(t)  for  the  same  four  runs.   Over  the  range  of 
times  measured,  it  is  apparent  that  the  data  for  stoichiometric  film  are  not 
consistent  with  L  «  log(t),  while  the  Cu  rich  film  data  are  consistent  with  a 
logarithmic  growth  law.   Since  it  is  energetically  favorable  for  ordered  domains  to 
have  CugAu  stoichiometry  the  extra  Cu  in  film  B  can  be  viewed  as  a  diffusing 
impurity. 

Recent  simulations^^^'^^^)  of  diffusing  impurities  in  an  Ising  model  with 
Kawasaki  dynamics  found  that  domain  coarsening  begins  with  an  apparent  t^/^ 
regime  which  saturates  at  an  impurity  concentration  dependent  value.   At  late 
times  impurities  concentrate  along  the  domain  walls,  eventually  restraining  waU 
motion.   The  order  of  the  ground  state  degeneracy  and  the  symmetry  of  the  order 
parameter^      ''  involved  do  not  change  these  results.. 

Although  the  Cuq^qAuq  21  data  do  not  display  an  obvious  early  time  t^/^ 
regime,  the  general  behavior  is  consistent  with  that  observed  in  the  simulations.  In 
particular,  the  ongoing  concentration  of  extra  Cu  atoms  along  domain  walls  is  a 
physicaUy  attractive  and  reasonable  scenario.  (^^)  Further,  large,  massive  clusters 
of  Cu  atoms  wiU  have  less  mobility  than  a  single  atom,  hence,  at  late  times  the 
situation  may  resemble  a  model  with  quenched  impurities. 

It  has  been  predicted^^^)  that  an  Ising  system  with  quenched,  random 
impurities  wiU  show  domain  coarsening  of  the  form  of  equation  2.12.   Figure  29 
shows  a  plot  of  the  log(L)  vs  log[log(t)].   The  straight  line  fits  give  the  exponents  m 
listed  in  table  3.  It  is  interesting  to  note  that  the  best  fit  value  m  «  2  for  the 
non-stoichiometric  film  is  essentially  equal  to  the  first  order  e  expansion  result 
(reference  80)  m  =  1.82  obtained  for  a  three  dimensional  Ising  model  with 
impurities.   Unfortunately,  the  available  range  of  data  is  probably  too  small  to 
extract  a  definitive  value  for  m.   It  is  evident  from  figure  26  that  the  film  B  data 
could  also  be  consistent  with  m  =  1. 
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FIG.  28.  L  vs.  Log,r^(t)  for  the  runs  summarized  in  table  3.   The  solid  lines  are 
the  best  linear  fits  to  the  data  for  each  run. 
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FIG.  29.  Log,Q(L)  vs.  Log,Q[Log,Q(t)]  for  the  runs  summarized  in  table  3.  The 

solid  lines  are  linear  least  square  fits  to  the  last  half  of  the  data 
for  each  run. 
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It  might  be  concluded  from  the  data  on  the  stoichiometric  sputtered  film 
that  there  is  no  change  due  to  sample  thickness  down  to  7300X.   This  would  not  be 
a  valid  conclusion.  The  time  dependence  of  the  domain  size  in  the  stoichiometric 
film  was  found  to  be  L=t  '        '  ,  consistent  with  L=t  '  .   However,  the  correction 
used  in  association  with  the  sputtered  films  was  a  probable  upper  bound  to  the  true 
resolution.   This  tended  to  make  the  Log,p|(L)  vs.  Log,p,(time)  plot  into  a  straighter 
line.   In  analyzing  the  data  on  the  4500A  M.B.E.  film  it  was  discovered  that 
significant  over  correction  of  the  data  was  necessary  to  force  the  Log,p,(L)  vs. 
Login(time)  plot  to  be  linear.   Moreover,  the  slope  of  such  an  over  corrected  line 
was  0.4,  a  value  close  to  that  obtained  for  the  sputtered  film.   Fortunately,  the 
off— stoichiometric  film  was  thicker  (10,000A)  than  its  counter  part,  and  having 
analyzed  the  two  sputtered  films  in  the  same  way,  means  that  the  difference  must 
be  due  to  the  composition  not  thickness.   These  ambiguities  due  to  the  correction 
serve  to  stress  the  importance  of  measuring  the  resolution  very  carefully. 

The  M.B.E.   Films 

Experimental  Details 

The  M.B.E.  films  were  grown  on  a  sapphire  (ALO2)  substrate  with  a  buffer 
layer  of  Niobium  at  the  University  of  Illinois  by  J.  A.  Dura  and  C.  P.  Flynn. 
Niobium  (Nb)  is  bcc  with  a  lattice  constant  of  3.3oX.   As  with  the  sputtered  films 
the  [1,1,1]  direction  is  normal  to  the  film  surface.   The  M.B.E.  films  are  single 
crystals  ,  therefore,  the  [1,0,-1]  is  in  the  plane  of  the  film  but  the  [0,1,0]  is  not.   The 
orientation  of  the  crystal  was  determined  from  the  [1,1,1],  [1,0,-1]  (or  the  [2,0,-2]), 
and  the  [0,2,0].   The  lattice  constants  for  the  three  films  varied  less  than  0.3  percent 
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between  different  film  thicknesses  and  less  than  0.5  percent  between  different 
directions  in  each  film.   The  in  plane  lattice  constant  is  3.770±.006X,  which  is  the 
same  value  found  for  the  bulk  sample.  The  out  of  plane  lattice  constant  is 
3.753±.006A.   This  difference  in  the  lattice  constant  strongly  suggests  that  strain  is 
important  in  these  systems. 

Shortly  after  the  films  were  prepared  Rutherford  Back  Scattering  was  used  to 
determine  film  thickness  and  sample  stoichiometry.   The  atomic  percentages  of  Au 
in  the  4500A,  TloX,  and  26oX  films  were  26.4±0.4,  27.2±0.4,  and  26.0±1.0 
respectively.  The  results  on  the  4500±100^  film  indicated  that  the  interface 
between  the  Cu^Au  layer  and  the  Nb  buffer  layer  was  sharp.   The  corresponding 
interfaces  for  the  710±20X  and  260^  film  overlapped.  After  the  time  dependent 
X— ray  data  had  been  collected  Auger  electron  spectroscopy  was  used  to  test  for 
impurities  in  the  4500a  and  260a  films  and  Rutherford  Back  Scattering  was  again 
used  for  the  710A  sample.  The  results  indicated  that  the  4500X  and  71oA  films  were 
free  of  impurities  and  that  the  260  a  film  had  slight  traces  of  oxygen  (2±1  atomic 
percent)  and  niobium  (3±1  atomic  percent)  diffused  throughout  the  film.   For  the 
thinnest  film  oxygen  was  also  present  in  the  niobium  buffer  layer  with  a  peak  in 
concentration  at  the  Nb/Cu^Au  interface  which  could  be  explained  by  oxygen 
diffusing  through  to  the  buffer  layer  along  grain  boundaries.   The  niobium 
concentration  decreases  to  3  atomic  percent  in  approximately  70R  for  the  260X 
sample.   For  the  4500a  film  the  niobium  concentration  dropped  to  zero  in 
approximately  170a.  For  all  the  films  there  was  of  course  surface  carbon  and 
oxygen. 

The  data  taken  on  the  4500a  film  was  treated  similarly  to  that  of  the  bulk 
single  crystal.   Radial  scans  were  collected  with  the  exact  same  experimental  set  up 
used  to  collect  the  bulk  radial  data.   The  only  change  for  the  26oA  and  710A  films 
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was  in  the  heating  element.   Due  to  the  lower  thermal  conduction  in  the  thinner 
films  it  was  necessary  to  improve  the  heating  element.   To  increase  the  temperature 
stability  a  larger  thermal  mass  was  used.   The  larger  mass  made  quenching  more 
difficult.   Therefore,  a  cooling  line  was  incorporated  into  the  heating  element.   The 
details  of  the  heating  element/sample  mount  are  shown  in  figure  30. 

The  use  of  the  wound  wire  resistors  as  shown  in  figure  30  significantly 
improved  the  furnace.   Using  two  resistors,  one  above  and  one  below  the  sample 
minimized  temperature  gradients.   Also  the  old  heating  element  was  prone  to 
burn  out,  where  its  life  time  was  dependent  on  the  vacuum  maintained  in  the 
furnace.   The  resistors  also  had  a  tendency  to  burn  out.   However,  the  life  time  of 
the  resistors  is  easily  increased  by  annealing  them  before  placing  them  in  the  copper 
mount.   The  resistors  were  heated  by  running  DC  current  through  them,  in  air,  till 
the  ceramic  casting  just  started  to  melt.  This  procedure  took  less  than  a  minute  per 
resistor.   Resistors  treated  in  this  way  have,  so  far,  not  burnt  out. 

Improved  quench  rates  were  achieved  by  flushing  nitrogen  through  the 
cooling  lines  in  the  mount.   Further  improvement  in  cooling  rates  can  be  achieved 
by  using  helium.   However,  the  film  temperature  lagged  that  of  the  mount  making 
the  thermal  conductivity  of  the  sample  the  limiting  factor  in  the  quench.   The  new 
sample  mount  could  be  quenched  45°C  in  10  seconds  using  nitrogen  gas.   Due  to  the 
lag  in  temperature  the  films  were  quenched  4  to  45  C  in  10  to  25  seconds  depending 
on  depth.   The  quenches  for  the  4500A  film  using  the  old  heating  element  (shown  in 
figure  9)  took  45  to  90  seconds.   The  sample  temperature  was  measured  by  a  type 
K-thermocouple  glued  to  the  back  of  the  ALO2  substrate.   A  second  thermocouple 
was  placed  on  the  copper  mount  were  the  lower  resistor  is  clamped.   This  second 
thermocouple  was  used  as  a  reference  for  temperature  control.   Quenches  where 
performed  in  a  similar  way  as  before,  only  now  a  digital  delay  generator  was  used  to 
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cooling   line 


FIG.  30.  Heating  element/sample  mount  used  for  the  260a  and 

710^  M.B.E.  grown  fils  of  Cu^Au.   The  Cooling  lines  were 

used  to  quench  the  sample.   The  placement  of  the  heating 
resistors  above  and  below  the  sample  position  minimizes 
temperature  gradients. 
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control  the  flow  of  nitrogen  and  start  the  scans.  The  Micristar  temperature 
controller  was  used  in  the  same  manner  as  described  in  connection  with  the  bulk 
sample.   Since  this  data  was  collected  the  entire  procedure  has  been  automated. 

Analysis  of  the  M.B.E.   Film  Data 

The  shape  of  the  [1,0,-1]  peak  for  the  4500X  film  is  LSQ,  therefore,  the 
convolution  fitting  routine  was  used,  as  it  was  with  the  bulk  sample.   A  scan  of  the 
[1,1,1]  of  Silicon  was  used  to  represent  instrumental  resolution.   This  scan  was 
convolved  with  a  gaussian  representing  finite  size  and  strain.   The  shape  of  the 
[1,0,-1]  peak  for  both  the  71oX  and  260X  films  is  gaussian.   Hence,  the  scans  were 
fit  by  straight  gaussians  and  then  the  effective  resolution  was  subtracted  in 
quadrature  as  was  done  with  the  sputtered  films. 

Only  one  fundamental  peak  in  the  [1,0,-1]  direction,  namely  the  [2,0,-2],  is 
observable  with  CuK*^  radiation.   For  the  bulk  data  on  the  [0,1,0]  two  fundamental 
peaks  were  used  to  separate  finite  size  from  strain  effects.   Recalling  that,  in  units  of 
inverse  distance  (X     )  finite  size  should  be  a  constant  and  strain  should  go  as  the 
magnitude  of  the  peak,  the  true  correction  at  the  [1,0,— 1]  must  be  between  one  half 
and  one  times  the  width  of  the  [2,0,-2],  where  instrumental  resolution  is  subtracted 
and  added  back  in  separately.   Although  finite  size  effects  for  peaks  out  of  the  plane 
of  the  film  are  important  the  crystal  size  in  the  plane  would  only  have  to  be  three 
times  larger  than  film  thickness  in  order  to  make  the  contribution  from  finite  size 
small  compared  to  that  of  strain.   Furthermore,  from  the  difference  in  measured 
lattice  constants  relative  to  the  film  surface  strain  should  be  substantial.   The  large 
contribution  to  the  width  from  strain  has  also  been  reported'^       '  to  exist  in  one  of 
the  M.B.E.  Cu„Au  films  prepared  at  the  University  of  Illinois.   Therefore,  the 
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widths  were  corrected  by  assuming  finite  size  and  strain  contribution  was  all  due  to 
strain.   This  correction  is  a  lower  bound  to  the  true  contibution.   An  upper  bound 
would  be  calculated  by  assuming  finite  size  broadening  dominates  over  strain.  This 
upper  bound  is  roughly  the  width  of  the  [2,0,-2]  as  listed  in  table  4.   Table  4  also 
lists  the  corrections  used  for  the  three  films. 

The  other  work^^^^^  on  a  2000X  M.B.E.  Cu^Au  film  investigated  surface 
ordering  with  time  resolved  grazing  incidence  x-ray  experiments.  In  these  grazing 
incidence  experiments  small  disordered  clusters  were  found  to  persist  near  the 
surface  as  the  bulk  ordered.   Domain  growth  was  seen  to  be  depth  independent. 
However,  the  average  order  parameter  increases  faster  in  the  bulk  than  at  the 
surface. 


Table  4 
M.B.E.  Film  Corrections 

M.B.E. 
film 

1  Instrumental 
Resolution 
|at  [2,0,-2; 

]  Width 

of 
![2,0,-2] 

]  Width  of 
[2,0,-2] 

!   -IR 

[Instrumental 
Resolution 
'at  [1,0,-1 

iFinitesize 
and  strain 
'correction 

jTotal 
Cor- 
'rection 

260X 

1       2.70 

j      17.9 

1      15.0 

1       3.17 

1        7.50 

[9.38 

710X 

1       2.70 

1       8.95 

1       7.28 

1       3.17 

j       3.64 

|5.32 

4500R 

1       2.09 

!       6.52 

]       5.26 

1       2.45 

1       2.63 

[3.49 

all  widths  (xlO"*^  X"^) 


Results  from  the  M.B.E.  films 


In  the  following  pages  are  plots  of  L  vs.  time,  Log.p(L)  vs.  Log,p(time),  and 

L  vs.  Log,Q(time)  for  all  three  films.  In  each  case  the  growth  is  slower  than 

1/2 
Cahn— Allen  curvature  driven  growth  (t  '   ).   These  plots  also  show  a  temperature 

dependence  for  the  thinner  films.  In  the  three  plots  for  the  4500 A  film,  shown  in 
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figures  31  through  33,  only  one  temperature  has  been  shown,  in  part  because  all  the 
quenches  display  similar  behavior  and  in  part  for  clarity.   The  plot  of  L  vs.  time  for 
the  4500a  film  is  similar  to  the  equivalent  plot  for  the  bulk  samples.   However,  the 
plot  of  Log,Q(L)  vs.  log,Q(time)  for  this  film  is  different  from  the  bulk.   At  times 
shortly  after  the  quench  the  domains  grow  with  at'    time  dependence.   The  power 
of  the  time  dependence  then  continuously  decreases  so  that  at  the  end  of  the 
experiment  the  power  of  the  time  dependence  was  0.2.  It  is  difficult  to  determine 

from  the  data  if  0.2  is  a  limiting  value  or  not.   However  the  range  of  data  that 

0  2 
follows  a  power  law  of  t  '    is  also  consistent  with  a  logarithmic  time  dependence. 

In  a  plot  of  L  vs.  Log,Q(time),  (figure  33),  one  can  see  that  at  late  times  the  data  is 

almost  a  straight  line.  In  fact  if  a  slightly  larger  resolution  correction  is  used  the 

slight  curvature  at  late  times  disappears  and  the  late  time  data  clearly  indicate  a 

logarithmic  time  dependence  for  L.   Recalling  that  the  minimum  value  for  the  finite 

1/2 
size  and  strain  correction  was  used  it  is  evident  that  L  grows  as  t  '     at  early  times, 

then  crosses  over  to  a  logarithmic  or  lower  power  law  at  late  times  in  the  4500^ 

film.   Monte  Carlo  simulations  have  shown  a  slowing  of  the  growth  due  to  finite  size 

effects,^     *  but  other  effects  such  as  vacancies^      '',  strain  and  impurities  (or 

incorrect  stoichiometry)  induced  pinning  are  more  likely  to  be  the  dominate  cause 

for  the  slower  growth. 

A  crossover  from  t  '    to  t  '   ,  as  discussed  in  chapter  2,  has  been  predicted 

for  late  times. *>     ^  The  late  time  power  of  0.2  is  probably  not  evidence  for  a  t  ' 

1/4 
crossover.   The  t  '    late  time  power  law  is  for  pure  systems,  here  the  small 

deviation  from  critical  stoichiometry,  among  other  things,  could  easily  explain  the 

slower  growth  at  late  times.  The  data  on  the  710R  and  260X  films,  discussed  below, 

do  not  favor  at'    late  time  power  law  suggesting  that  the  4500^  film  data, 

although  consistent  with  this  prediction  is  not  related  to  it. 
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FIG.  31.  Size  of  a  domain  L=2r/r  vs.  time  for  the  450oA 
M.B.E.  film  of  CugAu. 
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FIG.  32.  Log,^[L=2T/r]  vs.  Log,p,[time]  for  the  450oR 


'10 
M.B.E.  film  of  CUjAu. 
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FIG.  33.  Size  of  a  domain  L=2T/r  vs.  Log,Q[time] 
for  the  4500A  M.B.E.  film  of  CugAu. 
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The  most  interesting  M.B.E.  data  is  connected  with  the  710A  film.  A  plot  of 
L  vs.  Log,^(time)  is  shown  in  figure  34.   For  the  bulk  data  these  curves,  as  a 
function  of  temperature,  do  not  cross  each  other.   However,  for  the  710X  film  it 
seems  that  at  higher  temperatures  the  late  time  growth  is  slower.  In  fact,  it  appears 
as  if  growth  is  completely  stopping  as  the  transition  temperature  is  approached. 
Figure  35  is  a  plot  of  Log,r|(L)  vs.  Log,p(time)  for  the  same  data  shown  in  figure 
34.   The  Log— Log  plot  shows  a  slowing  of  the  growth  at  late  times  similar  to  that 
found  in  the  4500a  film.   However,  here,  the  late  time  power  law  for  shallow 
quenches  continuously  decreases  all  the  way  to  a  power  of  zero,  not  to  0.2  as  in  the 
bulk.   In  the  plot  of  L  vs.  Log,^(time)  (figure  36),  the  deeper  the  quench  the 
straighter  the  curve  suggesting  slower  growth  with  increasing  temperature.  The 
growth  is  also  slower  at  later  times  for  all  temperatures,  as  seen  from  the  curves 
being  concave  down.  It  appears  that  a  logarithmic  behavior  is  approached,  or  that 
the  period  over  which  such  a  behavior  is  obeyed  is  increased,  with  decreasing 
temperature.   As  discussed  in  chapter  two  logarithmic  time  dependence  is  expected 
in  systems  with  energy  barriers  to  growth  that  are  a  function  of  the  size  of  the 
domains  L.   Moreover,  these  barriers  are  assumed  to  decay  or  be  over  come 
exponentially.   In  the  data  from  the  710X  film  it  is  likely  that  there  are  barriers  to 
growth  relative  to  the  bulk.   In  fact,  for  shallow  quenches  these  barriers  seem  almost 
stable  at  late  time  giving  rise  to  the  near  freezing  of  the  growth.   As  the 
temperature  is  lowered  there  is  a  greater  driving  force  and  the  barriers  become  less 
stable  and  begin  to  decay.   However,  there  is  a  temperature  dependence  on  the 
decay  rate.   A  recent  Monte  Carlo  simulation  considers  ordering  in  a 
two— dimensional  lattice  with  annealed  vacancies  and  a  NCOP>      '  This  two 
dimensional  model  produces  a  crossover  from  Cahn— Allen  to  logarithmic  like 
growth,  even  for  half  a  percent  concentration  of  vacancies.   The  authors  also 
mention  a  "slab  effect"  in  which  the  domains  become  pinned. 
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FIG.  34.  Size  of  a  domain  L=2t/T  vs.  time  for  the  710X 
M.B.E.  film  of  Cu-^Au.  Five  different  temperatures 

are  shown. 


105 


710AM.B.E.  Cu.Au 


2.7 

^2.5 
ti2.3 


to 
o 


2.1  - 


1.9 


— 1 — 

T 1 1 1 1 1 1 1 

T 1 1 1 1 1 1 1 

- 

Oi 

W^&^^ 

_ 

J^ 

^ 

- 

' '  '^ 

T=407.5'C 

- 

fa  1 

- 

o  403.5°C 

~ 

- 

/a  / 

A  400.l''C 

- 

■ 

1 

°  395 .9°C 

~ 

1 

l>%5,tQ 

1 

1 

364.3°C 

1 

/ 

— — ^  - 

1 

3         4         5 
Logio[time  (sec)] 


FIG.  35.  Log^Q[L=2T/r]  vs.  Log^Q[time]  for  the  710X 

M.B.E.  film  of  Cu^Au.   Early  time  slopes  are  about 
0.4.   Late  time  slopes  are  between  0.0  and  0.2. 
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FIG.  36.  Size  of  a  domain  L=2t/T  vs.  Log,p,[time] 


for  the  710X  M.B.E.  film  of  CugAu. 


107 


The  effect  of  these  barriers  is  clearly  seen  in  the  260X  film.   Figure  37  is  the 
L  vs.  time  plot  for  the  two  relatively  deep  quenches  done.   At  late  times  there  is 
little  change  in  the  size  of  the  domains.   Figures  38  and  39  show  Log,p,(L)  and  L  as 
functions  of  Log,  Q(time)  respectively,  in  both  cases  at  late  times  there  is  a  freezing 
of  the  growth. 

This  freezing  of  the  growth  in  the  260°  film  and  in  the  710a  film  at 
temperatures  close  to  T  ,  might  be  due  to  domain  walls  becoming  pinned  at  the 
surfaces  or  may  be  connected  with  the  strain  in  the  crystal.^     '      '  The  strain  is 
largest  in  the  260a  film  and  is  a  function  of  temperature.  Changes  in  strain  might 
explain  the  710a  films  temperature  dependence. 

By  using  the  minimum  resolution  correction  the  slowing  down  is  emphasized. 
The  slowing  down  is  still  present  if  a  larger  correction  is  used,  however,  the  absolute 
freezing  is  not  as  apparent.   The  raw  uncorrected  widths  for  the  quenches  on  the 
260a  film  are  still  slightly  changing  at  the  end  of  the  experiment. 

Scaling  and  System  Dimension 

From  equation  2.20  it  can  be  seen  that  the  scaling  law  for  S(Q)  can  indicate 
the  dimensionality  of  the  system.  To  determine  the  effective  dimension  for  growth 
in  the  260X  and  710A  films  one  would  plot  S(Q)-  (t'^  vs.  Q/cr  for  d=2  and  3  for 
several  times,  then  the  plot  that  showed  scaling  would  indicate  the  dimension.   A 
consideration  of  the  instrumental  resolution  shows  that  an  uncorrected  a  should  be 
used.   However,  there  is  still  one  effect  of  instrumental  resolution  that  must  be 
considered,  whether  or  not  the  instrumental  resolution  integrates  over  one  or  more 
dimensions.   In  the  following  it  has  been  assumed  that  the  instrumental  resolution 
normal  to  the  scattering  plane  integrates  over  that  dimension.   Estimates  for  the 
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FIG.  37.  Size  of  a  domain  L=2r/r  vs.  time  for  the  260R 


M.B.E.  filmof  CugAu. 
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FIG.  38.  Log^Q[L=2T/r]  vs.  Log^g[time]  for  the  260.^ 
M.B.E.  Clm  of  Cu^Au.   Early  time  slope  is  0.35. 
Late  time  slopes  are  slightly  above  zero. 
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FIG.  39.  Size  of  a  domain  L=2t/T  vs.  Log.p,[time] 


for  the  260X  M.B.E.  film  of  Cu^Au. 
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260X  film  indicate  that  at  least  at  late  times  the  assumption  is  valid.  The  width 
due  to  strain  being  smaller  for  the  710a  film  strengthens  the  assumption  for  that 

system. 

2 
Figures  40  and  41  show  S(Q)-(7-    and  S(Q)-cr  vs.  Q/a  respectively  for  a 

quench  to  T  — 21.5°C  on  the  710X  film.   Figure  40  clearly  indicate  three 

dimensional  behavior  relative  to  the  2  dimensional  scaling  attempt  shown  in  figure 

41.   The  corresponding  plots  for  a  quench  to  T  -31.5°C  on  the  26oA  film  are  in 

figures  42  and  43.   For  comparison  the  raw  unsealed  scans  are  shown  in  figure  44. 

An  initial  observation  suggests  that  the  growth  in  the  260a  film  is  2  dimensional. 

However,  there  are  several  factors  that  make  the  difference  unclear.   First,  the 

spread  of  the  scaled  plots  at  the  peak  position  is  smaller  for  the  3  dimension  scaling. 

Second,  if  only  times  after  half  an  hour  are  considered  the  3  dimensional  plot  looks 

considerably  better  which  might  be  due  to  the  assumption  about  integrating  over 

one  dimension  only  being  true  at  later  times.   Three,  any  error  in  the  fitted  width 

would  more  important  in  the  3  dimensional  case. 

Figures  45  and  46  show  weak  scaling  (dividing  counts  by  fitted  amplitude)  as 

well  as  data  scaled  by  the  corrected  widths.   These  plots  might  incorrectly  be 

interpreted  as  indicating  two  dimensional  growth  in  the  thinnest  film.   However  the 

difference  probably  arises  from  instrumental  resolution  and  strain  (see  table  4). 

From  the  weak  scaling  plots  it  can  be  concluded  that  the  background  subtraction 

that  was  made  in  all  the  scaling  plots  does  not  contribute  difficulty  in  deciding  on 

the  dimension  for  the  26oX  film.   The  fact  that  the  scaling  in  the  710a  film  is 

clearly  3  dimensional  compared  to  being  unclear  in  the  260a  film  suggests  that 

there  is  a  difference  in  the  dimensionality  of  the  scaling  between  the  two. 

Additional  support  for  the  two  dimensional  case  is  found  in  L(t)  which  for  the  260a 

film,  unlike  in  the  others,  is  larger  than  the  film  thickness  at  late 
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FIG.  40.  3D  scaled  counts  vs.  scaled  position  (x~Xo)/'''mpasiirpd 

(counts -(T^)  for  S(Q)  of  the  710^  M.B.E.  film  at  a  few 
selected  times. 
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FIG.  41.  2D  scaled  counts  vs.  scaled  position  (x~Xo)/'''measured 

(counts-  (t)  for  scans  of  S(Q)  on  the  710X  M.B.E.  film  for 
a  few  selected  times..   Plot  shows  2D  scaling  does  not  work. 
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times  (see  figure  37).   It  is  reasonable  to  assume  that  at  early  times  the  scaling  in 
the  260R  is  three  dimensional  and  that  at  some  late  time  it  crosses  over  to  two  a 
dimensional  system.  This  may  account  for  some  of  the  confusion  in  the  scaling. 
Unfortunately  the  larger  effective  resolution  lessens  the  change  in  the  uncorrected 
widths  so  that  at  late  times  the  raw  widths  do  not  change  nearly  as  much  as  at 
early  times. 

A  Brief  Look  at   the  Early  Time  Behavair  of  the  Bulk  Single  Crystal 

Figure  47  shows  some  raw  scans  of  the  bulk  sample  of  Cu«Au  taken  at 
National  Synchrotron  Light  Source  (NSLS)  at  Brookhaven  National  Laboratory 
(BNL).   The  scans  were  collected  using  a  silicon  diode  array  PSD  made  by 
Princeton  Applied  Research  (PAR).   The  experimental  set— up  was  similar  to  that 
used  for  the  radial  data  collected  at  The  University  of  Florida  as  depicted  in  figure 
10,  only  the  synchrotron  replaces  the  rotating  anode  as  the  x— ray  source.  The 
instrumental  resolution  associated  with  this  early  time  data  is  much  less  than  the 
width  of  the  peaks  for  all  of  the  data  shown  here.   The  quench  took  15  seconds  using 
the  automated  quenching  described  in  connection  with  the  M.B.E.  films.   One 
second  scans  were  taken  for  500  seconds  starting  at  the  beginning  of  the  quench.   It 
was  expected  that  these  early  time  scans  would  be  gaussian,  however,  they  are  LSQ. 
The  LSQ  shape  may  be  due  to  the  sample  having  been  ordered  for  a  week  before 
going  to  BNL.   Although  the  history  dependence  of  the  line  shape  crossover  was 
already  known,  it  was  not  known  that  the  effect  could  be  eliminated  to  this  extent. 

A  plot  of  Amplitude  vs.  time  is  shown  in  figure  48.   At  about  125  seconds  the 
amplitude  rapidly  increases  in  a  linear  fashion.   A  plot  of  the  Log,  j^( width)  vs. 
Log,p(time),  shown  in  figure  49  displays  a  sharp  change  in  slope  from  nucleation  to 
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growth  also  at  125  seconds.   Figure  50  is  a  plot  of  the  integrated  intensity  (amp-T). 
It  is  some  what  surprising  that  this  curve  is  concave  down  after  the  onset  of  critical 
nucleation.   From  similar  plots  showing  the  late  time  data  (figure  18)  it  is  known 
that  the  integrated  intensity  saturates  at  a  much  later  time,  indicating  that  the  up 
turn  at  the  end  of  the  plot  is  real.   This  effect  is  an  indication  that  two  different 
types  of  nucleation  are  likely  to  be  involved. 


The  integrated  intensity  is  expected  to  follow  the  rule^  ' 


1  =  1 


1— exp 


-[(t-to)/r] 


m 


(4.5) 


where  I    is  the  maximum  obtainable  intensity,  r  is  a  rate  constant,  t    represents  a 

nucleation  time,  and  m  is  an  integer  that  depends  on  the  type  of  nucleation 

involved. 

4       ho  mogeneous 
m=    3        impurites 
2        s  u  r  face 

If  two  types  of  nucleation  are  involved  then  it  may  be  possible  that  the  initial 
increase  in  intensity  is  due  to  surface  (or  impurity  site)  nucleation  with  a  short 
nucleation  period  t  .  This  early  nucleation  would  soon  be  swamped  by 
homogeneous  nucleation.  This  scenario  would  lead  to  an  intensity  corresponding  to 
the  sum  of  two  equations  of  the  form  of  4.5,  each  with  different  values  for  I  ,  t^, 
and  T. 

Another  unexpected  discovery  made  at  NSLS  is  the  existence  of  satellites  at 
[0,0.9,0]  and  [0,1.1,0].  These  satellites  are  shown  in  figure  51.   They  are  found  in  the 
radial  direction,  but  not  in  the  transverse  direction,  suggesting  that  they  are 
connected  with  type  2  domain  walls.   Similar  satellites  have  been  found  in  the  CuAu 
system  which  also  has  an  order-disorder  transition^       ''  The  satellites  in  Cu^Au 
behave  similar  to  those  in  CuAu.  The  satellites  appear  and  then  grow  in  intensity 
after  a  quench  as  the  main  peak  grows,  but  as  the  integrated  intensity  of  the  main 
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peak  (the  [1,0,0])  saturates  the  satellites  decay.   The  satellites  do  not  appear  in 
quenches  to  20°C  or  more  below  T  .  The  satellites  in  CuAu  result  from  domain 
walls  forming  periodically  every  5  lattice  constants,  creating  a  super  domain  of 
length  10  lattice  constantsA^^^^   The  satellites  in  Cu^Au  may  be  due  to  the  same 
phenomena,  arising  from  a  AuCu  type  structure  forming  in  type  2  walls  or  at  the 
surface. ^^^^'^^^^   Satellites  have  also  been  observed  in  alloys  of  composition  near 
Cu  Au^^'^'*'^^^'^^^^  which  are  believed  to  arise  from  a  modulation  of  the  density 
inside  each  domain.   Furthermore,  the  satellites  seen  in  off-stoichiometric  samples 
are  found  in  a  direction  which  would  associate  them  with  type  1  walls.   The 
existence  of  satellites  in  off-stoichiometric  Cu^Au  suggests  that  it  is  the 
offstoichiometric  nature  of  type  2  walls  that  produces  our  satellites.   Our  satellites 
association  with  type  2  walls  has  interesting  implications  to  the  early  time  behavior 
for  example,  they  strengthen  the  assertion  that  scaling  is  not  valid  between  early 
and  late  times.   A  more  detailed  analysis  of  this  synchrotron  data  is  planed  for  the 
near  future  and  is  mentioned  here  in  part  as  an  indication  of  future  directions,  and 
in  part  for  completeness. 


CHAPTER  5 
CONCLUSIONS 


Current   Results 


The  development  of  order  in  the  bulk  single  crystal  Cu,,  Au  sample  after  a 
thermal  quench,  from  the  disordered  state  exhibits  three  distinct  regimes, 
nucleation,  ordering,  and  coarsening.   The  first  is  characterized  by  a  gaussian 
structure  factor  associated  with  short— range  order.   Afthe  beginning  of  the  ordering 
regime  the  shape  is  gaussian  but  quickly  crosses  over  becoming  rapidly  more 
lorenzti an— squared  with  time.   At  the  onset  of  coarsening  the  structure  factor  is 
very  well  described  by  a  lorenztian— squared  function.   During  coarsening,  domain 
growth  is  consistent  with  a  curvature— driven  growth  law  for  both  type— 1  and 
type— 2  walls.   Moreover,  the  power  "a"  associated  with  this  growth  has  been  found 
to  be  0.5  to  a  high  degree  of  accuracy  in  the  radial  direction  and  to  a  lesser  degree  in 
the  transverse  direction,  implying  that  scaling  during  coarsening  is  valid  in  an 
anisotropic  system.   However,  the  change  in  Line  shape  clearly  shows  that  S(q,t)  can 
only  be  scaled  during  coarsening  and  not  during  ordering.   This  fact  is  easily  seen  by 
viewing  the  disordered  phase  as  a  thick  domain  wall,  in  which  case  the  size  of  the 
wall  relative  to  the  ordered  domain  size  would  be  different  for  different  times  so  that 
scaling  would  not  exactly  recreate  the  state  of  the  system.   Our  data  is  qualitatively 
described  by  the  theory  recently  put  forth  by  Z.  Lai  based  on  a  three  component 
order  parameter.   This  theory  produces  a  slightly  modified  lorenztian-squared 
structure  factor  at  late  times  and  a  gaussian  shape  at  early  times.   Although  a 
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change  in  the  line  shape  exists  in  Lai's  theory  the  history  dependence  of  this  change 
is  not  explained.   However,  the  level  of  short— range  order  and  its  effect  on 
nucleation  is  likely  to  play  an  essential  role  in  explaining  the  history  dependence.   It 
is  possible  that  a  previously  well  ordered  sample  which  has  been  disordered  close  to 
the  transition  temperature  (within  a30  C)  has  a  more  homogeneous  nucleation.  If 
nucleation  takes  place  throughout  the  sample  with  a  high  density  of  nucleation  sites, 
then  the  system  would  reach  coarsening  very  rapidly  and  have  a  short  ordering 
regime. 

Ordering  in  an  off— stoichiometric  sample  is  slower  than  in  a  stoichiometric 
one.  The  extra  atoms  in  an  off— stoichiometric  sample  behave  like  diffusive 
impurities,  giving  rise  to  a  logarithmic  like  time  dependence. 

It  was  hoped  that  the  data  on  the  M.B.E.  films  would  show  a  clear  crossover 
to  two  dimensional  growth.   However,  it  is  unclear  whether  an  effective  two 
dimensional  stage  was  reached  or  not.   The  physics  involved  in  the  M.B.E.  films  is 
probably  dominated  by  the  strain  which  gives  rise  to  a  slower  type  of  growth  at  late 
times.  This  slower  growth  rate  in  the  M.B.E.  films  has  a  temperature  dependence 
such  that  at  late  times  the  power  of  the  time  dependence  is  lower  for  higher 
temperatures.   In  fact,  the  power  of  the  time  dependence  at  late  times  becomes 
almost  zero  for  the  710^  film  as  the  temperature  approaches  T  .   This  temperature 
dependence  could  be  related  to  strain  in  the  film  being  greater  at  higher 
temperatures.  Another  trend  is  that  the  thinner  the  film  the  lower  the  power  of  the 
time  dependence  for  a  given  temperature,  but  this  too  is  probably  due  to  greater 
strain  for  thinner  films. 
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Future  Directions 


As  discussed  at  the  end  of  chapter  6,  the  history  dependence  of  the  line  shape 
crossover  is  not  understood.   Its  likely  connection  with  the  nucleation  process  is  an 
interesting  problem  that  may  broaden  our  understanding  of  how  short  range  order 
effects  the  transition.   In  addition,  there  seems  to  be  a  void  of  information  regarding 
the  ordering  stage.   Nucleation  in  general  has  been  studied  in  many  systems,  as  well 
as  domain  coarsening,  yet  it  is  the  ordering  stage  that  does  not  obey  scaling  laws. 
Further  investigation  into  this  intermediate  stage  would  be  useful.   An  investigation 
of  the  satellites  found  in  the  radial  direction  about  the  [0,1,0]  may  provide  some 
insights.   The  effects  of  stoichiometry  in  a  bulk  single  crystal  on  both  types  of  walls 
could  also  have  some  useful  information. 


APPENDIX  A 
THE  EFFECT  OF  TYPE  1  WALLS  ON  SCATTERING 


The  width  of  a  superlattice  peak  is  of  course  greatly  affected  by  the  existence 
of  domain  walls.   Assuming  the  half  diagonal  glide  domain  boundaries  perpendicular 
to  each  of  the  axes  are  completely  independent  from  each  other  the  structure  factor 
for  the  m^m^m^  ordered  cell  is, ^     '' 

F(m,m2m„)  =  F  exp[i^(m,m2m-,)]  (A.l) 

where  F  is  from  3.25. 

F  =  f^^+f(-,^{exp[Ti(h+k)]+  exp[Tz(k+l)]+exp[Tz(l+h)]} 
The  phase  ^(m^m^m^)  accounts  for  cells  in  different  ground  states.   This  phase 
factor  has  four  possible  values,  [0,  r(h+k),   r(k+l),  T(l+h)],  related  to  the  four 
displacements  [0,  (a, +a2)/2,  (a2+a«)/2,  (a„+a,)/2].   For  the  structure  factor  F  the 
intensity  is  given  F  F.   After  summing  over  all  cells  one  finds, 

I=F   S        exp[i^(m,m2m„)]exp[i(A;-fc  )-(m,a,+m2a2+m-,a„)] 

^^2^3  (A.2) 

•F    S        exp[-i^(m:m2mx)]exp[— 2(fc-fc  )-(m:a,+m2a2+mAa^)] 
m^m^m^ 

Rewriting  {k-k  )  as  [(h+h^)b^+(k+h2)b2+(l+h2)b2]  where  [h^b^+h2b2+h2b2] 

is  the  departure  of  the  tip  of  the  diffraction  vector  from  the  point  [h,k,l].   The 

intensity  is  now  written. 
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I=F*F    S 


S        exp{z^[(m,m„in„)-(m'mAmx)]} 
m^m2m3  m^m^m^  (A.3) 

•ex-p{2ti[(m^-m^)h^+{m^-m^)h^+{m^-m^)h^]}. 

Defining  nj=mj-mr  and  ^(njn2n3)  =  i^{m^m^m^)-(l){m^m^m0  one  of  the  sums 

may  be  eliminated  by  combining  terms  with  the  same  value  of  n,  n^n^,  along  with 

introducing  an  average  phase  factor  and  multiplying  by  the  number  of  terms, 

N(nj^n2n3). 


-,2 


I(hjh2h3)  =  F^  S    S    S  N(n^n2n3)<exp(i^(n^n2n3) 


> 


^1  "2  "3 


exp[2Ti(njhj+n2h2+n3h3)  (A.4) 

The  three  sets  of  perpendicular  anti-phase  boundaries  are  independent, 
hence,  <exp[i^(n^n2n3)]>  =  <exp[i>(nj)]>  <exp[z^(n2)]>  <exp[z^(n3)]>.   The 
phase  difference  between  distant  cells  is  the  sum  of  the  phases  between  all  the  cells 
on  a  line  connecting  them.   Mathematically, 

^(n)  =  ^(l^)+^(2-l)+...  +  ^[n-{n-l)]  (A.5) 

Therefore,  <exp[z>(n)]>  =  <exp[i^(l)]>  1^1.   There  are  only  two  possible  values 
for  ^(1)  along  the  a^  direction,  (ii(l)  =  T(k+1)  with  probability  a  and  ^(1)  =  0  with 
probability  1-a.   The  probability  a  is  that  of  crossing  an  anti-phase  domain 
boundary  in  a  distance  |  a| .  Therefore,  along  the  a,  direction, 
<exp[i^(l)]>  =  (l-a)+a-exp[Ti(k-|-l)] 

=  l-2a  sinV(k+l)]  (A.6) 

Summing  over  n^  the  phase  factor  takes  on  two  values,  1  if  k+1  is  even  and 
(l-2a) '  °1 1  if  k+1  is  odd.   The  same  logic  applies  to  n2  with  1+h  and  n^  with  h+k. 
For  fundamental  reflections  the  sums  (h+k),  (k+1),  and,  (1+h)  are  even  and 
no  effect  is  produced.   For  supedattice  peaks,  two  of  the  indices  [h,k,l]  are  of  the 
same  parity  and  the  other  one  is  different.    For  convenience  the  convention  that  all 
superlattice  peaks  be  written  such  that  h  and  1  have  the  same  parity  will  be 
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adopted.   The  intensity  with  this  convention  simplifies  to, 

I(hjh2h3)  =  f2  S    S    S  N(n^n2n3)(l-2a)l°ll  (l-2a)l°3l 


°1  °2  °3 


exp[2Ti(njh^+n2h2+n3h3)].  (A.7) 

For  very  smaU  njn2n3  ,  N(njn2n3)  %  ^^^2^2'  ^^^  ^slcIots  (l-2a)l°ll  and 
(l-2a)'°3l  make  only  small  values  of  |nj  and  \n^\  important.   Consequently 
^(°1°2°3)  can  be  approximated  by  N^(N2-|  n2 1  jNg.^^^^  Equation  A.7  is  now 
expressed  as  three  independent  products. 

I(hjh2h3)  =  f2[n  S  (l-2a)l°ll  exp(2T2n,hJ 


^ 


S  (N2-I  n2 1 )  exp(2Tzn2h2)  (A.8) 


°2 


N3E  (l-2a)l'^3'  exp(2T2n3h3) 


^3 


Now  consider  the  sum  over  n2.  This  sum  is  equivalent  to  a  product  of  sums. 
+N2 

E    (N2-|n2|)exp(2Tzn2h2) 
°2="^2 

+N— 1  +N-1 

S  exp(2Tzm„h„)        S  exp(-2TzmAh„)  (A. 9) 

m2=0  ^  ^   m^=0  ^  ^ 

The  sums  over  m2  and  m^  have  a  geometric  progression  of  the  form, 

a  +  ar+ar2  +  ...  +1='-^  (A.IO) 

Summing  over  m2  and  mx  in  this  way  yields, 

exp(2TtN2h2)-l       exp(-2TiN2h2)-l 

exp(2Tzh2)  -  1       exp(-2Tzh2)  -  1 

=  2-2cos(2TN2h2)       sin2(TN2h2)  (A.ll) 

2-2cos(2Th2)  sin^(Th2) 
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Returning  to  the  sums  over  n,  and  n„,  both  of  which  are  of  the  form, 
£  (l-2a)l^lexp(-2mh) 


n=— <B 

CD 


=  1  +    S  [(l-2a)exp(2Tzh)]°  +    E  [(l-2o)exp(-2Tzh)]'^ 
n=l  n=l 

(l-2a)exp(2Tih)  (l-2a)exp(-2Tih) 

=  1  + + 


l-(l-2a)exp(2Tzh)      l-(l-2a)exp(-2Tzh) 

a(l-Q) 

=  -5 2 .  (A.12) 

a'^+(l-2a)sin''(Th) 

Making  the  final  approximation  that  h,,  h^,  and  h„  are  small  the  intensity  is  given 

2 

.,     N,a         sin  (rN^h^)        N„a 

I(h  h  h  )  =  F^-^-t 2 H T^ 2-        (^-^3) 

The  intensity  for  a  superlattice  peak  as  given  above  is  roughly  a  thin  disk  lying  in 
the  h,  h«  plane  with  a  thickness  determined  by  the  crystal  dimension  N^a^.   Type  2 
walls  would,  of  course  thicken  the  disk.   Figure  8  shows  the  orientations  of  theses 
disks  for  one  plane  of  reciprocal  space.   Also  this  equation  shows  that  the  radius  of 
the  disk  increases  with  the  probability  a  of  anti-phase  boundaries. 

Normally  the  condition  ^h^h^h^)  equal  (l/2)I(h,h2h„)  would  be  used  to 
define  the  width  of  the  peak  in  terms  of  a  (a  a  L).   However,  the  approximation  of 
h^hoh^  being  small  lead  to  problems.   Therefore,  one  starts  with  a  more  general 
equation  for  the  width. '^     '' 

/   I    dt 

FWHM«^5_| (A.14) 

max 

Although  the  integrated  intensity  is  greatly  affected  by  the  amount  of  order  vs. 

disorder  it  is  not  affected  by  the  existence  of  sharp  domain  walls.   Hence,  the 

2 
integrated  intensity  must  be  F  N,N2N2.   Considenng  I(h,  h2h2)  as  defined  in 
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equation  A. 8  with  the  same  approximations  for  n^  as  before  we  have, 
I(hjh2h3)  =  F^[Nj  S    (l-2a)l°ll  exp(2rzn^h^) 


n^=0 


.  N3  S    (l-2a)l^3l  exp(2Tzn3h3)  (A.  15) 


'03=0 


sin  (TN2h2)  -. 


The  intensity  is  a  maximum  at  h,=h„=h„=0.  So  that  A. 15  becomes, 

I(h^h2h3)  =  f2NjN2N3  S    (l-2a)l^ll    E    (l-2o)l'^3l    (A.16) 


nj=0 


n3=0 


where, 


lim     sin  (>N^h^)^N. 
Expressing  n,  a,  and  n„a«  in  terms  of  a  distance  (s)  along  the  disk, 

Vi  =  2^^h/y  h^+k^i^ 

^3^3 
we  can  write. 


n.a.  =  2Tsl/y  t,2^j^2^^2 


|n^l+|n3|=l^lH      (IH  +  I^I)    = 


a  ,/  .2  ,  ,_2  ,  ,  2 


Converting  the  sums  over  n,  and  n„  into  an  integral  over  (s)  gives, 


I 


max 


=  F^N^N2N3j'(l-2a)^l^lds 


(A.17) 


(5.18) 


(A.19) 


(A.20) 


Substitution  of  A.20  into  A.  14  shows, 
FWHMoc 


^f\^^n^ 


F^N^N2N3, 


-Pln{l-2a) 


/(l-2fl)^l^lds  2 


(5.21) 


For  small  0,  /7i(l— 2fl)  ~  -2a.   The  final  result  with  (a/a)  =  L,  the  size  of  a  domain 
i,,(90) 


FWHM  .  a/J  =  (2t/L1     (  I  >■  I  +  I  I*  I  ) 


/h^k^l^ 


(A.22) 


135 


This  result  is  for  type  1  walls,  however,  type  2  walls  also  produce  broadening  with  a 
2t/L  dependence.   The  dependence,  /?  for  the  [0,1,0]  type  peak  is  one,  and  for  the 
[1,0,T]  ^  is  v5.  For  a  given  functional  form  equation  A.  14  can  be  used  to  relate  a  the 
corresponding  width  to  L.  For  a  gaussian  width  a  we  have, 

'  =  (^/U     'lll-^.l'l'  (A.23) 

where  as  for  a  lorentzian  squared  width  F, 

r  =  (v^/L)_(^i+l_ii)     .  (A.24) 

/h^  +  k^l'^ 
The  reader  should  be  aware  of  the  above  corrections  to  2.21  for  type  1  walls. 

However,  throughout  this  work  it  is  assumed  that  2.21  is  correct  for  both  types  of 

walls. 


APPENDIX  B 
FUNCTIONAL  DEFINITIONS  USED  IN  FITTING 


In  the  following,  x  is  the  position  coordinate,  and  x    is  the  position  of  the 
peak.   Here  and  throughout  the  text  u  is  used  to  represent  a  gaussian  width  and  F  a 
width  for  a  lorentzian  to  any  power.   The  following  definitions  are  such  that  the 
widths  will  approximately  equal  a  HWHM.  The  amplitude  is  "A". 

Gaussian, 


Lorentzian, 


F(x)=A.exp[[(x-x^)/2(r]^ 
F(x)=  - 


Lorentzian  squared, 


F(x)= 


HT  Theory 
2 


(8) 


1+ 


^-\ 


(x-XQ)/r 


(x-XQ)/r 


(B.l) 


(B.2) 


-2 


(B.3) 


lW=!(fAu-fCu)  . A,    .  r.2r    ri3.[l+cos(Tk)-2cos(Tkj)]       (B.4) 


where. 


r.    =(l-i^)/[i+2;2_  2i    cos(2Tk.)] 
im    ^       m''  '■        m        m      ^       r^ 

i2=(l-2;7-27),   x^={\- ^) 


7  is  the  probability  of  finding  a  type  1  wall  between  two  unit  cells. 
Iri  is  the  probability  of  finding  a  type  2  waU  between  two  unit  cells. 
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At  the  [0,1,0]  a  transverse  direction  is  (x,0,0)  and  k,=2-  ,  ^2=1  ,  k„=0 
Equation  B.4  then  reduces  to. 


where, 


Amp{l+C[l-cos(a(x-XQ))]}   ^  (B.5) 


C=[0.5-{rj+7f]l{rj+^) 


are  fit  parameters. 

In  the  radial  direction  k^=0,  ^2       It'  ^3~^* 
This  leads  to; 

where, 

k„=2.4878856-sin(0.413643rad+(x-x^).2.1816616xl0^rad/PSD^hannel) 
X  is  in  PSD  channels,  0.4136  is  the  position  of  the  [0,1,0]  in  radians,  and  2.48. ..is  a./X 
O.J.K.  theory,^'^^) 


T      X 

J 

0 


3 


l{x)=li^      fdRR^J3/2(^-R)  (B.7) 

{i-R[exp(R2)-l]}^/2 


where  J^ /o  is  a  Bessel  function.   In  the  fitting  routine  i=(x-x  )/r  is  approximated 
as  10.47(1-1  )  if  the  absolute  value  of  i  <  10 

Copies  of  the  fitting  and  numerical  integration  routines  used  can  be  found  in 
Numerical  Recipes  in  Fortran  The  Art  of  Scientific   Computing,  Cambridge 
University  Press  (1985). 


APPENDIX  C 

CALIBRATION  OF  BRAUN  OED-50-M 

POSITION  SENSITIVE  METAL  WIRE  DETECTOR 


After  making  all  the  four  wire  and  two  gas  line  connections  (if  they  fit  they 
are  correct),  open  p-10  gas  bottle  and  set  pressure  at  6  bar.  Next,  slowly  turn  valve 
at  fount  of  "gas  control"  to  open  and  adjust  flow  rate  to  0.8  N//hour  (0.8  units  on 
scale).   The  pressure  should  still  read  6  bar  at  both  the  "gas  control"  and  regulator. 
At  this  point  one  should  always  check  for  leaks,  this  is  easily  done  with  soapy  water 
and  a  syringe.   There  are  two  connections  at  the  back  of  the  gas  control  as  well  as 
the  ones  at  the  detector. 

After  any  leaks  have  been  fixed,  the  detector  and  gas  lines  need  to  be  purged. 
To  purge  the  system  simply  open  the  valve  at  the  back  of  the  gas  control  and  leave 
it  open  for  at  least  3  minutes  (I  recommend  5  minutes).   If  an  expensive  gas  is  being 
used  in  place  of  p-10  it  is  possible  to  connect  a  ruffing  pump  to  the  exhaust  at  the 
back  of  the  gas  control  and  evacuate  the  system  before  opening  the  valve  at  the 
front  of  the  gas  control.   If  this  method  is  used  one  evacuates  the  system  then  closes 
the  valve  at  the  back  before  opening  the  valve  at  the  front.   After  the  system  has 
been  back  filled  the  valve  at  the  fount  is  closed  and  the  system  is  evacuated  and 
back  filled  a  second  time  before  turning  the  high  voltage  on. 

After  closing  the  valve  at  the  back  the  power  can  safely  be  turned  on  (x-rays 
must  not  be  impinging  on  the  detector).  If  the  high  voltage  has  not  been  previously 
adjusted  start  with  it  at  zero  and  slowly  increase  it  to  S.OkV.  Now,  with  four 
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0.002iii.  Ni  foils  or  equivalent  attenuators  in  the  bean,  open  the  shutter.   If  rate 
meter  does  not  show  counts,  raise  high  voltage  (HV),  but  do  not  exceed  3.8kV. 
After  impulses  have  been  seen  on  the  rate  meter  connect  the  PSM  energy  output  to 
the  ADC  input.   Now  the  energy  peak  can  be  displayed.  You  should  be  using  1024 
channels  and  have  the  ADC  gain  at  Ik.   Adjust  the  HV  so  that  the  peak  of  the 
energy  spectrum  is  at  channel  300.   I  have  also  been  told  to  set  it  at  channel  600, 
however  I  believe  this  contributed  to  the  dead  spots  on  the  wire.   The  energy  peak  is 
really  set  with  an  oscilloscope,  however.  Innovative  Technology  assumed  that  no  one 
would  have  one  so  they  gave  the  approximate  energy  that  corresponded  to  the 
correct  settings.   It  is  a  good  idea  to  ask  Innovative  Technology  for  directions  on 
how  one  would  use  an  oscilloscope  to  adjust  the  HV. 

Now  connect  the  PSM's  SCA  output  to  ADC's  gain  input  and  switch  ADC 
to  coinc.   Leaving  the  windows  on  the  ADC  open,  lower  upper  level  and  raise  the 
lower  level  of  the  PSM's  windows  till  approximately  20%  of  the  energy  spectrum, 
without  the  escape  peak,  is  left.   The  escape  peak  builds  up  below  the  main  peak. 
In  practice  I  try  to  keep  as  much  of  the  peak  as  possible  but  make  what  is  left  after 
removing  the  escape  peak  symmetric.   To  display  the  position  spectrum  disconnect 
both  BNC  cables,  switch  ADC  from  coinc.  to  anti.  and  connect  the  PSM's  position 
output  to  the  ADC's  input. 

The  wider  the  windows  the  better  the  counting  rate  but  the  linearity  is 
affected.   The  ends  of  the  wire  will  not  count  as  much  as  the  middle  with  wider 
windows.   I  suggest  that  if  the  peak  is  to  be  in  the  center  of  the  detector   to  keep 
the  count  rate  up  and  if  linearity  becomes  a  problem  assume  a  quadratic 
background. 


APPENDIX  D 
SPECTROMETER  ALIGNMENT 


In  order  to  align  the  scope  on  the  Huber  one  must  first  put  the  tip  of  the 
"pin"  at  the  center  of  rotation.   Then  it  is  a  simple  matter  of  adjusting  the  scope  so 
that  the  cross— hairs  are  on  the  tip  of  the  pin.   To  put  the  tip  of  the  pin  at  the  center 
of  rotation  when  the  scope  is  not  properly  adjusted,  first  adjust  all  the  goniometer 
translations  and  rotations  to  zero.   Next  adjust  the  scope  so  that  the  pin  is  some 
where  in  view.   Now,  separately  for  both  translations  rotate  Phi  180  degrees  noting 
the  position  through  the  scope  that  the  pin  is  rotating  about.  In  other  words,  note 
the  middle  of  the  two  extremes.   Put  the  pin  at  this  location  and  again  rotate  Phi  as 
before  and  re-adjust  the  pin.  Continue  to  rotate  and  re— adjust  until  the  tip  of  the 
pin  does  not  move.   The  height  adjustment  is  easer  but  is  done  the  same  way  only 
Chi  is  rotated  instead  of  Phi. 

Having  put  the  pin  at  the  center  of  rotation  the  scope  can  now  be  aligned. 
By  tightening  and  loosening  the  bolts  that  hold  the  mount  for  the  scope  to  the  Chi 
circle  you  can  adjust  the  scope  up  and  down.   After  the  height  is  correct  loosen  the 
bolts  on  the  bottom  of  the  scope's  mount  and  adjust  the  scope  right  and  left.   Of 
course  you  may  need  to  continue  adjusting  the  scope  by  the  above  method  until  one 
adjustment  does  not  affect  the  other. 

It  is  important  to  first  calculate  the  scattering  angle  for  the  monochromator 
and  wavelength  being  used,  and  carefully  mount  the  monochromator  in  its  housing 
at  the  approximate  position.   The  closer  your  first  guess  is,  the  less  time  it  will  take 
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to  line  up  the  monochromator.   Because  it  is  easiest  to  judge  0  or  90  degrees  one 
should  set  the  monochromator  at  zero  and  adjust  the  translations  while  looking 
down  with  a  ruler  across  the  top  of  the  housing,  then  rotate  the  monochromator  to 
the  calculated  angle.   The  Huber  should  also  be  moved  to  the  calculated  value  of  20. 
To  do  this  loosen  the  set  screw  on  the  under  side  of  the  monochromator  housing, 
and  make  sure  the  beam  path  from  the  monochromator  housing  is  connected  to  the 
Huber  so  that  the  angle  of  the  Huber  can  be  read  as  it  is  moved. 

Now  put  the  detector  up  close  to  the  monochromator  housing  with  the 
window  horizontal  in  order  to  cover  as  large  a  solid  angle  as  possible.   Next,  with  all 
the  slits  open,  recursively  adjust  ^^^^  ,  the  translations,  and  arcs  of  the 
monochromator  until  it  is  tweaked  up.   If  you  do  this  inside  the  x-ray  shielding  the 
X— ray  power  should  be  at  a  minimum.   Moreover  this  operation  is  best  done  with 
two  people,  one  to  adjust  the  monochromator  and  one  to  watch  the  rate  meter. 

Now  that  the  monochromator  is  lined  up  it  is  time  to  adjust  the  fine  position 
of  the  Huber.   The  easiest  way  to  determine  the  correct  position  is  to  put  a 
fluorescent  pad  behind  a  pin  that  is  at  the  center  of  rotation  and  with  out  foils  and 
the  lights  out  look  at  the  florence  relative  to  the  pin.   After  a  ruff  adjustment  using 
the  air  pads  to  move  the  Huber,  a  scope  is  mounted  outside  of  the  shielding  with  a 
view  of  the  pin.   The  scope  allows  more  precise  estimates  of  the  position  of  the  beam 
relative  to  the  pin.   Shining  a  flashlight  near  the  pin  will  allow  one  to  see  the  pin 
while  still  seeing  the  x-ray  fluorescence.   After  the  Huber  is  in  place  an  arm  zero  is 
done  to  position  the  detector  arm,  then  slitting.   If  an  analyzer  is  to  be  used  slitting 
should  be  done  with  the  analyzer  in  place. 

To  line  up  the  analyzer  one  once  again  start  by  carefully  estimating  the 
correct  zeros.   It  is  best  to  put  the  detector  up  close  to  the  analyzer  with  the 
window  horizontal  to  again  get  the  widest  solid  angle  as  possible,  as  seen  from  the 
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analyzer.   With  the  detector  at  the  approximate  position,  adjust  the  analyzer  angle, 
then  move  the  detector  back  and  put  its  window  vertical  so  that  the  detector  angle 
can  be  set.   Be  careful  to  line  up  on  k,  and  not  k^-   The  only  adjustment  of  the 
goniometer  holding  the  analyzer  that  is  important  is  the  forward  and  back  arc 
because  this  adjusts  the  height  of  the  beam  at  the  detector.   The  other  adjustments 
should  be  left  alone,  the  adjustments  change  26,  and  it  is  better  to  simple  readjust 
26. 

To  put  a  sample  at  the  center  of  rotation  of  the  Huber  first  translate  the 
scattering  surface  such  that  the  part  of  the  sample  one  wishes  to  scatter  from  is  on 
the  cross— hairs  of  the  scope.   Next  while  looking  at  the  side  of  the  sample  through 
the  scope  adjust  the  relevant  arc.   The  other  arc  can  be  adjusted  as  desired.   The 
translation  that  moves  the  scattering  face  in  a  direction  normal  to  the  scattering 
surface  can  be  adjusted  in  two  ways.   The  first  is  to  alternately  look  at  the  two  sides 
of  the  sample  and  adjust  its  position  until  rotating  from  one  side  to  the  other 
rotates  Phi  180  degrees.  The  second  method  which  will  work  on  irregular  surfaces  is 
to  watch  the  face  of  the  sample  while  rocking  Phi  back  and  forth  and  adjusting  the 
translation  till  a  spot  on  the  sample  stays  at  the  cross— hairs.   If  the  spot  one  finds  is 
not  the  one  desired  adjust  the  other  translation  and  then  re— adjust  the  second. 
After  one  has  gained  some  experience  with  this  method  it  becomes  easier.   The 
direction  of  apparent  movement  relative  to  the  cross-hairs  is  related  to  the 
direction  the  sample  must  be  moved.   However  once  the  position  is  close,  the  change 
in  magnitude  of  the  apparent  displacement  can  be  used  to  zero  in  on  the  center  of 
rotation.   In  practice  the  first  method  is  usually  a  good  start  and  the  second  method 
is  a  good  check  and  or  is  used  to  fine  tune  the  alignment. 


APPENDIX  E 
DETECTOR  DEAD  TIME  CORRECTION 


To  correct  for  the  non-linearity  of  the  detector  at  high  counting  rates  the 
detector  dead  time  must  be  known.   The  detector  dead  time  correction  is  then'^      '', 

^rue"%bserved/(^"^eff  %bserved^  (^-^^ 

where  n,       is  the  corrected  value  of  n  u-^j.  pj  and  t  cr  is  the  effective  detector 

dead  time.  The  detector  dead  time  can  be  measured  with  an  oscilloscope  or  by  a 

multiple— foil  method.  In  the  multiple— foil  method  n,       is  estimated  from  the 

attenuation  factor  of  a  foil  placed  in  the  beam  path  multiplied  be  the  number  of 

attenuated  counts,   t  rr  is  then  calculated  from  n.        and  n  ij-pypfj  by  equation  E.l. 

Here  a  similar  method  was  used  only  the  power  of  the  x— ray  source  was  lowered 

instead  of  placing  a  foil  in  the  beam.   The  value  of  n,       was  estimated  from  the 

count  rate  at  low  power,  n,      ,  the  monitor  counts  at  low  power,  m.      ,  and  the 

monitor  counts  at  high  power,  °^f,koprved' 

^rue=^low-(°'observed/"'low)  (^-2) 

The  effective  dead  time  was  then  calculated  by  the  relation  described  in  equation 

E.l  and  found  to  be  5.562x10      sec/count. 
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